Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 512,643 A→G 45.8% 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,6430AG45.8% 48.0 / 252.3 193intergenic (+711/‑6052)recR/speA_2Recombination protein RecR/Arginine decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (51/52);  new base G (41/47);  total (94/99)
Fisher's exact test for biased strand distribution p-value = 7.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 204 total aligned reads displayed.

GGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATAA  >  USA300TCH1516_ALE/512504‑512783
                                                                                                                                           |                                                                                                                                            
ggggtAATTTTTACTGTTGTTAGGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTtaa                                                                                                                                             <  2:739575/138‑1 (MQ=255)
ggTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTtaa                                                                                                                                             <  1:637906/141‑1 (MQ=255)
 gtgtAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTtaat                                                                                                                                            <  1:445243/141‑1 (MQ=255)
 gtgtAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTtaat                                                                                                                                            <  1:278999/141‑1 (MQ=255)
   gtAATTTTTACTGTTGTTAAGTAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTtaataa                                                                                                                                          <  1:236603/141‑1 (MQ=255)
   gtAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGAAATTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAGATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTtaataa                                                                                                                                          >  1:711502/1‑141 (MQ=255)
   gtAATTTTTACTGTTGTTAAGCAGGTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAGGTTATTGACTTCAATGTtaataa                                                                                                                                          >  1:401959/1‑141 (MQ=255)
      aTTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAAt                                                                                                                                       >  1:564662/1‑141 (MQ=255)
      aTTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAAt                                                                                                                                       >  1:442137/1‑141 (MQ=255)
       ttaagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggaagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATg                                                                                                                                      >  1:788469‑M1/110‑141 (MQ=255)
          agctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGtat                                                                                                                                   <  1:195091‑M1/35‑1 (MQ=255)
              tgttgtTAATCAGTTTGAAATCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaatt                                                                                                                               <  2:442816/141‑1 (MQ=255)
              tgttgtTAAGCAGTTTTAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACATGAACGTTAATAAAATGTAtaatt                                                                                                                               <  2:883125/141‑1 (MQ=255)
                  agcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATaattaatt                                                                                                                           <  2:245671‑M1/43‑1 (MQ=255)
                    ggtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCt                                                                                                                         >  1:731933‑M1/97‑141 (MQ=255)
                       gCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGt                                                                                                                      >  1:649327/1‑141 (MQ=255)
                        gatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATCAAATGTATAATTAATTCTTGTc                                                                                                                     <  2:325325‑M1/49‑1 (MQ=255)
                           ggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaATAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGt                                                                                                                  <  2:624082‑M1/52‑1 (MQ=255)
                           ggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGt                                                                                                                  <  1:911220‑M1/52‑1 (MQ=255)
                            gtagtcgaacttacgttccgctagagtagaacgttgcgaggcagttttttaatcaaattttggttaaaaaataaaatagacaagataaAGAAAGTTATTGACTTAAATGTTGATAAAATGTATAATAAATTCTTGTCGGTa                                                                                                                 <  2:209842‑M1/53‑1 (MQ=255)
                            gtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTa                                                                                                                 <  1:880029‑M1/53‑1 (MQ=255)
                              agtcgaacttacgttccgctagagtagaacgttgccaggcagttctttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAg                                                                                                               <  2:486338‑M1/55‑1 (MQ=255)
                               aaaGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGa                                                                                                              <  2:147082/141‑1 (MQ=255)
                                aaGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGaa                                                                                                             >  2:23190/1‑141 (MQ=255)
                                aaGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGaa                                                                                                             <  1:158027/141‑1 (MQ=255)
                                aaGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATATATAATTAATTCTTGTGGGTAAGaa                                                                                                             >  1:892502/1‑141 (MQ=255)
                                    acttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAAt                                                                                                         <  2:887429‑M1/61‑1 (MQ=255)
                                         taGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACa                                                                                                    <  1:335063/141‑1 (MQ=255)
                                         taGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACa                                                                                                    <  2:347968/141‑1 (MQ=255)
                                         cgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACa                                                                                                    <  2:899446‑M1/66‑1 (MQ=255)
                                           gTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTACCGAAAATTAATTTTAAAAAGTTATTGACTTCAATGTTAATAACATGTATAATTAATTCTGGTCGGTAAGAAAAATGAACAtt                                                                                                  >  1:344131/1‑141 (MQ=255)
                                           gTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGACGGTAAGCAAAATGAACAtt                                                                                                  >  1:630088/1‑141 (MQ=255)
                                            tccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACattt                                                                                                 <  2:179104‑M1/69‑2 (MQ=255)
                                            tccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTg                                                                                                 <  1:298685‑M1/69‑1 (MQ=255)
                                             atttatttGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGa                                                                                                >  2:464400/1‑141 (MQ=255)
                                                ctagagtaggacgttgccaggcaggtttttaatcaaatttaggttaaaaactaaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATtgaaca                                                                                             <  1:667419‑M1/73‑3 (MQ=255)
                                                 tagggtagaacgttgccaggcagttttttaattaaattttggttaaaaaataaaatagaaaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAAc                                                                                            >  2:338589‑M1/68‑141 (MQ=255)
                                                   ttgttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTg                                                                                          >  1:376614/1‑141 (MQ=255)
                                                     gttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGaa                                                                                        >  2:608864/1‑141 (MQ=255)
                                                       tgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATg                                                                                      >  2:768457/1‑141 (MQ=255)
                                                       tgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATg                                                                                      >  2:668064/1‑141 (MQ=255)
                                                       agaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATg                                                                                      >  1:602041‑M1/62‑141 (MQ=255)
                                                        gaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGa                                                                                     <  2:364695‑M1/81‑1 (MQ=255)
                                                           cgttgccaggcagttttttaatcaaattttggttaaaaaataaaatcgacaagataaAAAAAGTTATTGACTTAACTGTTGATAAAATGTATAATTAATTCTTGTCGGTCAGATAAACGAACATTGAAAACTGAATGACaa                                                                                  >  2:645395‑M1/58‑141 (MQ=255)
                                                             aaacaaaacaACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAata                                                                                >  2:334675/1‑141 (MQ=255)
                                                             aaacaaaacaACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAata                                                                                <  2:131028/141‑1 (MQ=255)
                                                              aaaaaaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAatat                                                                               <  1:166767/138‑1 (MQ=255)
                                                               acaaaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATg                                                                              >  1:724984/1‑141 (MQ=255)
                                                               acaaaacaaCTCAACTTAAGAAATAACTTGAATTACTAAAGAAAATTAATTTTTACTAGTTATTGACTTATATGTTCATATACTGAACAAGTAATACAATTCCGTATGAAAAATGAACATGAAAAACTGAATGACAATATg                                                                              >  2:78468/1‑141 (MQ=255)
                                                               acaaaacaaCTCAACCTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATg                                                                              >  1:675885/1‑141 (MQ=255)
                                                                   gacagttttttaataaaattttggttaaaaaataaaaagaaaaaaaaaaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCaa                                                                          >  2:427739‑M1/50‑141 (MQ=255)
                                                                       actcaactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGtt                                                                      >  1:282318/1‑141 (MQ=255)
                                                                       actcaactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGtt                                                                      >  2:109670/1‑141 (MQ=255)
                                                                             ctTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAATAGTTATTTACTTAAATGTTAATAAAATGTATAATTATTTTTTGTCGGTAAGAAAAATGAAGATTGAAAACTGAATGACATTATGTCAACGTTAATTcc                                                                <  2:39132/141‑1 (MQ=255)
                                                                             ctTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTcc                                                                <  2:515833/141‑1 (MQ=255)
                                                                                tcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaa                                                             >  1:485658‑M1/37‑141 (MQ=255)
                                                                                 aGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaa                                                            <  1:760510/141‑1 (MQ=255)
                                                                                 aGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaa                                                            <  1:23190/141‑1 (MQ=255)
                                                                                        tggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTaa                                                     <  1:822142‑M1/113‑1 (MQ=255)
                                                                                         ggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAAc                                                    <  2:274681‑M1/114‑1 (MQ=255)
                                                                                         ggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAAc                                                    <  1:308639‑M1/114‑1 (MQ=255)
                                                                                         ggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAAc                                                    <  2:847959‑M1/114‑1 (MQ=255)
                                                                                         ggttaaaaaataaaatgaacaagataaAAAAAGTTATTGACTTAAATGGTGATAAAATGTATAATTAATTCTTGTCGGTAAGACAAATCAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAAc                                                    >  2:643486‑M1/28‑141 (MQ=255)
                                                                                         cTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAAc                                                    <  1:320961/141‑1 (MQ=255)
                                                                                          tttaaaaaataaaatggacaagataaAAAAAGTTATCGACTTACCAGTTCAAAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACt                                                   >  1:518859‑M1/27‑141 (MQ=255)
                                                                                          ttGAATGACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACt                                                   <  2:290171/141‑1 (MQ=255)
                                                                                          gttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACgtaaat                                                   <  2:586393‑M1/115‑3 (MQ=255)
                                                                                          gttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACt                                                   >  1:585283‑M1/27‑141 (MQ=255)
                                                                                          gttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACCTTAATTCCAAAAAACGTAACt                                                   >  1:473711‑M1/27‑141 (MQ=255)
                                                                                             aaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACtata                                                >  1:354752‑M1/24‑141 (MQ=255)
                                                                                              aTTACTAACGAAAATTAATTTTAAAAAGTTATTGACTAAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATaa                                               <  2:381526/141‑1 (MQ=255)
                                                                                              aTTACTAAAGAAAATTTATTTTAAAAAGTTTGTGTCCAAAAGGTAAATACTATGTATAATTAACTCTTGTCAGTAAGAAAAATGAACATTGAAAACTGTATGGCAAGAGTTCAACGTTAATTCCAAAAAACGTAACTATaa                                               <  2:707154/141‑1 (MQ=16)
                                                                                               ttACTAATGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAg                                              <  2:72396/141‑1 (MQ=255)
                                                                                               ttACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAg                                              >  2:730383/1‑141 (MQ=255)
                                                                                                 aCTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGtt                                            <  2:59888/141‑1 (MQ=255)
                                                                                                  ataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTa                                           <  1:765785‑M1/123‑1 (MQ=255)
                                                                                                    aaaatgaacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACa                                         >  2:762916‑M1/17‑141 (MQ=255)
                                                                                                    aaaaaagaaaagaaaaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGAAAAAACCTAAAAACGTAACTATAAGTTACa                                         >  2:664843‑M1/17‑141 (MQ=255)
                                                                                                    aaCGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAACAACGTAACTATAAGTGACa                                         >  2:57774/1‑141 (MQ=255)
                                                                                                      aatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaaa                                       >  1:14284‑M1/15‑141 (MQ=255)
                                                                                                      aatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAACTTAaaaa                                       <  2:26466‑M1/127‑1 (MQ=255)
                                                                                                          gacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACatt                                   <  1:693478‑M1/131‑1 (MQ=255)
                                                                                                           attaattTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACatta                                  >  1:839098/1‑141 (MQ=255)
                                                                                                            ttaattTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAAAATTGAAATCTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACattat                                 >  2:399976/1‑141 (MQ=255)
                                                                                                              agataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTAAAAACAGAATGACAATATGTCAAAAAAAAAAACAAAAAACGTAACTATAAGTTACAAACATTAttt                               >  2:687820‑M1/7‑141 (MQ=255)
                                                                                                              aattTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAATTTACAAACATTAttt                               >  1:129389/1‑141 (MQ=255)
                                                                                                                ttttAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAg                             <  1:287418/141‑1 (MQ=255)
                                                                                                                 tttAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGt                            <  2:762132/141‑1 (MQ=255)
                                                                                                                  ttAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTa                           <  2:863234/141‑1 (MQ=255)
                                                                                                                  ttAAAAAGTTATTAACTTAAATGTTAATAAAATTTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTa                           <  1:415126/141‑1 (MQ=255)
                                                                                                                  aaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTa                           <  1:331058‑M1/139‑1 (MQ=255)
                                                                                                                   aAAAAAGTTATTGACTTAAAGGAGGAGAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACGTTGAAAACTGAATGACAATATGTCCACGTTAATTCCAAAAAACCTAACTATAAGTTACAAACAGTATTTAGTAt                          <  2:550473‑M1/140‑1 (MQ=255)
                                                                                                                    aaaaaGTAATTGACTTAAATGTTGATAACATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTAtt                         >  1:575662/1‑141 (MQ=255)
                                                                                                                     aaaaGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTAttt                        <  1:608864/141‑1 (MQ=255)
                                                                                                                        aGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATg                     >  1:718466/1‑141 (MQ=255)
                                                                                                                              tGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAATATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTaa               <  1:688609/141‑1 (MQ=255)
                                                                                                                              tGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGGGCTaa               >  1:116051/1‑141 (MQ=255)
                                                                                                                                   tAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTCCAAACATTATTTAGTATTTATGAGCTAATCaaa          >  2:756304/1‑141 (MQ=255)
                                                                                                                                       tGTTGATAAAATGTATAATTAATGCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcatc      <  2:448827/141‑1 (MQ=255)
                                                                                                                                         ttAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcatcat    >  1:529396/1‑141 (MQ=255)
                                                                                                                                           |                                                                                                                                            
GGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATAA  >  USA300TCH1516_ALE/512504‑512783

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: