Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 512,643 A→G 48.1% 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,6430AG48.1% 10.5 / 101.9 86intergenic (+711/‑6052)recR/speA_2Recombination protein RecR/Arginine decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (25/19);  new base G (22/19);  total (47/39)
Fisher's exact test for biased strand distribution p-value = 8.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCAT  >  USA300TCH1516_ALE/512508‑512781
                                                                                                                                       |                                                                                                                                          
tAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAAAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAATTATTGACATAAATGATAATaaa                                                                                                                                       >  2:473362/1‑141 (MQ=255)
aaaTTTTTAATGTTGTTAAGCAGTATGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTAATAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATaaa                                                                                                                                       <  2:478509/140‑1 (MQ=255)
       gctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGtata                                                                                                                                <  1:62310‑M1/36‑1 (MQ=255)
        aCTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTAATTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaa                                                                                                                               <  2:144053/141‑1 (MQ=255)
          tgtagtTAAGCAGGATGAAAGCGTGTAAAGTATTGATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaatt                                                                                                                             <  2:143795/141‑1 (MQ=255)
            ttgttAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaattaa                                                                                                                           <  1:367283/141‑1 (MQ=255)
              tttAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTGGAAGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATaattaatt                                                                                                                         <  1:425288/140‑1 (MQ=255)
               gcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttgattaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATgtataattaactt                                                                                                                        <  2:390574‑M1/44‑4 (MQ=255)
                       ggtagtcgaacttacgttccgctagagtagaacgttgcctggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAACAGTTATTGACTTAACTGTTGATAAAATGTAAAATTAATTCTtgtcagt                                                                                                                <  2:290614‑M1/52‑4 (MQ=255)
                       ggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGt                                                                                                                <  1:377259‑M1/52‑1 (MQ=255)
                       ggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGt                                                                                                                <  1:324211‑M1/52‑1 (MQ=255)
                       ggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGt                                                                                                                <  2:378747‑M1/52‑1 (MQ=255)
                        gtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTACTTCTTGTCGGTa                                                                                                               <  2:230309‑M1/53‑1 (MQ=255)
                          agtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAg                                                                                                             >  1:23815‑M1/87‑141 (MQ=255)
                           aaaGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATAAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGa                                                                                                            >  2:85843/1‑141 (MQ=255)
                            tcaaacgtacattccgctagagtagaacggtggcaggcagttttgtaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGaa                                                                                                           >  2:408326‑M1/85‑141 (MQ=255)
                            aaGCCTGTATAGTATTTATTTGTTGAGGCAAACACAACAACTCAACTTAAGAAATAACTTGAATTACTACCGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGaa                                                                                                           >  2:406681/1‑141 (MQ=255)
                              gaacataagttctgataaagtagaacgtagacaagcagttttttaaacaaattttagttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCATGTCTGTAAGaaaa                                                                                                         >  2:436844‑M1/83‑141 (MQ=255)
                                  gTATAGTATTTATTTGCTGAGGCAAACAAAACCACTCAACTTCAGAAATAACAAGAACTACTAAAGAAACTTAATCTTAAAAAGTTATTGACTTCAATGTTAATAAAATGTATAATTAATTCTTACAGGTAACAAAAATGa                                                                                                     >  2:84702/1‑141 (MQ=255)
                                    acgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATgaaa                                                                                                   <  1:241981‑M1/65‑2 (MQ=255)
                                     cgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAacatggagg                                                                                                  <  1:489718‑M1/66‑5 (MQ=255)
                                       gTATTTATTTGTTGAGGCAAACAAAACAACTCGACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTGAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTATGAAAAATGAACAtt                                                                                                >  2:439437/1‑141 (MQ=255)
                                        tccgctagagtagaacgttgacaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTg                                                                                               >  1:425586‑M1/73‑141 (MQ=255)
                                         atttatttGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGTTAAGACAAATGAACAtttg                                                                                              >  1:40551/1‑139 (MQ=255)
                                          cgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGaa                                                                                             <  2:186744‑M1/71‑1 (MQ=255)
                                            ctagagtagaacattgccaggcagtttattaatcaaatattggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGaaaa                                                                                           >  2:18335‑M1/69‑141 (MQ=255)
                                              agagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGACAACt                                                                                         <  1:24899‑M1/75‑1 (MQ=255)
                                              agagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACt                                                                                         >  1:170998‑M1/67‑141 (MQ=255)
                                               ttgttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAACAGTTATTGCCTTAAACGTAAATACAATGTATAATTAATTCTTGTCGCAAAGAAAAATGAACATGGAAAACtc                                                                                        >  2:128135/1‑140 (MQ=255)
                                                tgttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGa                                                                                       >  1:431577/1‑141 (MQ=255)
                                                tgttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGa                                                                                       >  1:495586/1‑141 (MQ=255)
                                                  tagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAAt                                                                                     <  1:198835‑M1/79‑1 (MQ=255)
                                                   tgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAACATGAACATTGAACGCTGAATg                                                                                    >  1:54604/1‑141 (MQ=255)
                                                   tgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATg                                                                                    >  1:146280/1‑141 (MQ=255)
                                                   tgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATg                                                                                    <  1:452857/141‑1 (MQ=255)
                                                      ggCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAACACTGAATGACa                                                                                 >  1:76835/1‑141 (MQ=255)
                                                       cgttgctaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACaa                                                                                >  1:157388‑M1/58‑141 (MQ=255)
                                                       cgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACaa                                                                                >  1:153570‑M1/58‑141 (MQ=255)
                                                        gttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGacatt                                                                               <  2:448727‑M1/85‑3 (MQ=255)
                                                        gttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAat                                                                               <  1:452060‑M1/85‑1 (MQ=255)
                                                          tgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATTTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAatat                                                                             >  1:191965‑M1/55‑141 (MQ=255)
                                                          tgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAatat                                                                             >  1:154580‑M1/55‑141 (MQ=255)
                                                             aaagcaaATCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTc                                                                          <  2:440946/141‑1 (MQ=255)
                                                                     tcaactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAACTGAACATTGAAAACTGAATGACAATATGTCAACGTTaa                                                                  >  2:205197/1‑141 (MQ=255)
                                                                      ttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATAATTGTCGGTAAGAAAAATCAACATTGAAAACTGAATGAAAATATGTCAACGTTAAt                                                                 <  2:327502‑M1/99‑1 (MQ=255)
                                                                      caactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCCACGTTAAt                                                                 >  1:154432/1‑141 (MQ=255)
                                                                      caactTAAGAAATAACTTGAATTACTAACGAAAAATAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAt                                                                 <  1:311729/141‑1 (MQ=255)
                                                                          ttAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCa                                                             >  1:70914/1‑141 (MQ=255)
                                                                           tAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaa                                                            <  1:7834/141‑1 (MQ=255)
                                                                           atcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaa                                                            <  1:73940‑M1/104‑1 (MQ=255)
                                                                            tcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaa                                                           >  1:277615‑M1/37‑141 (MQ=255)
                                                                            aaaaagtgtggttaaaaaataaaatggagaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaa                                                           >  2:26361‑M1/37‑141 (MQ=255)
                                                                             aGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaa                                                          >  2:344489/1‑141 (MQ=255)
                                                                             aGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaa                                                          <  1:126291/141‑1 (MQ=255)
                                                                               aaaTAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTCCTCGTTTTTTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaaaa                                                        <  1:234629/141‑1 (MQ=25)
                                                                                 ttttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACg                                                      >  1:201883‑M1/32‑141 (MQ=255)
                                                                                 ttttggtgaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACg                                                      >  1:280862‑M1/32‑141 (MQ=255)
                                                                                           aaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAAGTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAg                                            >  1:150469‑M1/22‑141 (MQ=255)
                                                                                             aCTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTTTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGtt                                          >  2:83488/1‑141 (MQ=255)
                                                                                             aCTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGtt                                          >  2:114301/1‑141 (MQ=255)
                                                                                              cTAACGAAAATTAATTTTTTAGAGGAACTGTTTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTa                                         >  2:37108/1‑141 (MQ=21)
                                                                                              cTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTa                                         >  2:215146/1‑141 (MQ=255)
                                                                                              cTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTa                                         >  2:410424/1‑141 (MQ=255)
                                                                                                aaCGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACa                                       <  2:146280/141‑1 (MQ=255)
                                                                                                aaCGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACa                                       >  2:287275/1‑141 (MQ=255)
                                                                                                 aCGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAATTGAATGACAATATGTCAACGTTAATTCCAAAAAACGAAAATATACGTTCCaa                                      >  2:48750/1‑141 (MQ=37)
                                                                                                   atggacaagataaAAAAAGTTATTGACTTACATGTTGATAAAATTTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAAc                                    <  1:343823‑M1/128‑1 (MQ=255)
                                                                                                     ggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACat                                  >  1:449731‑M1/12‑141 (MQ=255)
                                                                                                     ggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACat                                  >  2:427466‑M1/12‑141 (MQ=255)
                                                                                                           attTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTa                            <  1:315331/141‑1 (MQ=255)
                                                                                                            ttttAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTATTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCAAAAAAACGTAACTATAAGTTAAAAACATTATTTAg                           >  2:152986/1‑141 (MQ=255)
                                                                                                              ttAAATAGTTATTGACTTAAATGTTAATAAAAAGTATAAATAATTATTGTCAGTAAAAAAAATGAACATTGAAAACTGAATGATAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTa                         <  2:17304/141‑1 (MQ=21)
                                                                                                              aaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTACACAGAACTTAACTATAAGTTACAAACATTATTTAGTa                         >  1:96785‑M1/3‑141 (MQ=255)
                                                                                                               tAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTAt                        <  1:288411/141‑1 (MQ=255)
                                                                                                               aAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTAt                        >  1:251269‑M1/2‑141 (MQ=255)
                                                                                                               aAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTAt                        >  1:275889‑M1/2‑141 (MQ=255)
                                                                                                                  aaaGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATGTAGTAGTTa                     >  1:254304/1‑141 (MQ=255)
                                                                                                                   aaataaTTGAATGAAATGTTAATAAAATATATAATTAATTCTTGTCAGTAAGAAAAATAAACATTGAAAACTGAATAACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTAt                    <  2:233265/136‑1 (MQ=18)
                                                                                                                     gTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAATTTACAAACATTATATAGTATTTATGa                  >  1:207533/1‑141 (MQ=255)
                                                                                                                      ttATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAg                 <  1:487528/141‑1 (MQ=255)
                                                                                                                        agtGACTTAAATGTTAATAAAATGTATAATTAATGCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCt               <  2:281965/139‑1 (MQ=255)
                                                                                                                               tAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTACAAAAAACGTAACTATAAGTTACACACATTATTTAGTATTTATGAGCTAATCaaa        >  2:74625/1‑141 (MQ=255)
                                                                                                                                   tGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcatc    <  2:198835/141‑1 (MQ=255)
                                                                                                                                   tGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAACCAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcatc    <  1:446716/141‑1 (MQ=255)
                                                                                                                                   tGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAACTCCCAATAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcatc    <  1:446630/141‑1 (MQ=255)
                                                                                                                                     ttAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcatcat  <  2:418785/141‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                          
TAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCAT  >  USA300TCH1516_ALE/512508‑512781

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