Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 512,643 A→G 40.0% 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,6430AG40.0% 81.2 / 138.9 136intergenic (+711/‑6052)recR/speA_2Recombination protein RecR/Arginine decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (40/40);  new base G (24/30);  total (66/70)
Fisher's exact test for biased strand distribution p-value = 5.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 144 total aligned reads displayed.

AGGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCAT  >  USA300TCH1516_ALE/512503‑512781
                                                                                                                                            |                                                                                                                                           
 ggTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTtaa                                                                                                                                            >  2:448187/1‑141 (MQ=255)
 ggTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTtaa                                                                                                                                            <  2:503402/141‑1 (MQ=255)
 ggTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTAAAAAAGTTATTGACTTAAATGTgaa                                                                                                                                            >  2:255666/1‑141 (MQ=255)
  gtgtAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATATTAAAAAGTTATTGACTTACATGTtaac                                                                                                                                           >  2:163213/1‑140 (MQ=255)
   gaagttaagctccttagcgtcgatggtagtcgaacttacgttccgccagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaaataaaatggacaagataaAAAAAGTTCTTGACTTAAATGTtgata                                                                                                                                          <  2:127203‑M1/27‑1 (MQ=255)
    gaagttaagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTtgataa                                                                                                                                         <  2:620734‑M1/28‑1 (MQ=255)
     tAATTTTTACTGTTGTTAAGCAGTTCGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAAATTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATaaa                                                                                                                                        <  1:612252/141‑1 (MQ=255)
     aagttaagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATaaa                                                                                                                                        <  1:143289‑M1/29‑1 (MQ=255)
     aagttaagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggtgaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATaaa                                                                                                                                        >  2:397889‑M1/113‑141 (MQ=255)
      agttaagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATaaaa                                                                                                                                       <  1:135051‑M1/30‑1 (MQ=255)
       aTTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAAt                                                                                                                                      >  1:723122/1‑141 (MQ=255)
         ttttAGTGTGGTTAAGCAGTTTGAAAGACTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGt                                                                                                                                    <  2:586374/141‑1 (MQ=255)
            gctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGtata                                                                                                                                 <  2:318866‑M1/36‑1 (MQ=255)
            gcgccttagtgtcgatagtagtcgaacgtacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGtata                                                                                                                                 >  1:67716‑M1/106‑141 (MQ=255)
              tccttagaggagatggtagacgaacttacgttacgctagagtagaacgttgccaggaagttttttaatcaaattttgtttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTAtaat                                                                                                                               >  2:708134‑M1/104‑141 (MQ=255)
                        gCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGt                                                                                                                     >  1:285962/1‑141 (MQ=255)
                        cgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGt                                                                                                                     <  2:742628‑M1/48‑1 (MQ=255)
                           tggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCgg                                                                                                                  >  1:345610‑M1/91‑141 (MQ=255)
                            ggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGt                                                                                                                 <  1:64525‑M1/52‑1 (MQ=255)
                             gtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTa                                                                                                                <  1:31476‑M1/53‑1 (MQ=255)
                               gAAAGCCTGTATAAGATTTATCTGTAGAGCAAAACAAAACAACTAAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAg                                                                                                              <  1:752188/141‑1 (MQ=255)
                               agtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAg                                                                                                              >  1:288041‑M1/87‑141 (MQ=255)
                               agtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAg                                                                                                              <  1:163352‑M1/55‑1 (MQ=255)
                                gtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTGGGTAAGa                                                                                                             <  2:444856‑M1/56‑1 (MQ=255)
                                           aGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACAt                                                                                                 >  2:79037/1‑141 (MQ=255)
                                                 tatttGTTGAGGCTTACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGTAAAATGAACATTGaaaa                                                                                           >  2:98375/1‑141 (MQ=255)
                                                 tatttGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGaaaa                                                                                           <  2:723122/141‑1 (MQ=255)
                                                  atttGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTATCGGTAAGAAAAATCAACATAGaaaaa                                                                                          >  2:628817/1‑140 (MQ=255)
                                                   tttGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAAATAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTTTTGTCGGTAAGTAAAATGAACATTTAAAACt                                                                                         >  2:616097/1‑141 (MQ=255)
                                                   agagtagaacgttgccaggcagttttttaagcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACt                                                                                         >  1:20167‑M1/67‑141 (MQ=255)
                                                   agagcagaacgttaccaagcagttttttaatcaaaatttagttaaaaaataaaatgcacaagataaAAAAAGTTATTGAATTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACt                                                                                         >  2:388090‑M1/67‑141 (MQ=255)
                                                     tgttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAACTGAACATTGAAAACTGa                                                                                       >  1:116953/1‑141 (MQ=255)
                                                     tgttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGa                                                                                       >  2:379885/1‑141 (MQ=255)
                                                     tgttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAAGTTTAAAAAGTTATTGACTTACATGTTCAGAGATAAAAGAACTAATTCTGGGCTTGAAGAAAAATGAACATTGAAAACTGa                                                                                       >  1:771429/1‑141 (MQ=255)
                                                      gttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAAATAATAAAGAAAATTAATTTTAAAAAGTTATTGACTTAAATGATAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGaa                                                                                      >  2:46049/1‑141 (MQ=255)
                                                        tgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTAAAATGTTAATAAAATGTATAATTAATTATTGTCGGTAAGAAAAATGAACATAGAAAACTGAATg                                                                                    >  2:228186/1‑141 (MQ=255)
                                                           ggCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACa                                                                                 <  1:194033/141‑1 (MQ=255)
                                                           ggCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTATCGGTAAGAAAAATGAACATTGAAAACTGAATGACa                                                                                 >  2:698022/1‑141 (MQ=255)
                                                            gCAAACAAAACAACTAAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACaa                                                                                <  2:150359/141‑1 (MQ=255)
                                                            agttgccaggaagttttttaatcaaattttagttaaaaaataaaatggacaatataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACaa                                                                                >  2:472259‑M1/58‑141 (MQ=255)
                                                              ttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAata                                                                              >  1:176649‑M1/56‑141 (MQ=255)
                                                                gccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATg                                                                            <  2:660262‑M1/88‑1 (MQ=255)
                                                                  caggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGacaatatatc                                                                          <  1:455711‑M1/90‑4 (MQ=255)
                                                                   aggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCa                                                                         <  1:643280‑M1/91‑1 (MQ=255)
                                                                   aggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCa                                                                         <  1:244762‑M1/91‑1 (MQ=255)
                                                                   aaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCa                                                                         <  1:506056/141‑1 (MQ=255)
                                                                     acaacTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAAc                                                                       >  1:582066/1‑141 (MQ=255)
                                                                     acaacTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAAc                                                                       >  2:407714/1‑141 (MQ=255)
                                                                            aactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAtt                                                                <  1:228474/141‑1 (MQ=255)
                                                                             actTAAGAAATAAATTGAATTATTAAAGAAAAAAAAAAAAAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTc                                                               <  2:46261/141‑1 (MQ=37)
                                                                              ctTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTcc                                                              <  1:474382/141‑1 (MQ=255)
                                                                                tAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAAAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaa                                                            >  1:417767/1‑141 (MQ=255)
                                                                                atcagtttttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGACAAATGAACATTGAAACCCGAATGACAATATGTCCACGTTAATTCCaa                                                            >  1:348039‑M1/38‑141 (MQ=255)
                                                                                 aataaaTAACTTGAATTACTAAAAAAAATTAATTTTAAAAAGTTAATGACTTAAATATTAATAAAATGTATAATTAATTCTTGTAGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaa                                                           <  2:645652/138‑1 (MQ=37)
                                                                                   gAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTCATTCCcaaaa                                                         >  2:314660/1‑141 (MQ=255)
                                                                                      ttttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATGAATGATGGTTGTCGGTAGGAAAGATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACg                                                      >  1:775941‑M1/32‑141 (MQ=255)
                                                                                       tttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGt                                                     <  2:261594‑M1/111‑1 (MQ=255)
                                                                                          ggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAAc                                                  >  2:64525‑M1/28‑141 (MQ=255)
                                                                                          cTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAAc                                                  <  1:256678/141‑1 (MQ=255)
                                                                                          cTTGAATTACTAACGAAAATTAATTTTAAAAAGATATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTTGGTACGAAAAATCAACATTGAAAACTAAATGACAATATGTCAACGTTACTTCCAAAAAACGAAAc                                                  >  2:91147/1‑141 (MQ=38)
                                                                                          agtaaaaaaataaataggagaatataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAAc                                                  >  2:309254‑M1/28‑141 (MQ=255)
                                                                                              aaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATCAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACtata                                              >  1:615485‑M1/24‑141 (MQ=255)
                                                                                                  aataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGtt                                          <  2:170352‑M1/122‑1 (MQ=255)
                                                                                                  aCTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAACAAACGTAACTATAAGtt                                          >  2:76769/1‑141 (MQ=255)
                                                                                                  aCTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGtt                                          <  1:600056/141‑1 (MQ=255)
                                                                                                   cTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATCAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATCTGTCAACGTTAATTCCAAAAAACTTAACTATAAGttn                                         >  1:741492/1‑140 (MQ=255)
                                                                                                      aCGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGACAACGTTAATTCGAAAAAACTTAACTATAATTTACaa                                      >  1:107724/1‑141 (MQ=255)
                                                                                                       cGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaaa                                     >  2:611210/1‑141 (MQ=255)
                                                                                                            attaattTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAACTGAGTCAACGATACGAAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACatta                                >  1:774996/1‑141 (MQ=21)
                                                                                                            attaattTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACatta                                >  2:102573/1‑141 (MQ=255)
                                                                                                            attaattTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACatta                                >  1:142603/1‑141 (MQ=255)
                                                                                                               agataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTAttt                             >  1:281918‑M1/7‑141 (MQ=255)
                                                                                                               aattTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTAttt                             >  1:728835/1‑141 (MQ=255)
                                                                                                               aattTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAAAGTTATTTCTAAAAAACGTAACTATAAGTTACAAACATTAttt                             >  2:672914/1‑141 (MQ=255)
                                                                                                                gataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTa                            <  1:75579‑M1/136‑1 (MQ=255)
                                                                                                                 ttttAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAg                           <  2:254896/141‑1 (MQ=255)
                                                                                                                  tttAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGt                          >  1:558655/1‑141 (MQ=255)
                                                                                                                  aaaaaaaaaTTATTAACGTAAATGTTAAAAAAATGTATAATTAATTCATGTCGGTCAGAAGATTGGACATTGCAAACTGAATGACAATATGTCAACATTAAATCCCAACCACTCATTTATAAGTTACAAACATTATTTAGt                          <  2:323002/138‑1 (MQ=16)
                                                                                                                  aaaaaaaaGTTATTTAACAAAATGTTGATAAAATGTATAATTAATGCTTGTAGGTAAGAAAAATAAACATTGAAAACTGAATGACAATATGTCAACATAAAAAAAAAAAAACGTAACTATAAGTTACAAACATTATTTAGt                          <  2:241641/138‑1 (MQ=18)
                                                                                                                   aaAACAAGGTATTGCCTTAATTGTTGATAAAATGTATAATTAATCCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATCACAATATGTCCACGTTAATTCCAAAAAACGTAACTATAACTCACAAACATTCTTTAGTa                         <  2:26426‑M1/139‑1 (MQ=255)
                                                                                                                    tAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAAGAGAAAGGCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTAt                        >  1:533597/1‑141 (MQ=25)
                                                                                                                    tAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACACTATTTAGGAt                        >  1:688484/1‑141 (MQ=255)
                                                                                                                    tAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAAGTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTAt                        <  2:389894/141‑1 (MQ=255)
                                                                                                                      aaaaGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTAttt                      <  1:494369/141‑1 (MQ=255)
                                                                                                                          ggtATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGa                  <  1:526658/139‑1 (MQ=255)
                                                                                                                           ttATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAg                 <  1:427097/141‑1 (MQ=255)
                                                                                                                               taaCTTAAATTTTAATAAAATGTATAATTAATTCTTGTCGGTAATAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTaa             <  2:401686/139‑1 (MQ=255)
                                                                                                                               tGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTaa             <  1:318879/141‑1 (MQ=255)
                                                                                                                               tGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTaa             <  1:318866/141‑1 (MQ=255)
                                                                                                                                gaattAAATGTTAATAAAATGTAGAATTAATTCTTGTCGGTAAGAAAAATGAACATTAAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAAt            <  2:116953/138‑1 (MQ=255)
                                                                                                                                gACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAAt            <  1:379885/141‑1 (MQ=255)
                                                                                                                                 aCTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTTATTATTAGCTAAta           >  2:739836/1‑140 (MQ=255)
                                                                                                                                 aCTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAATTACAAACATTATTTAGTATATATGAGCTAATc           >  2:449923/1‑141 (MQ=255)
                                                                                                                                   ttAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTCATCaa         >  1:594918/1‑141 (MQ=255)
                                                                                                                                   ttAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCaa         >  2:34197/1‑141 (MQ=255)
                                                                                                                                   ttAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCaa         >  1:577116/1‑141 (MQ=255)
                                                                                                                                         gTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcatca   >  2:340693/1‑141 (MQ=255)
                                                                                                                                         gTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcatca   <  1:723197/141‑1 (MQ=255)
                                                                                                                                         gTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCAAAAAAACGTAACTATAAGTTACAAACAATATTTAGTATTTATAAGCTAATCAAAcatca   >  2:291764/1‑141 (MQ=255)
                                                                                                                                          ttAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcatcat  <  1:285112/141‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                           
AGGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCAT  >  USA300TCH1516_ALE/512503‑512781

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