Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,448,437 G→A 5.4% intergenic (‑14/+15) USA300HOU_RS07195 ← / ← mhqA_2 hypothetical protein/Putative ring‑cleaving dioxygenase MhqA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,448,4370GA 5.4% 284.9 / 3.5 92intergenic (‑14/+15)USA300HOU_RS07195/mhqA_2hypothetical protein/Putative ring‑cleaving dioxygenase MhqA
Reads supporting (aligned to +/- strand):  ref base G (44/43);  new base A (2/3);  total (46/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.76e-01

TCGTTGATAACATCATATTGTTTGCCAACCGTTAACATGTCAGAAACTGTAAATTTCTTTGCATCTGTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCTCAAGCGTTTGATATGGTATTTCCACTTCATCTAAGCGCGATTTCGCTTGAT  >  USA300TCH1516_ALE/1448299‑1448576
                                                                                                                                          |                                                                                                                                           
tCGTTGATAACATCATATTGTTTGCCAACCGTTAACATGTCAGAAACTGTAAATTTCTTTGCATCTGTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCt                                                                                                                                           <  1:817051/141‑1 (MQ=255)
  gTTGATAACATCATATTGTTTGCCAACCGTTAACATGTCAGAAACTGTAAATTTCTTTGCATCTGTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCACTAt                                                                                                                                         >  3:375735/1‑141 (MQ=255)
   ttGATAACATCATATTGTTTGCCAACCGTTAACATGTCAGAAACTGTAAATTTCTTTGCATCTGTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATc                                                                                                                                        >  2:954237/1‑141 (MQ=255)
   ttGATAACATCATATTGTTTGCCAACCGTTAACATGTCAGAAACTGTAAATTTCTTTGCATCTGTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATc                                                                                                                                        <  1:853009/141‑1 (MQ=255)
   ttGATAACATCATATTGTTTGCCAACCGTTAACATGTCAGAAACTGTAAATTTCTTTGCATCTGTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATc                                                                                                                                        <  1:440012/141‑1 (MQ=255)
    tGATAACATCATATTGTTTGCCAACCGTTAACATGTCAGAAACTGTAAATTTCTTTGCATCTGTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCg                                                                                                                                       <  3:77015/141‑1 (MQ=255)
           atcatATTGTTTGCCAACCGTTAACATGTCAGAAACTGTAAATTTCTTTGCATCTGTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTaaaaaa                                                                                                                                <  3:16954/141‑1 (MQ=255)
           atcatATTGTTTGCCAACCGTTAACATGTCAGAAACTGTAAATTTCTTTGCATCTGTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTaaaaaa                                                                                                                                >  1:761476/1‑141 (MQ=255)
             catATTGTTTGCCAACCGTTAACATGTCAGAAACTGTAAATTTCTTTGCATCTGTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAat                                                                                                                              <  4:452081/141‑1 (MQ=255)
                tttGTTTGCCAACCGTTAACATGTCAGAAACTGTAAATTTCTTGGCATGTGTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTTCAAGTTTCATGTAACATAAAGTTTGTCCTCGCTATCGTAAAAAAattat                                                                                                                           <  2:732796/140‑1 (MQ=255)
                 ttGTTTGCCAACCGTTAACATGTCAGAAACTGTAAATTTCTTTGCATCTGTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAattatt                                                                                                                          >  1:1015283/1‑141 (MQ=255)
                   gTTTGCCAACCGTTAACATGTCAGAAACTGTAAATTTCTTTGCATCTGTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTAtttt                                                                                                                        >  2:885172/1‑141 (MQ=255)
                    tttGCCAACCGTTAACATGTCAGAAACTGTAAATTTCTTTGCATCTGTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTAttttt                                                                                                                       <  1:857863/141‑1 (MQ=255)
                    tttGCCAACCGTTAACATGTCAGAAACTGTAAATTTCTTTGCATCTGTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTAttttt                                                                                                                       >  2:349153/1‑141 (MQ=255)
                     ttGCCAACCGTTAACATGCCAGAAACTGTAAATTTCTTTGCATCTGTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTAatttgg                                                                                                                      >  1:226951/1‑139 (MQ=255)
                        ccAACCGTTAACATGTCAGAAACTGTAAATTTCTTTGCATCTGTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTaa                                                                                                                   >  1:518555/1‑141 (MQ=255)
                          aaCCGTTAACATGTCAGAAACTGTAAATTTCTTTGCATCTGTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAg                                                                                                                 <  4:297579/141‑1 (MQ=255)
                                tAACATGTCAGAAACTGTAAATTTCTTTGCATCTGTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAata                                                                                                           >  1:589904/1‑141 (MQ=255)
                                tAACATGTCAGAAACTGTAAATTTCTTTGCATCTGTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAata                                                                                                           <  3:169522/141‑1 (MQ=255)
                                  aCATGTCAGAAACTGTAAATTTCTTTGCATCTGTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAACTTATTTTTGTAAAGTGAAtata                                                                                                         >  2:611384/1‑141 (MQ=255)
                                  aCATGTCAGAAACTGTAAATTTCTTTGCATCTGTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAAtata                                                                                                         <  1:168592/141‑1 (MQ=255)
                                     tGTCAGAAACTGTAAACTCCTTTGCATCTGTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAAc                                                                                                      <  1:332841/141‑1 (MQ=255)
                                        cataaaCTGTAAATTTCTTTGCATCTGTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGAt                                                                                                   <  3:209752/138‑1 (MQ=255)
                                          gAAACTGTAAATTTCTTTGCATCTGTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATaa                                                                                                 <  1:248729/141‑1 (MQ=255)
                                          gAAACTGTAAATTTCTTTGCATCTGTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATaa                                                                                                 <  2:755935/141‑1 (MQ=255)
                                                gTAAATTTCTTTGCATCTGTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACt                                                                                           >  1:613292/1‑141 (MQ=255)
                                                 tAAATTTCTTTGCATCTGTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTg                                                                                          >  2:819046/1‑141 (MQ=255)
                                                  aaaTTTCTTTGCATCTGTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGt                                                                                         >  2:46718/1‑141 (MQ=255)
                                                  aaaTTTCTTTGCATCTGTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAACCCACTGt                                                                                         >  4:224937/1‑141 (MQ=255)
                                                   aaTTTCTTTGCATCTGTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGtt                                                                                        >  2:9669/1‑141 (MQ=255)
                                                       tCTTTGCATCTGTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCa                                                                                    >  2:619683/1‑141 (MQ=255)
                                                       tCTTTGCATCTGTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCa                                                                                    >  2:619775/1‑141 (MQ=255)
                                                        cTTTGCATCTGTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAg                                                                                   >  3:271097/1‑141 (MQ=255)
                                                        cTTTGCATCTGTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAg                                                                                   <  3:4103/141‑1 (MQ=255)
                                                        cTTTGCATCTGTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAg                                                                                   <  3:380780/141‑1 (MQ=255)
                                                        cTTTGCATCGGTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAg                                                                                   >  1:955028/1‑141 (MQ=255)
                                                          ttGCATCTGTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTa                                                                                 <  2:814535/141‑1 (MQ=255)
                                                                  gTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAAt                                                                         <  1:790285/141‑1 (MQ=255)
                                                                  gTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAAt                                                                         <  1:982970/141‑1 (MQ=255)
                                                                  gTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAAt                                                                         >  1:229240/1‑141 (MQ=255)
                                                                    aTTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTg                                                                       >  1:619113/1‑141 (MQ=255)
                                                                     ttAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGa                                                                      <  1:796931/141‑1 (MQ=255)
                                                                      tAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGAt                                                                     <  4:453174/141‑1 (MQ=255)
                                                                        aTGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTc                                                                   >  3:411674/1‑141 (MQ=255)
                                                                           aCTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAat                                                                <  2:824083/141‑1 (MQ=255)
                                                                                tACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATatcat                                                           >  3:358265/1‑141 (MQ=255)
                                                                                tACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATatcat                                                           >  4:91387/1‑141 (MQ=255)
                                                                                tACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCACTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATatcat                                                           >  3:61795/1‑141 (MQ=255)
                                                                                 aCTTTTTTGATTTGTGTACTTTCTTTCCAATTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATatcatc                                                          <  2:1035263/141‑1 (MQ=255)
                                                                                     tCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGc                                                      <  1:1020343/141‑1 (MQ=255)
                                                                                     tCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGc                                                      >  4:124176/1‑141 (MQ=255)
                                                                                      cTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCt                                                     <  2:761476/141‑1 (MQ=255)
                                                                                      cTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCt                                                     >  2:475605/1‑141 (MQ=255)
                                                                                         gATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCTCaa                                                  >  4:185830/1‑141 (MQ=255)
                                                                                            tGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCTTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCTCAAGCg                                               <  4:425482/141‑1 (MQ=255)
                                                                                            tGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCTCAAGCg                                               >  2:157587/1‑141 (MQ=255)
                                                                                             ggTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTTTTTTTGCTATCGTTTAAAAATTGTCTTTGTAAAGTGAATATAAACGTTAATCCACTGTTTTCGGTAATTCTAATTTTTTCAATATCATCTTTCTCAAGCGt                                              >  2:446547/1‑141 (MQ=255)
                                                                                             ggTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCTCAAGCGt                                              >  1:92265/1‑141 (MQ=255)
                                                                                              gtgtACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCTCAAGCGtt                                             >  2:530100/1‑141 (MQ=255)
                                                                                               tgtACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCTCAAGCGttt                                            <  1:173061/141‑1 (MQ=255)
                                                                                                gtACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCTCAAGCGTTTg                                           >  2:299970/1‑141 (MQ=255)
                                                                                                gtACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCTCAAGCGTTTg                                           <  2:119190/141‑1 (MQ=255)
                                                                                                 tACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCTCAAGCGTTTGa                                          >  1:91491/1‑141 (MQ=255)
                                                                                                  aCGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCTCAAGCGTTTGat                                         >  4:291543/1‑141 (MQ=255)
                                                                                                          tCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCTCAAGCGTTTGATATGGTAtt                                 >  2:204339/1‑141 (MQ=255)
                                                                                                           ccAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCTCAAGCGTTTGATATGGTAttt                                <  1:788897/141‑1 (MQ=255)
                                                                                                           ccAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCTCAAGCGTTTGATATGGTAttt                                <  4:178835/141‑1 (MQ=255)
                                                                                                              aGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCTCAAGCGTTTGATATGGTATTTCCa                             >  2:531102/1‑141 (MQ=255)
                                                                                                              aGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCTCAAGCGTTTGATATGGTATTTCCa                             >  1:453660/1‑141 (MQ=255)
                                                                                                                 ttCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCTCAAGCGTTTGATATGGTATTTCCACtt                          <  1:1022011/141‑1 (MQ=255)
                                                                                                                 ttCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCTCAAGCGTTTGATATGGTATTTCCACtt                          <  3:121750/141‑1 (MQ=255)
                                                                                                                  tCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTATAGTTTATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCTCATGCGTTTGATATGGTATTTCCCCttt                         >  2:677823/1‑140 (MQ=255)
                                                                                                                      gTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCTCAAGCGTTTGATATGGTATTTCCACTTCATCt                     >  3:110928/1‑141 (MQ=255)
                                                                                                                           aTAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCTCAAGCGTTTGATATGGTATTTCCACTTCATCTAAgcg                <  2:548140/141‑1 (MQ=255)
                                                                                                                            tAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCCTCTTGCTCAAGCGTTTGATATGGTATTTCCACTTCATCTAAgcgc               <  1:792147/141‑1 (MQ=255)
                                                                                                                             aaaTTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCTCAAGCGTTTGATATGGTATTTCCACTTCATCTAAGcgcg              <  1:736659/141‑1 (MQ=255)
                                                                                                                               aTTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCTCAAGCGTTTGATATGGTATTTCCACTTCATCTAAGCGCGAt            <  1:678411/141‑1 (MQ=255)
                                                                                                                               aTTTTGTCCTCACTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCTCAAGCGTTTGATATGGTATTTCCACTTCATCTAAGCGCGAt            <  2:687428/141‑1 (MQ=255)
                                                                                                                                ttttGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCTCAAGCGTTTGATATGGTATTTCCACTTCATCTAAGCGCGAtt           <  4:52555/141‑1 (MQ=255)
                                                                                                                                 tttGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGATCAAGCGTTTGATATGTTATTTCTACTTCATCTAAGCAGGAttt          >  2:564905/1‑141 (MQ=255)
                                                                                                                                   tGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCTCAAGCGTTTGATATGGTATTTCCACTTCATCTAAGCGCGATTTCg        <  1:161808/141‑1 (MQ=255)
                                                                                                                                    gTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCTCAAGCGTTTGATATGGTATTTCCACTTCATCTAAGCTCGATTTCGc       >  2:566731/1‑141 (MQ=255)
                                                                                                                                    gTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCTCAAGCGTTTGATATGGTATTTCCACTTCATCTAAGCGCGATTTCGc       <  2:521118/141‑1 (MQ=255)
                                                                                                                                    gTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCTCAAGCGTTTGATATGGTATTTCCACTTCATCTAAGAGCGATTTTGc       >  2:762437/1‑141 (MQ=255)
                                                                                                                                    gTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCTCAAGCGTTTGATATGGTATTTCCACTTCATATAAGCGCGATTTCGc       <  1:404325/141‑1 (MQ=255)
                                                                                                                                    gTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCTCAAGCGTTTGATATGGTATTCCCACTTCATCTAAGCGCGATTTCGc       <  3:380957/141‑1 (MQ=255)
                                                                                                                                    gTCCTCACTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCTCAAGCGTTTGATATGGTATTTCCACTTCATCTAAGCGCGCTTTCGc       <  4:375735/141‑1 (MQ=255)
                                                                                                                                     tCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCTCAAGCGTTTGATATGGTATTTCCACTTCATCTAAGCGCGATTTCGCt      >  1:86769/1‑141 (MQ=255)
                                                                                                                                     tCCTCACTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCTCAAGCGTTTGATATGGTATTTCCACTTCATCTAAGCGCGATTTCGCt      <  1:505560/141‑1 (MQ=255)
                                                                                                                                      ccTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCTCAAGCGTTTGATATGGTATTTCCACTTCATCTAAGCGCGATTTCGCtt     >  1:619830/1‑141 (MQ=255)
                                                                                                                                         cGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCTCAAGCGTTTGATATGGTATTTCCACTTCATCTAAGCGCGATTTCGCTTGAt  >  1:493578/1‑141 (MQ=255)
                                                                                                                                         cGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCTCAAGCGTTTGATATGGTATTTCCACTTCATCTAAGCGCGATTTCGCTTGAt  <  3:77293/141‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                           
TCGTTGATAACATCATATTGTTTGCCAACCGTTAACATGTCAGAAACTGTAAATTTCTTTGCATCTGTATTAATGACTTCTACTTTCTTGATTGGTGTACGTTCTTTCCAAGTTTCATGTAACATAAATTTTGTCCTCGCTATCGTAAAAAAATTATTTTTGTAAAGTGAATATAAACGATAATCCACTGTTTTCAGTAATTCTAATTGATTCAATATCATCTTGCTCAAGCGTTTGATATGGTATTTCCACTTCATCTAAGCGCGATTTCGCTTGAT  >  USA300TCH1516_ALE/1448299‑1448576

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