Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,559,639 A→C 5.1% Y118D (TAC→GAC)  dosP ← oxygen sensor, c‑di‑GMP phosphodiesterase, heme‑regulated; cold‑ and stationary phase‑induced bioflim regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,559,6390AC 5.1% 513.1 / 18.4 176Y118D (TAC→GAC) dosPoxygen sensor, c‑di‑GMP phosphodiesterase, heme‑regulated; cold‑ and stationary phase‑induced bioflim regulator
Reads supporting (aligned to +/- strand):  ref base A (99/68);  new base C (3/6);  total (102/74)
Fisher's exact test for biased strand distribution p-value = 1.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

Only 100 of 183 total aligned reads displayed.

GTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAGCTCCCGACTCATCCCCTCAACACGCGCTTTACCGCCTTCACGGTTGTGACGAATGTATTCAGGATG  >  CP009273/1559494‑1559786
                                                                                                                                                 |                                                                                                                                                   
gTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGtaat                                                                                                                                                   <  2:1625729/148‑1 (MQ=255)
  aaaTGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGtaataa                                                                                                                                                 <  2:1548239/148‑1 (MQ=255)
    aTGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAAc                                                                                                                                               >  1:2617248/1‑148 (MQ=255)
      gcgcGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACtt                                                                                                                                             <  1:100621/148‑1 (MQ=255)
      gcgcGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACtt                                                                                                                                             <  1:1616913/148‑1 (MQ=255)
            tagcactgcacAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCt                                                                                                                                       <  2:213914/146‑1 (MQ=255)
                actgcacAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGgc                                                                                                                                   >  1:852570/1‑148 (MQ=255)
                    cacaATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTc                                                                                                                               >  2:2851431/1‑148 (MQ=255)
                     acaATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCa                                                                                                                              <  1:2141385/148‑1 (MQ=255)
                       aaTATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACt                                                                                                                            <  2:3216202/148‑1 (MQ=255)
                        atatGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACtt                                                                                                                           >  2:2800078/1‑148 (MQ=255)
                           tGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACATTCg                                                                                                                        >  1:2876068/1‑148 (MQ=255)
                            ggCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGa                                                                                                                       <  1:1736790/148‑1 (MQ=255)
                            ggCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGa                                                                                                                       >  2:1846405/1‑148 (MQ=255)
                            ggCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGa                                                                                                                       >  1:470700/1‑148 (MQ=255)
                             gCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGAt                                                                                                                      <  1:889696/148‑1 (MQ=255)
                              cGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATa                                                                                                                     >  2:2689344/1‑148 (MQ=255)
                              cGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATa                                                                                                                     >  2:1056898/1‑148 (MQ=255)
                                ttCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAgc                                                                                                                   >  2:752186/1‑148 (MQ=255)
                                 tCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAgcg                                                                                                                  >  1:1786744/1‑148 (MQ=255)
                                  ccGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAgcgc                                                                                                                 >  1:1023573/1‑148 (MQ=255)
                                  ccGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAgcgc                                                                                                                 >  2:2138531/1‑148 (MQ=255)
                                   cGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCa                                                                                                                <  1:97608/148‑1 (MQ=255)
                                     gATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCaaa                                                                                                              <  2:1165864/148‑1 (MQ=255)
                                     gATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCaaa                                                                                                              <  2:184236/148‑1 (MQ=255)
                                     gATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCaaa                                                                                                              <  2:3123796/148‑1 (MQ=255)
                                      aTCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCaaaa                                                                                                             >  2:1756972/1‑148 (MQ=255)
                                      aTCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCaaaa                                                                                                             <  2:603802/148‑1 (MQ=255)
                                      aTCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCaaaa                                                                                                             <  2:1485718/148‑1 (MQ=255)
                                      aTCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCaaaa                                                                                                             >  1:1519625/1‑148 (MQ=255)
                                       tCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAAc                                                                                                            >  1:1369622/1‑148 (MQ=255)
                                       tCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAAc                                                                                                            >  1:186129/1‑148 (MQ=255)
                                          aGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACggg                                                                                                         >  2:3516987/1‑148 (MQ=255)
                                            gACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACGTTCGATAGCGCAAAACGGGTc                                                                                                       >  1:3298704/1‑148 (MQ=255)
                                             aCAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTcc                                                                                                      >  2:3484964/1‑148 (MQ=255)
                                               aaTCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAg                                                                                                    >  2:2173161/1‑148 (MQ=255)
                                                aTCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGa                                                                                                   <  2:2058044/148‑1 (MQ=255)
                                                 tCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCCGAt                                                                                                  >  2:3371615/1‑148 (MQ=255)
                                                 tCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGAt                                                                                                  <  2:1016688/148‑1 (MQ=255)
                                                   aCCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGAttt                                                                                                <  2:1283536/148‑1 (MQ=255)
                                                     cggtcggtCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTa                                                                                              <  1:679677/148‑1 (MQ=255)
                                                      ggtcggtcGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTAc                                                                                             <  2:2005542/148‑1 (MQ=255)
                                                        tcggtcgAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTg                                                                                           >  1:3099606/1‑148 (MQ=255)
                                                           gtcgAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTCATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCg                                                                                        >  1:996150/1‑148 (MQ=255)
                                                            tcgAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGt                                                                                       >  1:2544487/1‑148 (MQ=255)
                                                            tcgAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGt                                                                                       >  1:761813/1‑148 (MQ=255)
                                                            tcgAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGt                                                                                       >  1:2485202/1‑148 (MQ=255)
                                                            tcgAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGt                                                                                       >  1:1191350/1‑148 (MQ=255)
                                                              gagaTGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCt                                                                                     <  2:3067370/148‑1 (MQ=255)
                                                                gaTGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCttt                                                                                   <  1:1058148/148‑1 (MQ=255)
                                                                    gTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTtc                                                                               >  2:641069/1‑148 (MQ=255)
                                                                    gTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTtc                                                                               >  2:638770/1‑148 (MQ=255)
                                                                      cAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTtctc                                                                             <  2:1893801/148‑1 (MQ=255)
                                                                      cAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGAATTTACTGCCGTCTTTTTtctc                                                                             >  2:3503331/1‑148 (MQ=255)
                                                                      aaaCGGCAATAATCAATTGTCTGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAATCTATACCCTCGGCGCTCACTATCGATAGAGCAAAAACGGTCCAGAATTTACTGCCGTCTTTTTtctc                                                                             >  2:3307390/2‑148 (MQ=255)
                                                                         cGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAg                                                                          >  2:1298915/1‑148 (MQ=255)
                                                                            cAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTg                                                                       <  1:2857500/148‑1 (MQ=255)
                                                                               taatCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAg                                                                    >  2:473571/1‑148 (MQ=255)
                                                                                  tCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAGCTc                                                                 <  2:1790581/148‑1 (MQ=255)
                                                                                           gggTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAGCTCCCGACTCAt                                                        >  1:2611252/1‑148 (MQ=255)
                                                                                           gggTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAGCTCCCGACTCAt                                                        <  2:226850/148‑1 (MQ=255)
                                                                                            ggTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAGCTCCCGACTCATc                                                       >  1:2308282/1‑148 (MQ=255)
                                                                                             gTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTATCCAGCTGCAGCACCCGACACATcc                                                      >  1:1670522/1‑148 (MQ=255)
                                                                                                   tCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAGCTCCCGACTCATCCCCACaa                                                >  1:2738721/1‑148 (MQ=255)
                                                                                                    cTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAGCTCCCGACTCATCCCCTCAac                                               >  2:482636/1‑148 (MQ=255)
                                                                                                      ttttGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAGCTCCCGACTCATCCCCTCAacac                                             >  2:21184/1‑148 (MQ=255)
                                                                                                         tGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAGCTCCCGACTCATCCCCTCAACAcgcg                                          <  1:1055293/148‑1 (MQ=255)
                                                                                                           cgcCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAGCTCCCGACTCATCCCCACAACACGCGCt                                        >  1:2548609/1‑148 (MQ=255)
                                                                                                             ccATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAGCTCCCGACTCATCCCCTCAACACGCGCttt                                      >  2:2454305/1‑148 (MQ=255)
                                                                                                             ccATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAGCTCCCGACTCATCCCCTCAACACGCGCttt                                      >  2:2454362/1‑148 (MQ=255)
                                                                                                              cATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAGCTCCCGACTCATCCCCTCAACACGCGCTTTa                                     >  2:135383/1‑148 (MQ=255)
                                                                                                               aTTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAGCTCCCGACTCATCCCCTCAACACGCGCTTTAc                                    >  2:2087904/1‑148 (MQ=255)
                                                                                                               aTTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAGCTCCCGACTCATCCCCTCAACACGCGCTTTAc                                    >  2:292345/1‑148 (MQ=255)
                                                                                                                 ttCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTATCCAGCTGCAGCTCCCGACTCATCCCCTCACCACGCGCGTTACCg                                  >  1:3011362/1‑148 (MQ=255)
                                                                                                                  tCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAGCTCCCGACTCATCCCCTCAACACGCGCTTTACCGc                                 <  2:3416152/148‑1 (MQ=255)
                                                                                                                  tCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAGCTCCCGACTCATCCCCTCAACACGCGCTTTACCGc                                 <  2:835611/148‑1 (MQ=255)
                                                                                                                   cTACGCTTGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACAGCCGTCTTTTTTCTCCAGCTGCAGCTCCCGACTCATCCCCTAAAAACGCGCTTTACCGcc                                >  2:547944/1‑148 (MQ=255)
                                                                                                                   cTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAGCTCCCGACTCATCCCCTCAACACGCGCTTTACCGcc                                >  2:1018972/1‑148 (MQ=255)
                                                                                                                     aCGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAGCTCCCGACTCATCCCCTCAACACGCGCTTTACCGCCtt                              <  2:3013070/148‑1 (MQ=255)
                                                                                                                     aCGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAGCTCCCGACTCATCCCCTCAACACGCGCTTTACCGCCtt                              >  2:2523397/1‑148 (MQ=255)
                                                                                                                      cggtGGCATCCGGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGAAAAAGGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAGCTCCCGACTCATCCCCTCAACACGCGCTTTACCGCCTTc                             <  2:1307909/145‑1 (MQ=255)
                                                                                                                      cGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCGTCGGCGGTCACTTTCGATAGCGCAAAACGGGTCCTGATTTTACTGCCGTCTTTTTTCTCCAGCTGCTGCTCCCGACTCATCCCCTCAACACGCGCTTTACCGCCTTc                             <  2:3077254/148‑1 (MQ=255)
                                                                                                                          ggCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAGCTCCCGACTCATCCCCTCAACACGCGCTTTACCGCCTTCACgg                         <  2:968166/148‑1 (MQ=255)
                                                                                                                          ggCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAGCTCCCGACTCATCCCCTCAACACGCGCTTTACCGCCTTCACgg                         <  1:599914/148‑1 (MQ=255)
                                                                                                                          ggCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAGCTCCCGACTCATCCCCTCAACACGCGCTTTACCGCCTTCACgg                         <  1:2858302/148‑1 (MQ=255)
                                                                                                                           gCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAGCTCCCGACTCATCCCCTCAACACGCGCTTTACCGCCTTCACGGt                        <  1:650710/148‑1 (MQ=255)
                                                                                                                            cATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAGCTCCCGACTCATCCCCTCAACACGCGCTTTACCGCCTTCACGGtt                       <  2:1852969/148‑1 (MQ=255)
                                                                                                                            cATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAGCTCCCGACTCATCCCCTCAACACGCGCTTTACCGCCTTCACGGtt                       <  2:932980/148‑1 (MQ=255)
                                                                                                                            cATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAGCTCCCGACTCATCCCCTCAACACGCGCTTTACCGCCTTCACGGtt                       <  1:2173161/148‑1 (MQ=255)
                                                                                                                                 cGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAGCTCCCGACTCATCCCCTCAACACGCGCTTTACCGCCTTCACGGTTGTGAc                  >  2:2213708/1‑148 (MQ=255)
                                                                                                                                   tACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAGCTCCCGACTCATCCCCTCAACACGCGCTTTACCGCCTTCACGGTTGTGACGa                >  2:3243982/1‑148 (MQ=255)
                                                                                                                                     ccAGCGCCAGGTCATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAGCTCCCGACTCATCCCCTCAACACGCGCTTTACCGCCTTCACGGTTGTGACGAAt              <  1:598126/148‑1 (MQ=255)
                                                                                                                                     ccAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAGCTCCCGACTCATCCCCTCAACACGCGCTTTACCGCCTTCACGGTTGTGACGAAt              <  2:3395625/148‑1 (MQ=255)
                                                                                                                                      cAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAGCTCCCGACTCATCCCCTCAACACGCGCTTTACCGCCTTCACGGTTGTGACGAATg             <  1:2366695/148‑1 (MQ=255)
                                                                                                                                           ccAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAGCTCCCGACTCATCCCCTCAACACGCGCTTTACCGCCTTCACGGTTGTGACGAATGTATTc        >  2:2400763/1‑148 (MQ=255)
                                                                                                                                             aGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAGCTCCCGACTCATCCCCTCAACACGCGCTTTACCGCCGTCACGGTTGTGACGAATGTATTCAg      >  2:1389408/1‑148 (MQ=255)
                                                                                                                                              ggTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACCGGTCCAGATTTTACTGGCGTCTTTTTTATCCAGCAGAAGCTCCCGACTCAAACCCACAACACCCGCCTTAACCCCATCACGGTTTTGACGAATTTATTAAgg     >  1:2757791/1‑148 (MQ=255)
                                                                                                                                               gTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAGCTCCCGACTCATCCCCTCAACACGCGCTTTACCGCCTTCACGGTTGTGACGAATGTATTCAGGa    >  1:3200747/1‑148 (MQ=255)
                                                                                                                                                 aataaACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAGCTCCCGACTCATCCCCTCAACACGCGCTTTACCGCCTTCACGGTTGTGACGAATGTATTCAGGATg  <  1:3226742/148‑1 (MQ=255)
                                                                                                                                                 aataaACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAGCTCCCGACTCATCCCCTCAACACGCGCTTTACCGCCTTCACGGTTGTGACGAATGTATTCAGGATg  <  1:1459020/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCGGCGCTCACTTTCGATAGCGCAAAACGGGTCCAGATTTTACTGCCGTCTTTTTTCTCCAGCTGCAGCTCCCGACTCATCCCCTCAACACGCGCTTTACCGCCTTCACGGTTGTGACGAATGTATTCAGGATG  >  CP009273/1559494‑1559786

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