Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 512,643 A→G 46.6% intergenic (+711/‑6052) recR → / → speA_2 Recombination protein RecR/Arginine decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE512,6430AG46.6% 26.6 / 180.0 138intergenic (+711/‑6052)recR/speA_2Recombination protein RecR/Arginine decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (42/31);  new base G (31/33);  total (74/64)
Fisher's exact test for biased strand distribution p-value = 3.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 142 total aligned reads displayed.

GTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCAT  >  USA300TCH1516_ALE/512507‑512781
                                                                                                                                        |                                                                                                                                          
gtAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAGTTAATTTTAGAAAGTTATTGACTTAAATGTtaataa                                                                                                                                        >  1:676298/1‑141 (MQ=255)
gaagttaagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaatttgggttaaaaaataaaatggacaagataaAAAAAGTTATTGCCTTAAATGTtgataa                                                                                                                                        <  1:551241‑M1/28‑1 (MQ=255)
  aaTTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATaaaa                                                                                                                                      >  1:556118/1‑141 (MQ=255)
   gttaagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAAt                                                                                                                                     <  1:361084‑M1/31‑1 (MQ=255)
   aTTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATACAAt                                                                                                                                     >  1:51672/1‑141 (MQ=255)
       agctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGtat                                                                                                                                 >  1:507468‑M1/107‑141 (MQ=255)
          cTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaat                                                                                                                              <  1:174138/141‑1 (MQ=255)
           ccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTAtaatt                                                                                                                             <  2:159283‑M1/39‑1 (MQ=255)
           ccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTAtaatt                                                                                                                             <  2:308029‑M1/39‑1 (MQ=255)
              tgttAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaattaat                                                                                                                          <  1:700825/141‑1 (MQ=255)
              tagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTAtaattaat                                                                                                                          >  1:159283‑M1/100‑141 (MQ=255)
                gcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTc                                                                                                                        >  1:131012‑M1/98‑141 (MQ=255)
                   tcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaaacaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATCTTGATAAAATGTATAATTAATTCTTg                                                                                                                     <  1:468119‑M1/47‑1 (MQ=255)
                     cAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTACTGATGATAGATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTc                                                                                                                   >  1:763288/1‑141 (MQ=255)
                      aGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCg                                                                                                                  >  1:359311/1‑141 (MQ=255)
                        tttGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGt                                                                                                                <  2:534762/141‑1 (MQ=255)
                         gtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGgtc                                                                                                               <  2:576083‑M1/53‑2 (MQ=255)
                         gtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggtaaaaaaataaaatgaacaagataaAAAAAGTTATTGACTTAAAAGTTGATAAAATGTATAATTAATACTTGtcggaa                                                                                                               <  2:276286‑M1/53‑3 (MQ=255)
                             tcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggataacaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGaa                                                                                                           <  2:688678‑M1/57‑1 (MQ=255)
                             ccgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGaa                                                                                                           >  1:253073‑M1/85‑141 (MQ=255)
                               gCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTAGGTAAGaaaa                                                                                                         >  2:116467/1‑141 (MQ=255)
                                    tataGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATTACTTGAATTTCTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAAGGTATAATTAATACTTGTCGGTAAGAAAAATGaa                                                                                                    >  2:272318/1‑141 (MQ=255)
                                     ataGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAAc                                                                                                   <  1:262992/141‑1 (MQ=255)
                                     acgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAAc                                                                                                   <  2:456583‑M1/65‑1 (MQ=255)
                                      tgttccgctaaagtagaacggtgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGAATTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACa                                                                                                  >  2:120985‑M1/76‑141 (MQ=255)
                                         tatttattTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCAATAAAAAAAATGAACATTg                                                                                               <  2:589766/141‑1 (MQ=255)
                                            gctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGaaa                                                                                            >  1:681609‑M1/70‑141 (MQ=255)
                                              ctttATTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGaaaaa                                                                                          >  2:475413/2‑140 (MQ=255)
                                                 tgttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGa                                                                                       >  2:530438/1‑141 (MQ=255)
                                                 tgttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAAGGAACATTGAAAACTGa                                                                                       >  1:738919/1‑141 (MQ=255)
                                                    tgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAAAGATACAAAAGTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATg                                                                                    <  1:215534/141‑1 (MQ=255)
                                                        gCAAACAAAACAACCCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACaa                                                                                <  2:213083/141‑1 (MQ=255)
                                                          ttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAata                                                                              <  2:697037‑M1/86‑1 (MQ=255)
                                                            acaaaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATg                                                                            <  1:63523/141‑1 (MQ=255)
                                                            acaaaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATg                                                                            >  1:609322/1‑141 (MQ=255)
                                                             caaaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATCAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAatacgc                                                                           >  2:76644/1‑138 (MQ=255)
                                                              aaaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTc                                                                          <  2:578456/141‑1 (MQ=255)
                                                               aggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCa                                                                         >  1:445479‑M1/51‑141 (MQ=255)
                                                               aaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCa                                                                         >  1:421206/1‑141 (MQ=255)
                                                                ggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCaa                                                                        >  1:226987‑M1/50‑141 (MQ=255)
                                                                ggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAAAATGTCaa                                                                        <  2:217485‑M1/92‑1 (MQ=255)
                                                                ggcagttttttaatcaaatttcggggaaataatgaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATgtcca                                                                        <  1:374504‑M1/92‑3 (MQ=255)
                                                                 gcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAAc                                                                       <  1:543284‑M1/93‑1 (MQ=255)
                                                                 acaacTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAAc                                                                       >  2:612260/1‑141 (MQ=255)
                                                                 acaacTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAAc                                                                       <  2:707224/141‑1 (MQ=255)
                                                                   agttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGt                                                                     >  1:309610‑M1/47‑141 (MQ=255)
                                                                       ttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAt                                                                 >  1:429984‑M1/43‑141 (MQ=255)
                                                                         ttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTc                                                               <  1:319537‑M1/101‑1 (MQ=255)
                                                                         actTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTc                                                               >  2:184365/1‑141 (MQ=255)
                                                                         actTAAGAAATAACTTAAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTc                                                               <  1:211801/141‑1 (MQ=255)
                                                                          ctTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACTATATGTCAACGTTAATTcc                                                              <  1:462853/141‑1 (MQ=255)
                                                                          ctTAAGAAACAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTcc                                                              <  2:611618/141‑1 (MQ=255)
                                                                              caaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaa                                                          <  1:13513‑M1/106‑1 (MQ=255)
                                                                               gAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaaa                                                         >  2:33253/1‑141 (MQ=255)
                                                                                aattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaaaa                                                        >  1:122046‑M1/34‑141 (MQ=255)
                                                                                    ttggttaaaaaataaaatggacgagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTa                                                    <  2:161795‑M1/112‑1 (MQ=255)
                                                                                    ttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTa                                                    >  1:602935‑M1/30‑141 (MQ=255)
                                                                                     tggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAGAGCTGAATGAAAATATGTCAACGTTAATTCCAAAAAACGTaa                                                   <  2:447480‑M1/113‑1 (MQ=255)
                                                                                      ggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAAc                                                  <  1:773225‑M1/114‑1 (MQ=255)
                                                                                           aaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATaa                                             <  1:37379‑M1/119‑1 (MQ=255)
                                                                                           aaaaataaaataaaaaaaaaaaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATaa                                             >  2:761658‑M1/23‑141 (MQ=255)
                                                                                             aaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAGAACGTAACTATAAGt                                           >  1:575853‑M1/21‑141 (MQ=255)
                                                                                               cTACCGCAAATTCATTTTAAAAAGTTATTGACTTAAATCTGCCTAAAATGGATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGACTGACAATATGTCAACGTTAATTGCAAAAAACGTAACTATAAGTTa                                         >  2:709082/1‑141 (MQ=21)
                                                                                               cTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTa                                         >  1:739369/1‑141 (MQ=255)
                                                                                               atAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTa                                         <  1:410121/140‑1 (MQ=255)
                                                                                                 aaCGAAAATTAATTTTAAAAAGTTATTTACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAAAAAAATGAACATTGAAAACTGAATGACAATATGTCAAAAAAAAAAAAAAAAAACGTAACTATAAGTTACa                                       <  2:522904/141‑1 (MQ=37)
                                                                                                  aCGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaa                                      <  1:78335/141‑1 (MQ=255)
                                                                                                    gAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAAc                                    >  1:631560/1‑141 (MQ=255)
                                                                                                     tggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACa                                   <  1:775966‑M1/129‑1 (MQ=255)
                                                                                                     tggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACa                                   <  2:28402‑M1/129‑1 (MQ=255)
                                                                                                     tggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGCCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACa                                   <  2:126462‑M1/129‑1 (MQ=255)
                                                                                                        attaattTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACatta                                <  1:696416/141‑1 (MQ=255)
                                                                                                        aaaaaataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACatta                                >  2:598591‑M1/10‑141 (MQ=255)
                                                                                                         ttaattTTAAAAAGTTATTGCCTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGACAACGTTAATTCCAAAAAACGTAACTATAGATTACAAACgttat                               >  2:397598/1‑141 (MQ=25)
                                                                                                           agataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTAttt                             <  2:9045‑M1/135‑1 (MQ=255)
                                                                                                           aattTTAAAAAGTTATTGAATTAAATGGTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACGATAAGTTACAAACATTAttt                             <  2:79678/141‑1 (MQ=255)
                                                                                                            attTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATCATTAATGCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTACCTATAAGTTACAAACATTATTTa                            >  2:502645/1‑141 (MQ=255)
                                                                                                             ataaTAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAg                           >  2:264936‑M1/5‑141 (MQ=255)
                                                                                                             aaaaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAAAGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAg                           >  2:649120‑M1/5‑141 (MQ=255)
                                                                                                              ttaAAAAAGTTAGTGACTTAAATGTTAATAAAATTTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGt                          >  1:583339‑M1/4‑141 (MQ=255)
                                                                                                                aAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTGTTTAGTAt                        <  1:568626‑M1/140‑1 (MQ=255)
                                                                                                                aAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTAt                        >  1:698558‑M1/2‑141 (MQ=255)
                                                                                                                aAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATAATTTAGTAt                        <  2:436494‑M1/140‑1 (MQ=255)
                                                                                                                 aaaaaGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTAtt                       >  1:359882/1‑141 (MQ=255)
                                                                                                                    aatttATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTAt                    <  2:327978/138‑1 (MQ=255)
                                                                                                                     aGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACTTTGAAAACTGAATGAACATATGTCAACATTACTTCCAAAAAACGTAACTAGAAGTTACAAACATTATTTAGTATTTATg                   <  2:165383/141‑1 (MQ=25)
                                                                                                                     aGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATg                   <  1:530438/141‑1 (MQ=255)
                                                                                                                      gTTATTGACGTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGa                  >  1:569201/1‑141 (MQ=255)
                                                                                                                       ttattaACTTAAATGATAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTGAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAg                 <  1:547019/141‑1 (MQ=255)
                                                                                                                        tATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGc                <  1:612072/141‑1 (MQ=255)
                                                                                                                           tGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTaa             <  2:173666/141‑1 (MQ=255)
                                                                                                                              cTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCa          >  2:41645/1‑141 (MQ=255)
                                                                                                                              cTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTACAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATAAGCTAATCa          >  2:363307/1‑141 (MQ=255)
                                                                                                                                tAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTCACTATAAGTTACAAACATTATTTAGTATTTATTAGCTAATCaaa        >  2:659053/1‑141 (MQ=255)
                                                                                                                                   aTGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcat     >  1:548570/1‑141 (MQ=255)
                                                                                                                                   aTGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATGCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTACTTATGAGCTAATCAAAcct     >  2:437087/1‑139 (MQ=255)
                                                                                                                                     gTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcatca   >  2:499304/1‑141 (MQ=255)
                                                                                                                                     gTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcatca   <  1:49999/141‑1 (MQ=255)
                                                                                                                                      ttAATAAAATGTATGATTAATTCTTGGCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCACCGTTAATTCGAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcatcga  >  2:621028/1‑139 (MQ=255)
                                                                                                                                      ttAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcatcat  >  2:701272/1‑141 (MQ=255)
                                                                                                                                        |                                                                                                                                          
GTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCAT  >  USA300TCH1516_ALE/512507‑512781

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