Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,848,983 C→T 5.1% intergenic (‑169/‑72) USA300HOU_RS14590 ← / → moeZ_3 Lactonase drp35/putative adenylyltransferase/sulfurtransferase MoeZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,848,9830CT 5.1% 472.0 / 7.0 147intergenic (‑169/‑72)USA300HOU_RS14590/moeZ_3Lactonase drp35/putative adenylyltransferase/sulfurtransferase MoeZ
Reads supporting (aligned to +/- strand):  ref base C (73/65);  new base T (2/6);  total (75/72)
Fisher's exact test for biased strand distribution p-value = 1.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 149 total aligned reads displayed.

ACTATCAATTTTATTTTAAACAAGATAACTATCAAATAAAATGATAACGGTTTCATCTATCTATTTTATCGGTCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTTTATATTATAAATATATGACGATTGATGACCCTGAACAGTTTGCTCAGGA  >  USA300TCH1516_ALE/2848843‑2849122
                                                                                                                                            |                                                                                                                                            
   aTCAATTTTATTTTAAACAAGATAACTATCAAATAAAATGATAACGGTTTCATCTATCTATTTTATCGGTCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACaaa                                                                                                                                           >  1:386061/1‑141 (MQ=255)
    tCAATTTTATTTTAAACAAGATAACTATCAAATAAAATGATAACGGTTTCATCTATCTATTTTATCGGTCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAat                                                                                                                                          <  3:15318/141‑1 (MQ=255)
        ttttattttaAACAAGATAACTATCAAATAAAATGATAACGGTTTCATCTATCTATTTTATCGGTCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATgg                                                                                                                                      >  2:83864/1‑140 (MQ=255)
              tttaAACAAGATAACTATCAAATAAAATGATAACGGTTTCATCTATCTATTTTATCGGTCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTg                                                                                                                                >  4:260782/1‑141 (MQ=255)
                 aaaCAAGATAACTATCAAATAAAATGATAACGGTTTCATCTATCTATTTTATCGGTCTAGTGGCTGATTTCAAGCTAGAAATATTGACTGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACg                                                                                                                             >  2:673484/1‑141 (MQ=255)
                   aCAAGATAACTATCAAATAAAATGATAACGGTTTCATCTATCTATTTTATCGGTCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCt                                                                                                                           <  1:226263/141‑1 (MQ=255)
                      aGATAACTATCAAATAAAATGATAACGGTTTCATCTATCTATTTTATCGGACTAGAGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGGTTTa                                                                                                                        >  2:950717/1‑141 (MQ=255)
                      aGATAACTATCAAATAAAATGATAACGATTTCATCTATCTATTTTATCGGTCTAGTGGCTGATATCAAGCTAGCAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATTCGTATTGaaaaaaaa                                                                                                                        <  2:450848/141‑8 (MQ=255)
                       gATAACTATCAAATAAAATGATAACGGTTTCATCTATCTATTTTATCGGTCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTAt                                                                                                                       >  1:195988/1‑141 (MQ=255)
                            cTATCAAATAAAATGATAACGGTTTCATCTATCTATTTTATCGGTCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTaa                                                                                                                 >  1:655171/1‑141 (MQ=255)
                               tCAAATAAAATGATAACGGTTTCATCTATCTATTTTATCGGTCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTaaaaa                                                                                                              <  1:85938/141‑1 (MQ=255)
                                   aTAAAATGATAACGGTTTCATCTATCTATTTTATCGGTCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTaa                                                                                                          >  1:357443/1‑141 (MQ=255)
                                   aTAAAATGATAACGGTTTCATCTATCTATTTTATCGGTCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTaa                                                                                                          >  1:358134/1‑141 (MQ=255)
                                     aaaaTGATAACGTTTTCATCTATCTATTTTATCGGTCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACa                                                                                                        <  2:178078/141‑1 (MQ=255)
                                     aaaaTGATAACGGTTTCATCTATCTATTTTATCGGTCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGATAAATATGCGTATTGACGCTTTATTTTAAAAATTAACa                                                                                                        <  3:190613/141‑1 (MQ=255)
                                      aaaTGATAACGGTTTCATCTATCTATTTTATCGGTCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACAt                                                                                                       <  4:117068/141‑1 (MQ=255)
                                       aaTGATAACGGTTTCATCTATCTATTTTATCGGTCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATg                                                                                                      >  4:424708/1‑141 (MQ=255)
                                          gATAACGTTTTCATCTATCTATTTTATCGGTCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCtt                                                                                                   <  2:392658/141‑1 (MQ=255)
                                          gATAACGGTTTCATCTATCTATTTTATCGGTCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGCCAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCtt                                                                                                   <  3:215277/141‑1 (MQ=255)
                                          gATAACGGTTTCATCTATCTATTTTATCGGTCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGATAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCtt                                                                                                   <  2:993712/141‑1 (MQ=255)
                                          gATAACGGTTTCATCTATCTATTTTATCGGTCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCtt                                                                                                   <  1:854739/141‑1 (MQ=255)
                                           aTAACGGTTTCATCTATCTATTTTATCGGTCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACCTGCTta                                                                                                  <  4:351331/141‑1 (MQ=255)
                                             aaCGGTTTCATCTATCTATTTTATCGGTCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTtata                                                                                                >  2:221358/1‑141 (MQ=255)
                                             aaCGGTTTCATCTATCTATTTTATCGGTCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTtata                                                                                                >  4:63725/1‑141 (MQ=255)
                                               cGGTTTCATCTATCTATTTTATCGGTCTAGTGGCTGATTTCAAGGTGGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAAc                                                                                              <  4:175323/141‑1 (MQ=255)
                                                  tttCATCTATCTATTTTATCGGTCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATg                                                                                           <  2:843412/141‑1 (MQ=255)
                                                  tttCATCTATCTATTTTATCGGTCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATg                                                                                           <  1:487220/141‑1 (MQ=255)
                                                  tttCATCTATCTATTTTATCGGTCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATg                                                                                           <  1:965997/141‑1 (MQ=255)
                                                  tttCATCTATCTATTTTATCGGTCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATg                                                                                           >  4:323229/1‑141 (MQ=255)
                                                  tttCATCTATCTATTTTATCGGTCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATg                                                                                           <  2:757899/141‑1 (MQ=255)
                                                    tCATCTATCTATTTTATCGGTCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGtt                                                                                         <  2:478322/141‑1 (MQ=255)
                                                      atctatctATTTTATCGGTCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTta                                                                                       >  2:282830/1‑141 (MQ=255)
                                                          atctatTTTATCGGTCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGa                                                                                   <  1:747361/141‑1 (MQ=255)
                                                            ctatTTTATCGGTCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAg                                                                                 <  2:1006124/141‑1 (MQ=255)
                                                            ctatTTTATCGGTCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAg                                                                                 <  1:194170/141‑1 (MQ=255)
                                                               ttttATCGGTCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGgag                                                                              >  2:62226/1‑141 (MQ=255)
                                                                   aTCGGTCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTa                                                                          >  1:249726/1‑141 (MQ=255)
                                                                    tCGGTCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTaa                                                                         >  1:556577/1‑141 (MQ=255)
                                                                    tCGGTCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTaa                                                                         >  1:164433/1‑141 (MQ=255)
                                                                      ggTCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAAcc                                                                       >  3:371394/1‑141 (MQ=255)
                                                                      ggTCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAAcc                                                                       >  2:117366/1‑141 (MQ=255)
                                                                       ggcTATTGGCAGACTTCAAGCTAGGAATATGTAATGACAATGCAACTCTGTTAAATTGATGGACGTAGACAAACATGCGGATGGACGCTTTATTTTAAAAATTAACATGATTATGACATGTTTATAGAAAGAGATTAACCt                                                                      >  2:413377/3‑141 (MQ=255)
                                                                        tCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGATAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTa                                                                     >  2:736795/1‑141 (MQ=255)
                                                                        tCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTa                                                                     >  2:114637/1‑141 (MQ=255)
                                                                          tAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATg                                                                   <  1:318011/141‑1 (MQ=255)
                                                                           aGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGa                                                                  <  2:1075602/141‑1 (MQ=255)
                                                                                cTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTTTTAAAATGATGGACGTAGACAAATATGCGTATCGACGCTTTATTGTAAAAATTAACAAGCTTATAACATTTTTATATAAGGAGATAAACCTATGAACTaa                                                             >  4:209174/1‑140 (MQ=255)
                                                                                cTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTAt                                                             >  1:395330/1‑141 (MQ=255)
                                                                                cTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTAt                                                             <  4:332371/141‑1 (MQ=255)
                                                                                cTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTAt                                                             >  2:757746/1‑141 (MQ=255)
                                                                                  gATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCa                                                           <  4:166740/141‑1 (MQ=255)
                                                                                     ttCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGt                                                        >  2:669867/1‑141 (MQ=255)
                                                                                     ttCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGt                                                        <  3:429430/141‑1 (MQ=255)
                                                                                     ttCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGt                                                        >  2:823647/1‑141 (MQ=255)
                                                                                     ttCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGACATTAACCTATGAACTATCAAGt                                                        >  2:823632/1‑141 (MQ=255)
                                                                                         aGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCtt                                                    >  1:705525/1‑141 (MQ=255)
                                                                                          gCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGATAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGga                                                                               >  3:63189/1‑113 (MQ=255)
                                                                                             aGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTttat                                                >  1:761535/1‑141 (MQ=255)
                                                                                              gAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTttat                                                >  4:278734/1‑141 (MQ=255)
                                                                                              gAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTttata                                               >  2:523053/1‑141 (MQ=255)
                                                                                              gAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTttata                                               >  1:878645/1‑141 (MQ=255)
                                                                                               aaaTATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTttatat                                              <  1:378078/141‑1 (MQ=255)
                                                                                                 atatTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTttatatta                                            >  2:325236/1‑141 (MQ=255)
                                                                                                      gAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTTTATATTATAAat                                       <  1:282830/141‑1 (MQ=255)
                                                                                                         tGACAATACAACTCTGTTAAAATGATGGACGTAGATAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTTTATATTATAAatata                                    <  4:449971/141‑1 (MQ=255)
                                                                                                          gACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTTTATATTATAAatatat                                   <  2:557068/141‑1 (MQ=255)
                                                                                                          gACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTTTATATTATAAatatat                                   <  4:66185/141‑1 (MQ=255)
                                                                                                              aTACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTTTATATTATAAATATATGACg                               <  2:164433/141‑1 (MQ=255)
                                                                                                              aTACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTTTATATTATAAATATATGACg                               <  2:694429/141‑1 (MQ=255)
                                                                                                                tcAATTCTGTGAAAATTTTGTTTGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTTTATATTATAAATATTTGACTAt                             >  4:190734/2‑141 (MQ=255)
                                                                                                                 cAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTTTATATTATAAATATATGACGAtt                            >  2:532794/1‑141 (MQ=255)
                                                                                                                 cAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTTTATATTATAAATATATGACGAtt                            >  2:449659/1‑141 (MQ=255)
                                                                                                                  aaCTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTTTATATTATAAATATATGACGATtg                           >  2:882562/1‑141 (MQ=255)
                                                                                                                   aCTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTTTATATTATAAATATATGACGATtga                          <  3:131234/141‑1 (MQ=255)
                                                                                                                    ctctGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTTTATATTATAAATATATGACGATtgat                         <  2:607956/141‑1 (MQ=255)
                                                                                                                      ctGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTTTATATTATAAATATATGACGATtgatga                       <  1:366244/141‑1 (MQ=255)
                                                                                                                        gTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTTTATATTATAAATATATGACGATTGATGAcc                     >  2:132768/1‑141 (MQ=255)
                                                                                                                         ttAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTTTATATTATAAATATATGACGATTGATGAccc                    >  3:50301/1‑141 (MQ=255)
                                                                                                                          tAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTTTATATTATAAATATATGACGATTGATGACCCt                   <  2:989388/141‑1 (MQ=255)
                                                                                                                           aaaaTGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTTTATATTATAAATATATGACGATTGATGACCCTg                  <  3:431033/141‑1 (MQ=255)
                                                                                                                           aaaaTGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTTTATATTATAAATATATGACGATTGATGACCCTg                  >  1:681552/1‑141 (MQ=255)
                                                                                                                           aaaaTGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTTTATATTATAAATATATGACGATTGATGACCCTg                  >  1:76870/1‑141 (MQ=255)
                                                                                                                            aaaTGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTTTATATTATAAATATATGACGATTGATGACCCTGa                 >  4:410118/1‑141 (MQ=255)
                                                                                                                              atgatgGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTTTATATTATAAATATATGACGATTGATGACCCTGAAc               >  3:339153/1‑141 (MQ=255)
                                                                                                                              atgatgGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTTTATATTATAAATATATGACGATTGATGACCCTGAAc               >  3:274457/1‑141 (MQ=255)
                                                                                                                               tgatgGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTTTATATTATAAATATATGACGATTGATGACCCTGAACa              <  1:926800/141‑1 (MQ=255)
                                                                                                                                gatgGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTTTATATTATAAATATATGACGATTGATGACCCTGAACAg             >  2:456236/1‑141 (MQ=255)
                                                                                                                                 atgGGCGGAGATAAGTAGGCTTTTTTATGCTTTTTTTTAATATTTAACATGGTTATAACATGTTTTTAGGAGTATATTTACCTATGAAGTGTCAAGTTCTTGTATGTTATAAATATATGACGATTGATTGTCCGGAACAGt            <  4:246468/141‑1 (MQ=255)
                                                                                                                                 atgGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTTTATATTATAAATATATGACGATTGATGACCCTGAACAGt            <  2:882266/141‑1 (MQ=255)
                                                                                                                                 atgGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTTTATATTATAAATATATGACGATTGATGACCCTGAACAGt            <  4:452168/141‑1 (MQ=255)
                                                                                                                                 atgGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTTTATATTATAAATATATGACGATTGATGACCCTGAACAGt            <  1:1061431/141‑1 (MQ=255)
                                                                                                                                  ttgACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTTTATATTATAAATATATGACGATTGATGACCCTGAACAGtt           <  3:152384/139‑1 (MQ=255)
                                                                                                                                  tgGACGTAGACAAATCTGCGTATTGACGATTTAATAAAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTTTATATTATAAATATATGACGATTGATGACCCTGAACAGtt           <  2:232985/141‑1 (MQ=255)
                                                                                                                                     aCGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTTTATATTATAAATATATGACGATTGATGACCCTGAACAGTTTgg        <  1:1032180/141‑2 (MQ=255)
                                                                                                                                     aCGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTTTATATTATAAATATATGACGATTGATGACCCTGAACAGTTTGc        >  3:76758/1‑141 (MQ=255)
                                                                                                                                        tAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTTTATATTATAAATATATGACGATTGATGACCCTGAACAGTTTGCTCa     <  2:933916/141‑1 (MQ=255)
                                                                                                                                        tAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTTTATATTATAAATATATGACGATTGATGACCCTGAACAGTTTGCTCa     >  1:52836/1‑141 (MQ=255)
                                                                                                                                         aGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTTTATATTATAAATATATGACGATTGATGACCCTGAACAGTTTGCTCAg    <  3:273445/141‑1 (MQ=255)
                                                                                                                                         aGACAAATATGCGTATTGACGCTTTATTTCAAAAATTAACATGCTTATAACATGTTGATAGAAGTAGATTAACCTATGAACTATCAAGTTCTTTTATATTATAAATATATGACGATTGATGACCCTGAACAGTTTGCTCAg    <  3:119872/141‑1 (MQ=255)
                                                                                                                                           aCAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTTTATATTATAAATATATGACGATTGATGACCCTGAACAGTTTGCTCAGGa  >  1:470809/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                            
ACTATCAATTTTATTTTAAACAAGATAACTATCAAATAAAATGATAACGGTTTCATCTATCTATTTTATCGGTCTAGTGGCTGATTTCAAGCTAGAAATATTGAATGACAATACAACTCTGTTAAAATGATGGACGTAGACAAATATGCGTATTGACGCTTTATTTTAAAAATTAACATGCTTATAACATGTTTATAGAAGGAGATTAACCTATGAACTATCAAGTTCTTTTATATTATAAATATATGACGATTGATGACCCTGAACAGTTTGCTCAGGA  >  USA300TCH1516_ALE/2848843‑2849122

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: