Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,526,934 C→A 10.2% intergenic (+195/‑189) cysK → / → ptsH cysteine synthase A, O‑acetylserine sulfhydrolase A subunit/phosphohistidinoprotein‑hexose phosphotransferase component of PTS system (Hpr)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,526,9340CA10.2% 1091.3 / 60.1 329intergenic (+195/‑189)cysK/ptsHcysteine synthase A, O‑acetylserine sulfhydrolase A subunit/phosphohistidinoprotein‑hexose phosphotransferase component of PTS system (Hpr)
Reads supporting (aligned to +/- strand):  ref base C (158/137);  new base A (15/19);  total (173/156)
Fisher's exact test for biased strand distribution p-value = 3.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 336 total aligned reads displayed.

TTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACT  >  CP009273/2526786‑2527082
                                                                                                                                                    |                                                                                                                                                    
     tgGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCtt                                                                                                                                                >  2:3603933/1‑150 (MQ=255)
     tgGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCtt                                                                                                                                                <  1:37640/150‑1 (MQ=255)
     tgGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCtt                                                                                                                                                >  2:2745239/1‑150 (MQ=255)
        tcTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAg                                                                                                                                             >  2:727434/2‑150 (MQ=255)
                aaaCTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATaa                                                                                                                                     >  1:1695983/1‑150 (MQ=255)
                 aaCTTTCGCCCCTCCTGTCTTTGATTCAGCCTGTCGGAACTGTTATTTAAACAGACTAATTATTTTGATGCTCGAAATTAATCGTTACAGGAAAATACAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCGTAAt                                                                                                                                    <  1:3020336/150‑1 (MQ=255)
                  aCTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATg                                                                                                                                   <  1:2564139/150‑1 (MQ=255)
                  aCTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATg                                                                                                                                   <  1:1920273/150‑1 (MQ=255)
                  aCTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTATTCCTTTAGCGGCATAATg                                                                                                                                   >  2:3354443/1‑150 (MQ=255)
                    tttCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGtt                                                                                                                                 >  1:4809744/1‑150 (MQ=255)
                     ttCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGttt                                                                                                                                <  2:5009702/150‑1 (MQ=255)
                     ttCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTATTCCTTTAGCGGCATAATGttt                                                                                                                                >  2:1512948/1‑150 (MQ=255)
                        gCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAAt                                                                                                                             <  2:4928268/150‑1 (MQ=255)
                        gCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAAt                                                                                                                             <  1:3684123/150‑1 (MQ=255)
                        gCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAAt                                                                                                                             <  1:100532/150‑1 (MQ=255)
                           cctcctGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGAc                                                                                                                          <  1:953141/150‑1 (MQ=255)
                              cctGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTa                                                                                                                       <  2:845021/150‑1 (MQ=255)
                                tGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACg                                                                                                                     >  2:3619686/1‑150 (MQ=255)
                                 ggCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGa                                                                                                                    >  1:30738/1‑150 (MQ=255)
                                  gCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGaa                                                                                                                   >  1:4710201/1‑150 (MQ=255)
                                    aTTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAAc                                                                                                                 >  2:4459530/1‑150 (MQ=255)
                                    aTTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAAc                                                                                                                 <  2:1359930/150‑1 (MQ=255)
                                    aTTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAAc                                                                                                                 <  2:2234351/150‑1 (MQ=255)
                                    aTTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTATTCCTTTAGCGGCATAATGTTTAATGACGTACGAAAc                                                                                                                 <  2:3961780/150‑1 (MQ=255)
                                      tGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGt                                                                                                               >  1:373083/1‑150 (MQ=255)
                                         ttCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAg                                                                                                            <  1:4427049/150‑1 (MQ=255)
                                          tCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGc                                                                                                           >  2:4665857/1‑150 (MQ=255)
                                              ccTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTc                                                                                                       <  2:165227/150‑1 (MQ=255)
                                               cTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCa                                                                                                      >  1:2295868/1‑150 (MQ=255)
                                               cTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTATTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCa                                                                                                      <  1:657040/150‑1 (MQ=255)
                                                tGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCaa                                                                                                     >  2:3428765/1‑150 (MQ=255)
                                                   cGGAACTGGTATTTAACCCGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAacac                                                                                                  >  1:1553684/1‑150 (MQ=255)
                                                      aaCTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCg                                                                                               <  2:1583576/150‑1 (MQ=255)
                                                        cTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGcc                                                                                             <  1:4476392/150‑1 (MQ=255)
                                                         tGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCa                                                                                            >  1:1747186/1‑150 (MQ=255)
                                                         tGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCa                                                                                            <  1:1867002/150‑1 (MQ=255)
                                                              tttAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAAt                                                                                       <  1:1401900/150‑1 (MQ=255)
                                                                tAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATgg                                                                                     >  2:2557941/1‑150 (MQ=255)
                                                                tAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATgg                                                                                     <  1:3222500/150‑1 (MQ=255)
                                                                  aCCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTATTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGAc                                                                                   <  1:4968186/150‑1 (MQ=255)
                                                                   ccAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACt                                                                                  >  2:3060204/1‑150 (MQ=255)
                                                                   ccAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACt                                                                                  <  2:4396816/150‑1 (MQ=255)
                                                                     aGACTAATTATTTTGATGCGCGAAATTGATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGt                                                                                >  2:2807988/1‑150 (MQ=255)
                                                                     aGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGt                                                                                >  2:1971372/1‑150 (MQ=255)
                                                                      gACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTa                                                                               >  2:1226519/1‑150 (MQ=255)
                                                                        cTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTAtt                                                                             >  2:2140392/1‑150 (MQ=255)
                                                                         tAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTg                                                                            <  1:4336336/150‑1 (MQ=255)
                                                                           attattTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTgcg                                                                          >  2:2507773/1‑150 (MQ=255)
                                                                             tattTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGct                                                                        >  1:712451/1‑150 (MQ=255)
                                                                             tattTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGct                                                                        <  1:2445278/150‑1 (MQ=255)
                                                                               ttttGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGctct                                                                      >  1:836752/1‑150 (MQ=255)
                                                                                tttGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCtt                                                                     <  2:1046501/150‑1 (MQ=255)
                                                                                tttGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCtt                                                                     <  1:2520720/150‑1 (MQ=255)
                                                                                   gATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGt                                                                  >  2:2926171/1‑150 (MQ=255)
                                                                                   gATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGt                                                                  >  2:2412726/1‑150 (MQ=255)
                                                                                   gATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGt                                                                  >  2:2197224/1‑150 (MQ=255)
                                                                                      gcgcGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTgcg                                                               <  2:3445454/150‑1 (MQ=255)
                                                                                       cgcgAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTgcgt                                                              >  2:4050400/1‑150 (MQ=255)
                                                                                         cgAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTAAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTgcgtcg                                                            >  1:2233183/1‑150 (MQ=255)
                                                                                           aaaTTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTATTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCCACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTgcgtcgcg                                                          >  1:1371067/1‑150 (MQ=255)
                                                                                            aaTTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTATTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTgcgtcgcgt                                                         >  1:1140668/1‑150 (MQ=255)
                                                                                              ttAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCt                                                       <  2:4221904/150‑1 (MQ=255)
                                                                                                 aTCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGtt                                                    >  2:4198565/1‑150 (MQ=255)
                                                                                                  tCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTa                                                   <  1:5059235/150‑1 (MQ=255)
                                                                                                    gTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTaaa                                                 >  1:2918210/1‑150 (MQ=255)
                                                                                                        cAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACt                                             <  1:3892510/150‑1 (MQ=255)
                                                                                                         aGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTg                                            >  2:553408/1‑150 (MQ=255)
                                                                                                         aGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTg                                            >  1:3373007/1‑150 (MQ=255)
                                                                                                          ggAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTATTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTgg                                           <  1:4866411/150‑1 (MQ=255)
                                                                                                           gAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGgc                                          >  1:5141329/1‑150 (MQ=255)
                                                                                                            aaaaGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGgcg                                         >  2:3857993/1‑150 (MQ=255)
                                                                                                            aaaaGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGgcg                                         >  1:2565024/1‑150 (MQ=255)
                                                                                                            aaaaGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGgcg                                         <  1:4058946/150‑1 (MQ=255)
                                                                                                               aGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTa                                      <  1:4665626/150‑1 (MQ=255)
                                                                                                                 ccAAAGCTGAATCGATTTTATGATTTGGTTCAATTATTCCATTAGCGGCATATTGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAAc                                    <  2:4180155/150‑1 (MQ=255)
                                                                                                                      gCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATAc                               <  2:4766040/150‑1 (MQ=255)
                                                                                                                        tGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAg                             <  1:311940/150‑1 (MQ=255)
                                                                                                                        tGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAg                             >  1:2506192/1‑150 (MQ=255)
                                                                                                                        tGAATCGATTTTATGATTTGGTTCAATTATTCCTTTAGCTGCATTATGTTTAATGACGTACGAAAGGTCAGCTGTCTACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACTATTCAg                             <  2:1470398/150‑1 (MQ=255)
                                                                                                                          aaTCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACTAAACGTCAGCGGACAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACAGGCGCTAACAAAACATGc                           >  1:1418799/1‑150 (MQ=255)
                                                                                                                           aTCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCt                          >  1:3816802/1‑150 (MQ=255)
                                                                                                                              gATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTaaa                       >  1:5156658/1‑150 (MQ=255)
                                                                                                                              gATTTTATGATGTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACACTACAGGCTcaa                       >  2:2031188/1‑150 (MQ=255)
                                                                                                                               aTTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAg                      <  1:1211104/150‑1 (MQ=255)
                                                                                                                               aTTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAg                      >  2:1763828/1‑150 (MQ=255)
                                                                                                                                  ttATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCg                   >  2:3245498/1‑150 (MQ=255)
                                                                                                                                  ttATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCg                   <  2:4317962/150‑1 (MQ=255)
                                                                                                                                   tATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGa                  >  1:3624251/1‑150 (MQ=255)
                                                                                                                                    aTGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGaa                 <  1:2743484/150‑1 (MQ=255)
                                                                                                                                      gATTAGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAAcc               <  2:4809708/150‑1 (MQ=255)
                                                                                                                                       aTTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCg              <  2:4537959/150‑1 (MQ=255)
                                                                                                                                          tGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCa           <  2:522954/150‑1 (MQ=255)
                                                                                                                                          tGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCa           <  1:3965694/150‑1 (MQ=255)
                                                                                                                                            gTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAgg         >  1:2054751/1‑150 (MQ=255)
                                                                                                                                            gTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAgg         <  2:5268191/150‑1 (MQ=255)
                                                                                                                                            gTTCAATTATTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAgg         <  2:5066219/150‑1 (MQ=255)
                                                                                                                                               cAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTa      <  2:3710591/150‑1 (MQ=255)
                                                                                                                                                aaTTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAg     >  2:1894060/1‑150 (MQ=255)
                                                                                                                                                 aTTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGa    >  1:4316504/1‑150 (MQ=255)
                                                                                                                                                   tcttcCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACt  <  2:5112406/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGGCATAATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACT  >  CP009273/2526786‑2527082

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