Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 512,643 A→G 44.5% 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,6430AG44.5% 62.8 / 155.5 142intergenic (+711/‑6052)recR/speA_2Recombination protein RecR/Arginine decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (43/34);  new base G (22/40);  total (68/74)
Fisher's exact test for biased strand distribution p-value = 2.58e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 144 total aligned reads displayed.

GGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCA  >  USA300TCH1516_ALE/512504‑512780
                                                                                                                                            |                                                                                                                                         
ggTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTta                                                                                                                                           <  1:537993/140‑1 (MQ=255)
      aTTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGAGTTACGTGTTCATaaaa                                                                                                                                     >  2:85883/1‑140 (MQ=255)
      aTTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATaaaa                                                                                                                                     >  2:105277/1‑140 (MQ=255)
      aTTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACACAACTTAAGAATTGACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTCTTGACTTAAATGTtaataaca                                                                                                                                     >  2:604719/1‑138 (MQ=255)
        taagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATg                                                                                                                                   <  2:337708‑M1/32‑1 (MQ=255)
         aagctccttagcgtcgatggtagtcgaacttacggtccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggaacagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGt                                                                                                                                  <  1:757947‑M1/33‑1 (MQ=255)
            aCTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTACATGTTTATAAAATGtaca                                                                                                                               >  1:619455/1‑138 (MQ=255)
            aCTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGtatc                                                                                                                               >  2:670211/1‑139 (MQ=255)
                   gcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATaattaatt                                                                                                                        <  1:534408‑M1/43‑1 (MQ=255)
                   gcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATaattaatt                                                                                                                        <  1:165382‑M1/43‑1 (MQ=255)
                   gcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATaattaatt                                                                                                                        <  1:493059‑M1/43‑1 (MQ=255)
                   gcgtagatgagagtcgaagttacgttccgctagagtagaacgttgccaggcagttgtttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATaattaatt                                                                                                                        >  2:710688‑M1/98‑140 (MQ=255)
                    cgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTc                                                                                                                       <  2:630748‑M1/44‑1 (MQ=255)
                    cgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcaggtttttaatcaacttttggttaaaaactaaaatggccaagataaAAAAAGTTATTGACTTAAATGTTGATAAAAGGTAtcattaattt                                                                                                                       <  1:271655‑M1/44‑2 (MQ=255)
                     aaGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCt                                                                                                                      <  2:687682/140‑1 (MQ=255)
                         aGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATGAATTTTAAAAAGTTATTGCCTTAAATGTTAATAAAATGTATAATTAATTCTTGTc                                                                                                                  >  2:92692/1‑140 (MQ=255)
                           ggtagtcgaacttacgttccgctagagtagaacgttgccagacagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTTTTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCgg                                                                                                                <  1:547527‑M1/51‑1 (MQ=255)
                            gtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGt                                                                                                               <  1:104246‑M1/52‑1 (MQ=255)
                            gtGAAAGCCTGGATAGTATTTATTTGTTCAGTCAAACAAAACAGCTCAACGTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGt                                                                                                               <  1:23259/139‑1 (MQ=255)
                              gAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGCaa                                                                                                             >  2:498110/1‑140 (MQ=255)
                                tcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGCCTTAATTGTTGATAAAATGTCTAATTAATTcttctcgccaaga                                                                                                           <  1:66511‑M1/56‑7 (MQ=255)
                                     cttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAAt                                                                                                      >  1:665343‑M1/80‑140 (MQ=255)
                                       tatcGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATCGACTTAACTGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGa                                                                                                    >  2:62780/1‑140 (MQ=255)
                                       taagttccgctagaagagaacgttgacaggcagttttttaatcaaatattggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGa                                                                                                    >  2:550286‑M1/78‑140 (MQ=255)
                                        acgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTAGGTAAAAAAAATGaa                                                                                                   <  2:20945‑M1/64‑1 (MQ=255)
                                         taGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAAc                                                                                                  >  2:597737/1‑140 (MQ=255)
                                            tccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACAtt                                                                                               <  1:674302‑M1/68‑1 (MQ=255)
                                            tatttattTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAAGTTGAGTTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACAtt                                                                                               <  1:218856/140‑1 (MQ=255)
                                            tatttattTGTTGAAGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACAtt                                                                                               <  2:358499/140‑1 (MQ=255)
                                             atttatttGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTg                                                                                              >  1:606577/1‑140 (MQ=255)
                                              tttatttGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACAttaa                                                                                             >  2:631370/1‑138 (MQ=255)
                                               gctagaggagaacgttgccaggcaggttattaatcaaattctggttaaaaaataaaatggacaagataaAAAACGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTAGTCGGTAAGAAAAATCAACATTGaa                                                                                            <  2:73134‑M1/71‑1 (MQ=255)
                                                cgagaatacaacgttgcaaagcagatctttaatcaaattttgattaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATCCTTATCGGAGAGAAAAATGAACATTGaaa                                                                                           >  2:660023‑M1/69‑140 (MQ=255)
                                                 tagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGaaaa                                                                                          <  2:298982‑M1/73‑1 (MQ=255)
                                                   ttgttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACt                                                                                        <  2:78121/140‑1 (MQ=255)
                                                    tgttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTGGTCGGTAAGAAAAATGAACATTGAAAAGTg                                                                                       >  1:27463/1‑140 (MQ=255)
                                                            gttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACaa                                                                               <  1:179761‑M1/84‑1 (MQ=255)
                                                            gttgccaggcagttttttaataaaattatgataaaagaataaaatggacaagataaAAAAAGTTATTGACTTAAGTGTTGATAAAATGTATAATTAATTCTTGTCGGTAGGAAAAATGAACATTGAAAACTGAATGACaa                                                                               <  2:205767‑M1/84‑1 (MQ=255)
                                                            ggtgccaggcagttttttaacactatatgggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACaa                                                                               >  2:126867‑M1/57‑140 (MQ=255)
                                                            caaacaaaACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACaa                                                                               >  2:105849/1‑140 (MQ=255)
                                                              aacaaaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAata                                                                             >  2:326718/1‑140 (MQ=255)
                                                               gccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGGCTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAatat                                                                            <  1:376424‑M1/87‑1 (MQ=255)
                                                               gccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTACGAAAAATGAACATTGAAAACTGAATGACAatat                                                                            <  2:35696‑M1/87‑1 (MQ=255)
                                                               acaaaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAatat                                                                            <  1:419246/140‑1 (MQ=255)
                                                                 caggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGt                                                                          >  1:512308‑M1/52‑140 (MQ=255)
                                                                 caggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGt                                                                          <  1:356932‑M1/89‑1 (MQ=255)
                                                                 aaaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAATATTTATTTTAAATAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGt                                                                          <  1:42523/140‑1 (MQ=255)
                                                                 aaaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACCTTGAAAACTGAATGACAATATGt                                                                          >  2:442816/1‑140 (MQ=255)
                                                                   ggcagttttttcatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCa                                                                        >  2:192083‑M1/50‑140 (MQ=255)
                                                                   aacaacTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCa                                                                        >  1:480267/1‑140 (MQ=255)
                                                                    gcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGCCTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAAGGAACATTGCAAACTGAATGACAATATgtcca                                                                       <  1:152965‑M1/92‑3 (MQ=255)
                                                                    gcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTACTTGAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGACCATTGAAAACTGCATCACAATATGTCaa                                                                       <  2:96411‑M1/92‑1 (MQ=255)
                                                                     caactcaactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGCAGGTAAGAAAAATGAAGATTGAAGACTGAATGACAATATGTCaaa                                                                      >  2:516368/1‑139 (MQ=255)
                                                                      aactcaactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACg                                                                     >  2:393812/1‑140 (MQ=255)
                                                                            actTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAtt                                                               >  2:768504/1‑140 (MQ=255)
                                                                             taatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCACCGTTAATTc                                                              <  1:264300‑M1/101‑1 (MQ=255)
                                                                              ttAAGAAATGACTTGATTTACTAACGAAAAGTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTcc                                                             <  2:471639/140‑1 (MQ=255)
                                                                               atcaaattttagttaaaaaataaaatggacaagataaAAAAAGTTATCGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCa                                                            >  1:494944‑M1/38‑140 (MQ=255)
                                                                                aaGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAATTGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaa                                                           <  2:400749/140‑1 (MQ=255)
                                                                                aaGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaa                                                           <  1:121482/140‑1 (MQ=255)
                                                                                 aGAAATAATTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaa                                                          <  1:323976/140‑1 (MQ=255)
                                                                                 aGAAATAACTTGAATTACTAACGAAAATTAATTATAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaa                                                          <  1:323748/140‑1 (MQ=255)
                                                                                    aaTAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaaaa                                                       <  1:102789/140‑1 (MQ=255)
                                                                                       aaCTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGt                                                    <  1:630515/140‑1 (MQ=255)
                                                                                       aaCTTGAATTACTAACGAAAATTAATTTGAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGt                                                    <  2:179773/140‑1 (MQ=255)
                                                                                         ggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTaa                                                  <  1:475466‑M1/113‑1 (MQ=255)
                                                                                         ggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTaa                                                  <  2:709992‑M1/113‑1 (MQ=255)
                                                                                         ggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTaa                                                  <  1:281411‑M1/113‑1 (MQ=255)
                                                                                         cTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTaa                                                  <  1:514816/140‑1 (MQ=255)
                                                                                          ttGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAAc                                                 <  2:374026/140‑1 (MQ=255)
                                                                                            gAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGCAAAAAAAACTCTGATAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACta                                               >  1:628387/1‑140 (MQ=38)
                                                                                            gAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGATAATATGTCAATGTTAATTCCAAAAAACGTAACta                                               >  1:88289/1‑140 (MQ=255)
                                                                                              aTTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACtata                                             >  2:380596/1‑140 (MQ=255)
                                                                                                 aCTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATaaat                                          >  2:262139/1‑138 (MQ=255)
                                                                                                  cTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGtt                                         >  1:692472/1‑140 (MQ=255)
                                                                                                    aaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATaagtaac                                       <  1:205641‑M1/124‑4 (MQ=255)
                                                                                                      caaaaaTTAATTGTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaa                                     <  1:145726/138‑1 (MQ=255)
                                                                                                         aaaTTAATTATAAAAAGTTATTGAGTGAAATGTTAATTAAATGTATAATTAATTCTGGGCGGTAAGAAAAATGAACATTGAAAAATGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACa                                  <  2:313365/140‑1 (MQ=21)
                                                                                                          aattaattTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACag                                 >  2:269416/1‑139 (MQ=255)
                                                                                                            caagtaaaAAAAAGTTATTGACGTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACatta                               >  2:70092‑M1/9‑140 (MQ=255)
                                                                                                             taactaaAAAAAGTTATTGCCTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAATATTAACATTGAAATCTGAATCACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACattat                              >  2:12947‑M1/8‑140 (MQ=255)
                                                                                                               tcgggaaaaaTTGATTGACTAAGAAGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATGCCAAAAAACGTAACTATAAGTTACAAACATTAttt                            <  2:773897/135‑1 (MQ=21)
                                                                                                               attTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTAttt                            <  1:288935/140‑1 (MQ=255)
                                                                                                                ataaAAAAAGTTATTGATTGAAAAGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTa                           >  2:45078‑M1/5‑140 (MQ=255)
                                                                                                                 tttAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAg                          <  1:749091/140‑1 (MQ=255)
                                                                                                                  ttAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGt                         <  1:202401/140‑1 (MQ=255)
                                                                                                                  gaaaaaaGATATTGACTTAAATGGCGATAAAATGTATAAGTAATTCTTGTCGGTGAGAAAAATGAACATTGAAAACTGAATGACAATAGGGCAACGTTCCTCCCATAATACGTCACTGTGAGATAGACACAGTGTTTAGt                         >  2:199395/3‑140 (MQ=16)
                                                                                                                  aaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGt                         >  1:709992‑M1/3‑140 (MQ=255)
                                                                                                                   aAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTa                        >  1:460499‑M1/2‑140 (MQ=255)
                                                                                                                        atttATTGACTTAAATGTTAATAAAATGTATTATTAATACAACTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTAt                   <  2:7238/138‑1 (MQ=18)
                                                                                                                        aGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGAAACTATAAGTTACAAACACTATTTAGTATTTAt                   >  2:130232/1‑140 (MQ=255)
                                                                                                                          ttATTGACGTAAATGTTGATAAAATGTATAATTAATTCTTGTCTGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGa                 <  1:312062/140‑1 (MQ=255)
                                                                                                                           tATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAg                >  2:685906/1‑140 (MQ=255)
                                                                                                                            aTTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTCTCTAGTATTTATGAGc               >  1:234223/1‑140 (MQ=255)
                                                                                                                              tGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTa             <  2:297054/140‑1 (MQ=255)
                                                                                                                               gACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTaa            <  1:89440/140‑1 (MQ=255)
                                                                                                                                aCTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAAt           <  1:266715/140‑1 (MQ=255)
                                                                                                                                  ttAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCa         <  2:169458/140‑1 (MQ=255)
                                                                                                                                     aaaTGTTGATATAATGTATAACTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCaaa       <  2:162968/140‑1 (MQ=255)
                                                                                                                                          ttAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcatca  >  2:557500/1‑140 (MQ=255)
                                                                                                                                            |                                                                                                                                         
GGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCA  >  USA300TCH1516_ALE/512504‑512780

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: