Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 512,643 A→G 54.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,6430AG54.5% 6.1 / 142.9 108intergenic (+711/‑6052)recR/speA_2Recombination protein RecR/Arginine decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (22/27);  new base G (29/30);  total (51/57)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.65e-02

Only 100 of 110 total aligned reads displayed.

AATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATA  >  USA300TCH1516_ALE/512509‑512782
                                                                                                                                      |                                                                                                                                           
aaTTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATaaa                                                                                                                                        <  1:113707/140‑1 (MQ=255)
 gttaagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATaaaa                                                                                                                                       >  1:235998‑M1/111‑140 (MQ=255)
    aagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGt                                                                                                                                    <  2:450685‑M1/33‑1 (MQ=255)
     agctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGta                                                                                                                                   <  1:606874‑M1/34‑1 (MQ=255)
      tACTGTTGTTAAACAAATTGAAAACCTGTAAAGTAATTATATGTCGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACTAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGtat                                                                                                                                  <  2:516578/140‑1 (MQ=255)
      gctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGtat                                                                                                                                  <  2:135017‑M1/35‑1 (MQ=255)
      gctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGtat                                                                                                                                  <  1:476084‑M1/35‑1 (MQ=255)
      gctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGtat                                                                                                                                  <  2:513726‑M1/35‑1 (MQ=255)
      gctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttgtttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTGAAATGTTGATAGaatgaat                                                                                                                                  <  2:563080‑M1/35‑4 (MQ=255)
      gctccttagcgtcgatagtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaatttgggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATGAaatggat                                                                                                                                  <  1:102416‑M1/35‑4 (MQ=255)
        cTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAGAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaa                                                                                                                                >  1:161138/1‑140 (MQ=255)
          gttgttAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaatt                                                                                                                              <  2:625662/140‑1 (MQ=255)
          cttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTAtaatt                                                                                                                              <  1:406918‑M1/39‑1 (MQ=255)
           ttgttAAGCAGTTTGGAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaatta                                                                                                                             <  2:308503/140‑1 (MQ=255)
           ttgttAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaatta                                                                                                                             >  1:401481/1‑140 (MQ=255)
           ttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTTGATGTTGATAAAATGTAtaatta                                                                                                                             <  1:186553‑M1/40‑1 (MQ=255)
            tgttAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaattaa                                                                                                                            <  2:588710/140‑1 (MQ=255)
            tagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTAtaattaa                                                                                                                            >  2:525190‑M1/100‑140 (MQ=255)
             agcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTAtaattaat                                                                                                                           <  2:182203‑M1/42‑1 (MQ=255)
                gtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCt                                                                                                                        <  2:703338‑M1/45‑1 (MQ=255)
                  gCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTg                                                                                                                      >  1:402222/1‑140 (MQ=255)
                     gTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCg                                                                                                                   >  1:537853/1‑140 (MQ=255)
                      tttGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCgg                                                                                                                  >  2:429464/1‑140 (MQ=255)
                      ggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTACATGTTGATAAAATGTATAATTAATTCTTGTCgg                                                                                                                  <  2:450706‑M1/51‑1 (MQ=255)
                       gtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGTTTAATTTTATAATTAATTATTGTCGGt                                                                                                                 <  2:278282‑M1/52‑1 (MQ=255)
                       gtagtcgaacttacgttacgctagagtagaacggtgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGt                                                                                                                 >  1:374830‑M1/89‑140 (MQ=255)
                         gAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTaa                                                                                                               >  2:406165/1‑140 (MQ=255)
                         agtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTaa                                                                                                               <  1:53704‑M1/54‑1 (MQ=255)
                          gtcgaactgacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAg                                                                                                              >  1:17692‑M1/86‑140 (MQ=255)
                          aaaGCCTGTATAGTATTTATTTGTCGAGGCAACCAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAg                                                                                                              >  2:613767/1‑140 (MQ=255)
                          aaaGCCGGTAGACTATGTATGTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAAGTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAAGTCTTGTCGGTAAg                                                                                                              <  1:526729/140‑1 (MQ=255)
                            aGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGaa                                                                                                            >  1:382142/1‑140 (MQ=255)
                                tGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAACAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAAt                                                                                                        <  2:11140/140‑1 (MQ=255)
                                 gTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATg                                                                                                       >  1:681869/1‑140 (MQ=255)
                                      gTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAATAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACAt                                                                                                  <  2:574964/140‑1 (MQ=255)
                                            tagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGaaaa                                                                                            >  2:240777‑M1/68‑140 (MQ=255)
                                            atttGTTGAGGAAAACAAAACAACGCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGaaaa                                                                                            <  2:238275/140‑1 (MQ=255)
                                            aagcgaagaaagttgccaggcagttgattaatcaaatttgggttaaaaaacaaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGaaaa                                                                                            >  2:322801‑M1/68‑140 (MQ=255)
                                              ttgttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTCATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGACGGTAAGAAAAATGAACATTGAAAACt                                                                                          >  1:329989/1‑140 (MQ=255)
                                              ttgttgAGGAAAACAAAACCGCGCGGCTCAAAAAATAACGTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACt                                                                                          <  1:126236/140‑1 (MQ=255)
                                               agtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTg                                                                                         <  1:668690‑M1/76‑1 (MQ=255)
                                                gttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGa                                                                                        >  2:490672/1‑140 (MQ=255)
                                                 ttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTTTCGGTAAGAAAAATGACCATTGAAATCTgga                                                                                       >  1:466015/1‑138 (MQ=255)
                                                       caaacaaaACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACaa                                                                                 <  1:498996/140‑1 (MQ=255)
                                                       caaacaaaACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACaa                                                                                 >  1:68835/1‑140 (MQ=255)
                                                          acaaaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAatat                                                                              >  1:247589/1‑140 (MQ=255)
                                                             aggcagttttgtaaacaaagtttagttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTc                                                                           >  2:468567‑M1/51‑140 (MQ=255)
                                                               gcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCaa                                                                         <  1:150661‑M1/92‑1 (MQ=255)
                                                                 agttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACg                                                                       >  1:68730‑M1/47‑140 (MQ=255)
                                                                 agtcttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACg                                                                       >  1:200755‑M1/47‑140 (MQ=255)
                                                                   ctcaactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGACAACGtt                                                                     >  2:434886/1‑140 (MQ=255)
                                                                     ttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTaa                                                                   <  2:194478‑M1/98‑1 (MQ=255)
                                                                        taatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTc                                                                <  2:477185‑M1/101‑1 (MQ=255)
                                                                        ctTAATAACTAACTACAATTACTAACTAAAATTAATTTTAAAAATTTATTGTCTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTTAAAACTGAATGACAATATGTCAACGTTAATTc                                                                <  2:77449/140‑1 (MQ=38)
                                                                           tcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAAGGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaa                                                             >  1:13119‑M1/37‑140 (MQ=255)
                                                                           acgAAATAACTTGAATTACTAACGAAAATTGATTTTAAAAAGTTATTGACTCAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaa                                                             <  1:263477/138‑1 (MQ=255)
                                                                           aaGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaa                                                             <  1:563672/140‑1 (MQ=255)
                                                                             gAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaa                                                           >  1:589052/1‑140 (MQ=255)
                                                                             acatattagttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTGAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaa                                                           >  2:100637‑M1/35‑140 (MQ=255)
                                                                              aatattagttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCCTGTAGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaaa                                                          >  2:499779‑M1/34‑140 (MQ=255)
                                                                               attttggtgaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaaaa                                                         >  2:366029‑M1/33‑140 (MQ=255)
                                                                               atcaaCTTTAATGACTAACGAATATTACTTTTAAAAAGTTATTGACTTAAATGTTAATAAAAAGTATAATTAATTATTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaaaa                                                         <  2:35896/137‑1 (MQ=37)
                                                                               aaTAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCTaaaaaa                                                         >  2:479617/1‑140 (MQ=255)
                                                                                ttttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAAc                                                        <  2:75046‑M1/109‑1 (MQ=255)
                                                                                 tttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACg                                                       >  1:477185‑M1/31‑140 (MQ=255)
                                                                                 tagcTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACg                                                       <  2:41992/137‑1 (MQ=255)
                                                                                        aaTTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACtat                                                <  1:210411/140‑1 (MQ=255)
                                                                                         aaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACtata                                               >  1:135017‑M1/23‑140 (MQ=255)
                                                                                          aaaataaaatggacaagataaAAAAAGTTATTGACGTAAAGGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAAGGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATaa                                              >  1:172762‑M1/22‑140 (MQ=255)
                                                                                          aaaataaaacggacaagagaaAAAAAGGTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATaa                                              >  1:383177‑M1/22‑140 (MQ=255)
                                                                                             ataaaatggacaagatcaAAAAAGTTATTGACTTAAATGTTGATAAAATGTCTAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAacgtaacttatagtt                                           >  2:186553‑M1/19‑133 (MQ=255)
                                                                                              taaaatggacaagaaaaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTa                                          <  2:696000‑M1/123‑1 (MQ=255)
                                                                                                 tgtggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAAAGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAATCGTAACTATAAGTTACaa                                       >  1:503178‑M1/15‑140 (MQ=255)
                                                                                                 aatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaa                                       >  1:317539‑M1/15‑140 (MQ=255)
                                                                                                 aagggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAAGTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaa                                       >  1:231146‑M1/15‑140 (MQ=255)
                                                                                                  gAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaaa                                      <  1:611093/140‑1 (MQ=255)
                                                                                                  gAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaaa                                      <  1:639175/140‑1 (MQ=255)
                                                                                                  gAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaaa                                      <  1:610743/140‑1 (MQ=255)
                                                                                                  gAAAATTAATTTTAAAAAGTCATTAACTTAAATGAAAAAAAAATAAATAATGAAGTCTTGTCGGAAAGAAAAATGAACATTGAAAAATGAATAAAAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaaa                                      <  2:350718/140‑1 (MQ=18)
                                                                                                  atggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaaa                                      >  1:703338‑M1/14‑140 (MQ=255)
                                                                                                      attaattTTAAAAAGTTATTGACTTCAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTAAAAGTAACAACCatt                                  >  2:534315/1‑140 (MQ=255)
                                                                                                      attaattTTAAAAAGTTATTGACTTAAATGTTAAGAAAATGCATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACatt                                  <  1:128432/140‑1 (MQ=255)
                                                                                                       ataattTTAAAAAGTTATTGACTAAAAAGTAAATAAAATGTATAAATAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACatta                                 <  2:402222/139‑1 (MQ=255)
                                                                                                        aagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTATTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACATTAATTCCAAAAAACGTAACTATAAGTTACAaacatgat                                <  2:162038‑M1/133‑4 (MQ=255)
                                                                                                         agataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAAAGTTAATTCCAAAAAACGTAACTATAAGTTACAAACattatt                               <  2:127039‑M1/134‑1 (MQ=255)
                                                                                                          attTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTAAAAAAAACGTAACTATAAGTTACAAACATTAttt                              >  2:232325/1‑140 (MQ=255)
                                                                                                             aaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAGTTCAAAAAAACGTAACTATAAGTTACAAACAttatttcga                           <  2:414379‑M1/138‑4 (MQ=255)
                                                                                                              aAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACAGTATTTAGTa                          <  1:312064‑M1/139‑1 (MQ=255)
                                                                                                              aAAAAAGTTATTGACTTAAATGTTGATAAAATGTAGAATCAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTa                          <  1:135623‑M1/139‑1 (MQ=255)
                                                                                                                 aaaGTTATTGGCTTAAATGTTGATAAAATGTATAATTACTTCTTATCAGTAAGAAAAATGAATATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTAttt                       <  1:226442/140‑1 (MQ=18)
                                                                                                                  aaGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTa                      <  1:319532/140‑1 (MQ=255)
                                                                                                                   aGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTAt                     <  2:510865/140‑1 (MQ=255)
                                                                                                                    gTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATg                    >  1:61702/1‑140 (MQ=255)
                                                                                                                       aTTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTTCACACATTATTTAGTATTAATGAGc                 >  2:457845/1‑140 (MQ=255)
                                                                                                                       aTTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATCTAGTATTTATGAGc                 >  2:139230/1‑140 (MQ=255)
                                                                                                                        ttGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGACAGCTGAATGATAATATGTCAACGTTAATTCGAACAACCGTAAGTATGAGGTACAAACGACATCTAGTATTTATgcgcc                >  2:614619/1‑139 (MQ=255)
                                                                                                                                aaTGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAc        <  2:61702/140‑1 (MQ=255)
                                                                                                                                   gTTGATAAAATGTATAATTAATTCTTGTCGTTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcatc     <  2:102416/140‑1 (MQ=255)
                                                                                                                                     taataaAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcatcat   <  1:376656/140‑1 (MQ=255)
                                                                                                                                      aataaAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACAACATa  >  2:330807/1‑140 (MQ=255)
                                                                                                                                      |                                                                                                                                           
AATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATA  >  USA300TCH1516_ALE/512509‑512782

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: