Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,860,183 G→T 100% R183S (CGC→AGC)  mipA ← scaffolding protein for murein synthesizing machinery

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,860,1830GT96.8% 754.4 / 14.0 252R183S (CGC→AGC) mipAscaffolding protein for murein synthesizing machinery
Reads supporting (aligned to +/- strand):  ref base G (0/8);  new base T (106/138);  total (106/146)
Fisher's exact test for biased strand distribution p-value = 2.24e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.40e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 253 total aligned reads displayed.

CATCGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCGAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAGGCCATATCCCAGACGATGCCGTTGCTGTTA  >  CP009273/1860037‑1860328
                                                                                                                                                  |                                                                                                                                                 
caTCGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAg                                                                                                                                                 >  1:3907962/1‑149 (MQ=255)
 aTCGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGc                                                                                                                                                <  2:4080528/149‑1 (MQ=255)
         gTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCtt                                                                                                                                        >  2:2056954/1‑149 (MQ=255)
           cAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTg                                                                                                                                      <  2:2313104/149‑1 (MQ=255)
            aGTAACTTCCTCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTgc                                                                                                                                     <  2:3790897/149‑1 (MQ=255)
            aGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTgc                                                                                                                                     >  1:1063694/1‑149 (MQ=255)
            aGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTgc                                                                                                                                     <  2:693612/149‑1 (MQ=255)
            aGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTgc                                                                                                                                     <  2:4043894/149‑1 (MQ=255)
            aGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTgc                                                                                                                                     >  2:2024609/1‑149 (MQ=255)
                aCTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTTCGCGa                                                                                                                                 >  1:628968/1‑149 (MQ=255)
                aCTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGa                                                                                                                                 >  2:1720728/1‑149 (MQ=255)
                 cTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGAt                                                                                                                                <  2:3219820/149‑1 (MQ=255)
                   tcatcaGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATAc                                                                                                                              >  1:3025032/1‑149 (MQ=255)
                    catcaGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACg                                                                                                                             >  1:1184036/1‑149 (MQ=255)
                      tcaGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGcc                                                                                                                           <  2:1386978/149‑1 (MQ=255)
                           cagacGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCataat                                                                                                                      <  2:1488731/149‑1 (MQ=255)
                           cagacGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCataat                                                                                                                      <  1:1050244/149‑1 (MQ=255)
                             gacGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAg                                                                                                                    <  2:2687890/149‑1 (MQ=255)
                             gacGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAg                                                                                                                    <  1:308459/149‑1 (MQ=255)
                              acGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGt                                                                                                                   <  1:1631259/149‑1 (MQ=255)
                               cGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTa                                                                                                                  <  1:2056999/149‑1 (MQ=255)
                                gggTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTAt                                                                                                                 <  2:1691593/149‑1 (MQ=255)
                                gggTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTAt                                                                                                                 <  1:3779579/149‑1 (MQ=255)
                                gggTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTAt                                                                                                                 <  1:3445607/149‑1 (MQ=255)
                                  gtgtAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTc                                                                                                               >  1:960807/1‑149 (MQ=255)
                                  gtgtAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCAGAATAGTATTc                                                                                                               >  1:31049/1‑149 (MQ=255)
                                    gtAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGt                                                                                                             >  2:1249232/1‑149 (MQ=255)
                                    gtAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGt                                                                                                             <  1:1029535/149‑1 (MQ=255)
                                       gCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCt                                                                                                          >  2:1319794/1‑149 (MQ=255)
                                         gTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTgg                                                                                                        >  1:2056317/1‑149 (MQ=255)
                                          tGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGt                                                                                                       <  1:2294476/149‑1 (MQ=255)
                                           gCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGtt                                                                                                      <  1:182196/149‑1 (MQ=255)
                                             ggTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGtttt                                                                                                    >  2:1400216/1‑149 (MQ=255)
                                             ggTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGtttt                                                                                                    >  1:3993071/1‑149 (MQ=255)
                                               tACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCg                                                                                                  >  1:1163453/1‑149 (MQ=255)
                                                aCCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGc                                                                                                 <  1:1827535/149‑1 (MQ=255)
                                                aCCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGc                                                                                                 <  1:3116411/149‑1 (MQ=255)
                                                aCCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGc                                                                                                 >  2:135322/1‑149 (MQ=255)
                                                  cGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGCCTCTTTGCGCGATACGCCCTAATAGTATTCGTTCTGGTTTTAGCTg                                                                                               >  1:2768271/1‑149 (MQ=255)
                                                    tAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGtt                                                                                             <  1:1578692/149‑1 (MQ=255)
                                                     aaaCACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTc                                                                                            <  2:3754204/149‑1 (MQ=255)
                                                     aaaCACTCCACTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTc                                                                                            <  2:2211660/149‑1 (MQ=255)
                                                      aaCACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTcc                                                                                           >  1:2360304/1‑149 (MQ=255)
                                                      aaCACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCGAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTcc                                                                                           <  1:3425700/149‑1 (MQ=255)
                                                         acTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACt                                                                                        <  1:3041132/149‑1 (MQ=255)
                                                          cTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTg                                                                                       >  2:649568/1‑149 (MQ=255)
                                                             cAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGcac                                                                                    <  1:996245/149‑1 (MQ=255)
                                                              aGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGcaca                                                                                   <  2:4014421/149‑1 (MQ=255)
                                                              aGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTAGCTGTTACACTGcaca                                                                                   >  1:69117/1‑149 (MQ=255)
                                                                  gCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCaa                                                                               <  1:592491/149‑1 (MQ=255)
                                                                  gCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCGAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCaa                                                                               <  1:479934/149‑1 (MQ=255)
                                                                    cGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATa                                                                             >  2:1202876/1‑149 (MQ=255)
                                                                         aaGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCgg                                                                        <  1:811276/149‑1 (MQ=255)
                                                                          aGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGa                                                                       <  2:514958/149‑1 (MQ=255)
                                                                           gTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAg                                                                      <  1:2931424/149‑1 (MQ=255)
                                                                           gTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAg                                                                      >  1:3699279/1‑149 (MQ=255)
                                                                           gTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAg                                                                      <  1:258227/149‑1 (MQ=255)
                                                                           gTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAg                                                                      >  1:3735530/1‑149 (MQ=255)
                                                                            ttGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCGGCTAGCGGACTCTTTGCGCGCTACGCCATAATAGTATTAGTTATGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGt                                                                     >  1:835986/1‑149 (MQ=255)
                                                                               tAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCAc                                                                  <  2:3523911/149‑1 (MQ=255)
                                                                               tAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCAc                                                                  <  1:1915407/149‑1 (MQ=255)
                                                                                  cTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTAGAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACGGCACACCAATACCCGGAGTCACGGt                                                               <  2:4078738/149‑1 (MQ=255)
                                                                                    ggCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCa                                                             <  2:1212493/149‑1 (MQ=255)
                                                                                         tCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCa                                                        <  2:351737/149‑1 (MQ=255)
                                                                                         tCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCa                                                        <  2:3337446/149‑1 (MQ=255)
                                                                                         tCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATAACCGGAGTCACGGTCAGGCCa                                                        >  1:3277554/1‑149 (MQ=255)
                                                                                            gCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCg                                                     >  2:2942106/1‑149 (MQ=255)
                                                                                                 aGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGgt                                                <  1:1231754/149‑1 (MQ=255)
                                                                                                 aGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGgt                                                >  2:2689065/1‑149 (MQ=255)
                                                                                                 aGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTAGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGgt                                                >  1:2136098/1‑149 (MQ=255)
                                                                                                  ggTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGgtg                                               <  2:3092852/149‑1 (MQ=255)
                                                                                                        ggCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGa                                         >  1:2590791/1‑149 (MQ=255)
                                                                                                          cTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATa                                       <  1:3065961/149‑1 (MQ=255)
                                                                                                             cAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACaa                                    <  1:2812773/149‑1 (MQ=255)
                                                                                                               gCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAAcc                                  >  2:3615592/1‑149 (MQ=255)
                                                                                                                 tGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAg                                <  1:905918/149‑1 (MQ=255)
                                                                                                                  gTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAgg                               <  1:1501190/149‑1 (MQ=255)
                                                                                                                  gTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAgg                               <  1:1286465/149‑1 (MQ=255)
                                                                                                                      ttCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAGGCCat                           <  2:328956/149‑1 (MQ=255)
                                                                                                                        cGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAGGCCatat                         <  2:2150209/149‑1 (MQ=255)
                                                                                                                          ggTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAGGCCATATcc                       <  2:2360304/149‑1 (MQ=255)
                                                                                                                          ggTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAGGCCATATcc                       >  2:3442826/1‑149 (MQ=255)
                                                                                                                           gTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAGGCCATATccc                      <  2:2730393/149‑1 (MQ=255)
                                                                                                                           gTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAGGCCATATccc                      <  2:548569/149‑1 (MQ=255)
                                                                                                                            ttATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAGGCCATATCCCa                     <  1:1400210/149‑1 (MQ=255)
                                                                                                                              ataGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAGGCCATATCCCAGa                   >  1:1976728/1‑149 (MQ=255)
                                                                                                                              ataGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAGGCCATATCCCAGa                   >  2:1912093/1‑149 (MQ=255)
                                                                                                                                 gCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAGGCCATATCCCAGACGa                <  2:1513588/149‑1 (MQ=255)
                                                                                                                                 gCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAGGCCATATCCCAGACGa                >  1:2748285/1‑149 (MQ=255)
                                                                                                                                 gCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAGGCCATATCCCAGACGa                >  2:756743/1‑149 (MQ=255)
                                                                                                                                    aCGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAGGCCATATCCCAGACGATGc             >  2:3573284/1‑149 (MQ=255)
                                                                                                                                      gCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAGGCCATATCCCAGACGATGCCg           <  1:3641701/149‑1 (MQ=255)
                                                                                                                                        aGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAGGCCATATCCCAGACGATGCCGtt         <  2:3649899/149‑1 (MQ=255)
                                                                                                                                         gACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAGGCCATATCCCAGACGATGCCGTTg        >  2:1830332/1‑149 (MQ=255)
                                                                                                                                         gACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAGGCCATATCCCAGACGATGCCGTTg        >  2:429022/1‑149 (MQ=255)
                                                                                                                                         gACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAGGCCATATCCCAGACGATGCCGTTg        >  2:1081098/1‑149 (MQ=255)
                                                                                                                                         gACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAGGCCATATCCCAGACGATGCCGTTg        >  1:3688368/1‑149 (MQ=255)
                                                                                                                                          aCCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAGGCCATATCCCAGACGATGCCGTTGc       >  2:247532/1‑149 (MQ=255)
                                                                                                                                              cTGCGAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAGGCCATATCCCAGACGATGCCGTTGCTGtt   <  2:454111/149‑1 (MQ=255)
                                                                                                                                               tGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAGGCCATATCCCAGACGATGCCGTTGCTGTTa  <  1:1488348/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CATCGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCGAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAGGCCATATCCCAGACGATGCCGTTGCTGTTA  >  CP009273/1860037‑1860328

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: