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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE512,6430AG43.0% 80.0 / 213.2 184intergenic (+711/‑6052)recR/speA_2Recombination protein RecR/Arginine decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (52/52);  new base G (43/36);  total (96/88)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

Only 100 of 188 total aligned reads displayed.

GTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATAA  >  USA300TCH1516_ALE/512505‑512783
                                                                                                                                           |                                                                                                                                            
gtgtAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTtaat                                                                                                                                            >  1:889986/1‑141 (MQ=255)
gtgtAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTtaat                                                                                                                                            >  1:575340/1‑141 (MQ=255)
    aagttaagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATaaa                                                                                                                                         <  1:692841‑M1/29‑1 (MQ=255)
     aTTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAAt                                                                                                                                       >  1:840347/1‑141 (MQ=255)
      tttttACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATg                                                                                                                                      >  1:135179/1‑141 (MQ=255)
       ttttACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGt                                                                                                                                     >  1:670282/1‑141 (MQ=255)
           gctccttagcgtcgatggtagtcgaacttacgttgcgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGtata                                                                                                                                  <  1:157165‑M1/36‑1 (MQ=255)
           gctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGtata                                                                                                                                  >  1:885633‑M1/106‑141 (MQ=255)
           aCTGTTGTTAAGCCGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaa                                                                                                                                 <  2:408041/141‑1 (MQ=255)
            ctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcaggtttttaaccaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGGTGATAAAATGGAtaa                                                                                                                                 <  2:418677‑M1/37‑1 (MQ=255)
            cTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAACACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaat                                                                                                                                >  1:995522/1‑141 (MQ=255)
            cTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaat                                                                                                                                >  1:752810/1‑141 (MQ=255)
              gttgttAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGTGGTAAACAAAACAACTTAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAATATTTTTTGACTTAAATGTTAATAAAATGTAtaatta                                                                                                                              <  2:331853/141‑1 (MQ=255)
              gttgttAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaatta                                                                                                                              >  2:781025/1‑141 (MQ=255)
               cttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTAtaatta                                                                                                                              <  2:194429‑M1/40‑1 (MQ=255)
                ttagcgtggatggtagtcgaacttacgttccgctagagtatatcgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTAtaattaa                                                                                                                             >  2:578953‑M1/101‑141 (MQ=255)
                tgagcatccaaggaagtcgaacttacgttccgctagagtagaacgatgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTAtaattaa                                                                                                                             >  2:522226‑M1/101‑141 (MQ=255)
                 tagcgtccatagtagtcgaactaacattccgctagagtagaacgatgccaagcagttttttaatcaaattttggttaaaaaataaaaaaaacaaaaaaaAAAAAGTTATTGACTTAAATGTTGATAAAATGTAtaattaat                                                                                                                            >  2:355950‑M1/100‑141 (MQ=255)
                 gttAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATaattaatt                                                                                                                           <  1:184816/141‑1 (MQ=255)
                  ttAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTc                                                                                                                          >  2:171844/1‑141 (MQ=255)
                    tgtcgatggtagtcgaacttacgttccgctagattagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCt                                                                                                                         >  1:366003‑M1/97‑141 (MQ=255)
                          tttGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGt                                                                                                                  >  2:799132/1‑141 (MQ=255)
                          tttGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGt                                                                                                                  >  2:268494/1‑141 (MQ=255)
                           ggtagtcgaacttacgttccgctagagtaggacgttgccaggcagttttttaaacaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTATTGTcggg                                                                                                                  <  2:647066‑M1/52‑2 (MQ=255)
                                  gaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGaaaa                                                                                                           >  1:263823‑M1/83‑141 (MQ=255)
                                    acttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAAt                                                                                                         >  1:618783‑M1/81‑141 (MQ=255)
                                      tataGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGaa                                                                                                      >  1:85946/1‑141 (MQ=255)
                                        acgtgccgccagagtagaaagatgccaggcagttttttaatcaacttttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAAc                                                                                                     >  2:543900‑M1/77‑141 (MQ=255)
                                         cgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACa                                                                                                    <  1:419001‑M1/66‑1 (MQ=255)
                                          gtcccgctagagtagaacgttgccaggcagttttttaataaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACAt                                                                                                   >  2:210589‑M1/75‑141 (MQ=255)
                                          gTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATCATTAATTCTTGTCGGTAAGAAAAATGAACAtt                                                                                                  >  1:123940/1‑141 (MQ=255)
                                           tatttattTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTTTCGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCTGTAAGAAAAATGAACATTg                                                                                                 <  1:637475/141‑1 (MQ=255)
                                            tccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTg                                                                                                 >  1:985107‑M1/73‑141 (MQ=255)
                                               tatttGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGaaaa                                                                                             <  2:725871/141‑1 (MQ=255)
                                                 tttGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACt                                                                                           <  2:808487/141‑1 (MQ=255)
                                                  agagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAGGAAAAATGAACTttgaaagct                                                                                           <  2:285991‑M1/75‑4 (MQ=255)
                                                   gagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTg                                                                                          <  2:516455‑M1/76‑1 (MQ=255)
                                                   aagtagaacgttgctaggcagttttttaataaaaaaaaaaaaaaaaaataaaatggacaagatagAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTg                                                                                          >  2:790820‑M1/66‑141 (MQ=255)
                                                      tagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAACTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAAt                                                                                       >  1:303784‑M1/63‑141 (MQ=255)
                                                      tagaaagttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAAt                                                                                       >  1:1063505‑M1/63‑141 (MQ=255)
                                                       gAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATCTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTCAGAAAACTGAACATCGAAAACTGAATGa                                                                                     >  2:977074/1‑141 (MQ=255)
                                                          gCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACaa                                                                                  >  1:234192/1‑141 (MQ=255)
                                                            aaacaaaacaACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAata                                                                                <  1:522628/141‑1 (MQ=255)
                                                               gccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATg                                                                              <  1:1063202‑M1/88‑1 (MQ=255)
                                                               aaaaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGt                                                                             <  1:501733/140‑1 (MQ=255)
                                                               aaaaacaaATCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGt                                                                             <  1:484692/140‑1 (MQ=255)
                                                                aaaacaaACTCAACCTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGTTAAGAAAAATGAACATTGAAAACTGAATGACAATATGt                                                                             <  2:833327/141‑1 (MQ=255)
                                                                 caggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTc                                                                            <  1:906619‑M1/90‑1 (MQ=255)
                                                                  aggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCa                                                                           >  1:413362‑M1/51‑141 (MQ=255)
                                                                   ggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAACTAATTCTTGTCGGTAAGAAAAATGACCATTGAACACTGAATGACAATATGTCaa                                                                          >  1:686641‑M1/50‑141 (MQ=255)
                                                                           tttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAtt                                                                  <  2:310007‑M1/100‑1 (MQ=255)
                                                                           aactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTAGAAAGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAtt                                                                  <  2:844616/141‑1 (MQ=255)
                                                                                aaGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaa                                                             <  2:1056437/141‑1 (MQ=255)
                                                                                 aGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAAGTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaa                                                            <  1:796731/141‑1 (MQ=255)
                                                                                   aattttggttaaaaaataaactgggccctataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaaaa                                                          <  2:1023472‑M1/108‑1 (MQ=255)
                                                                                   aaaTAACTTGAATTACAAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaaaa                                                          >  1:528286/1‑141 (MQ=255)
                                                                                    aaTAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAAc                                                         <  2:889493/141‑1 (MQ=255)
                                                                                    aaTAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAAc                                                         >  1:940384/1‑141 (MQ=255)
                                                                                        tggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTaa                                                     <  1:446415‑M1/113‑1 (MQ=255)
                                                                                         ggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACgtacc                                                    <  2:668611‑M1/114‑3 (MQ=255)
                                                                                         ggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGtaaa                                                    <  2:613936‑M1/114‑2 (MQ=255)
                                                                                         cTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAAc                                                    >  1:656808/1‑141 (MQ=255)
                                                                                          ttGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCAAAAAAACGTaaat                                                   >  2:367375/1‑139 (MQ=255)
                                                                                          gtaaaaaaataaaatggacaagaaaaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACt                                                   >  2:841629‑M1/27‑141 (MQ=255)
                                                                                           tGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACta                                                  >  1:804847/1‑141 (MQ=255)
                                                                                             aaTTACTAACGAAAAATAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAAAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACaata                                                >  1:564778/1‑141 (MQ=255)
                                                                                              aTTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATaa                                               <  1:73558/141‑1 (MQ=255)
                                                                                              aTTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATaa                                               <  1:696190/141‑1 (MQ=255)
                                                                                              aTTACTAACAAAAAAAAATATAAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACATTAAAAAAAAAAAACGTAACTATaa                                               <  2:452814/141‑1 (MQ=37)
                                                                                                tACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGt                                             <  2:257327/141‑1 (MQ=255)
                                                                                                tACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATAACAATATGTCAACGTTAATTACAAAAAACGTAACTATaatt                                             >  2:506518/1‑139 (MQ=255)
                                                                                                  cTAACGAAAATTAATTCTAAATAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGGACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACATAACTATAAGTTa                                           >  1:546474/1‑141 (MQ=255)
                                                                                                   taaaattaacaagataaAAAAAGTTATTAACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCTGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTAc                                          >  1:530894‑M1/18‑141 (MQ=255)
                                                                                                   aaaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTAc                                          >  1:292412‑M1/18‑141 (MQ=255)
                                                                                                    aatgAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAGAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACa                                         <  2:722495/138‑1 (MQ=255)
                                                                                                    aaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACa                                         >  1:518425‑M1/17‑141 (MQ=255)
                                                                                                    aaCGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACa                                         <  2:213663/141‑1 (MQ=255)
                                                                                                     aCGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaa                                        >  1:761024/1‑141 (MQ=255)
                                                                                                      cGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaaa                                       >  2:681313/1‑141 (MQ=255)
                                                                                                      aacggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaaa                                       >  1:516455‑M1/15‑141 (MQ=255)
                                                                                                       gAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAGAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAAc                                      <  1:31809/141‑1 (MQ=255)
                                                                                                             taattTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACattatt                                <  1:171844/141‑1 (MQ=255)
                                                                                                              aattTTAAACAGGTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAAGATTGAAAACTGAATGACAATATGTCAACGTTAATCCCAAAAACCGTAACTTTAAGTTACACAAATTAttt                               >  2:410660/1‑141 (MQ=25)
                                                                                                                ttttAAAAAGTTGTTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAg                             <  2:503994/141‑1 (MQ=255)
                                                                                                                ttttAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAg                             <  2:940384/141‑1 (MQ=255)
                                                                                                                  aaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTa                           <  1:352464‑M1/139‑1 (MQ=255)
                                                                                                                  aaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTa                           >  2:88133‑M1/3‑141 (MQ=255)
                                                                                                                  aaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTa                           <  1:603003‑M1/139‑1 (MQ=255)
                                                                                                                   tAAAAAGGTATTGGCTTAAATGTTAATAAAATGTATAATTAATTCGTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTAt                          <  2:103798/141‑1 (MQ=255)
                                                                                                                   aAAAAAGTTATTGAATTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATACGTCAAAGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTCGTAt                          >  2:645901‑M1/2‑141 (MQ=255)
                                                                                                                     aaaaGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTAttt                        >  2:493703/1‑141 (MQ=255)
                                                                                                                         gTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGa                    <  1:198880/141‑1 (MQ=255)
                                                                                                                         gTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGa                    <  1:969104/141‑1 (MQ=255)
                                                                                                                         gTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAGTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGa                    <  1:470812/141‑1 (MQ=255)
                                                                                                                            aTTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCt                 >  2:466791/1‑141 (MQ=255)
                                                                                                                               gACTGAAATGTTGATAAAAGGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTAAGAACTGAATGACAATATGTCAACGTTAATTACAAAAAACGTAACTATACGTTACAAACATTATTTATTACTTATGACCTAAt              >  2:594629/1‑141 (MQ=18)
                                                                                                                                 cTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTTAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCa            <  2:146181/141‑1 (MQ=255)
                                                                                                                                   aaaaTGTAAATAAAATGTATAATCAATTCTAGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCaaa          <  2:974343/140‑1 (MQ=255)
                                                                                                                                      aTGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcat       <  2:234192/141‑1 (MQ=255)
                                                                                                                                         ttAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcatcat    >  1:397209/1‑141 (MQ=255)
                                                                                                                                           |                                                                                                                                            
GTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATAA  >  USA300TCH1516_ALE/512505‑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: