Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,459,603 A→T 43.0% intergenic (+24/+5) gspO → / ← bfr bifunctional prepilin leader peptidase/ methylase/bacterioferritin, iron storage and detoxification protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,459,6030AT43.0% 103.8 / 300.1 214intergenic (+24/+5)gspO/bfrbifunctional prepilin leader peptidase/ methylase/bacterioferritin, iron storage and detoxification protein
Reads supporting (aligned to +/- strand):  ref base A (66/56);  new base T (46/46);  total (112/102)
Fisher's exact test for biased strand distribution p-value = 5.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-01

Only 100 of 218 total aligned reads displayed.

TTGGTTACGCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATG  >  CP009273/3459458‑3459751
                                                                                                                                                 |                                                                                                                                                    
ttGGTTACGCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACgg                                                                                                                                                   >  2:2585365/1‑149 (MQ=255)
ttGGTTACGCCATCGTTTCACAACTTATGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACgg                                                                                                                                                   >  2:2141674/1‑149 (MQ=255)
   gTTACGCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTc                                                                                                                                                <  2:3579800/149‑1 (MQ=255)
   gTTACGCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTc                                                                                                                                                <  2:1315547/149‑1 (MQ=255)
   gTTACGCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTc                                                                                                                                                <  2:3365355/149‑1 (MQ=255)
          cATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACcttc                                                                                                                                         <  1:631204/149‑1 (MQ=255)
           aTCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATGTGGCCTGGATCTCTTATTACCCGGGTGGGTTCTCCCTCCCGACACGGCTCAACcttct                                                                                                                                        <  1:698823/149‑1 (MQ=255)
           aTCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACcttct                                                                                                                                        <  2:1078734/149‑1 (MQ=255)
            tCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACcttctt                                                                                                                                       <  2:3185458/149‑1 (MQ=255)
                 tCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCgg                                                                                                                                  >  2:1880487/1‑149 (MQ=255)
                  cacaACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGa                                                                                                                                 >  1:4021338/1‑149 (MQ=255)
                    caACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATc                                                                                                                               <  1:1605432/149‑1 (MQ=255)
                    caACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATc                                                                                                                               <  1:3383529/149‑1 (MQ=255)
                    caACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATc                                                                                                                               <  1:2038639/149‑1 (MQ=255)
                         tCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGc                                                                                                                          <  1:624246/149‑1 (MQ=255)
                          cTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCt                                                                                                                         <  2:297222/149‑1 (MQ=255)
                          cTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCt                                                                                                                         <  2:598532/149‑1 (MQ=255)
                           tGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCtt                                                                                                                        >  2:2203642/1‑149 (MQ=255)
                           tGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTGATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCtt                                                                                                                        >  1:1505949/1‑149 (MQ=255)
                            gCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTg                                                                                                                       >  1:1494346/1‑149 (MQ=255)
                                cGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGTTCTCTTATTAACCGGTAGTGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGa                                                                                                                   <  1:1553139/149‑1 (MQ=255)
                                cGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGa                                                                                                                   <  1:1532867/149‑1 (MQ=255)
                                 gCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGAt                                                                                                                  >  1:178230/1‑149 (MQ=255)
                                   cGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATaa                                                                                                                <  2:2724449/149‑1 (MQ=255)
                                   cGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATaa                                                                                                                <  1:2172394/149‑1 (MQ=255)
                                    ggACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAAt                                                                                                               <  1:1423028/149‑1 (MQ=255)
                                    ggACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAAt                                                                                                               >  2:1997703/1‑149 (MQ=255)
                                    ggACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAAt                                                                                                               <  1:3935465/149‑1 (MQ=255)
                                     gACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAAtt                                                                                                              >  2:1779964/1‑149 (MQ=255)
                                       cTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAAtttt                                                                                                            <  1:2095608/149‑1 (MQ=255)
                                       cTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAAtttt                                                                                                            <  1:3177632/149‑1 (MQ=255)
                                           aaCCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAg                                                                                                        <  1:2982500/149‑1 (MQ=255)
                                            aCCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAga                                                                                                       <  1:503843/149‑2 (MQ=255)
                                             ccACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGc                                                                                                      >  2:3065560/1‑149 (MQ=255)
                                               aCTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGccc                                                                                                    >  1:1581218/1‑149 (MQ=255)
                                                cTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGACCa                                                                                                   >  1:2197208/1‑149 (MQ=255)
                                                     gCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCAtcttc                                                                                              <  1:3640096/149‑1 (MQ=255)
                                                       aTTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTg                                                                                            >  2:2931400/1‑149 (MQ=255)
                                                        ttCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGa                                                                                           >  1:1982218/1‑149 (MQ=255)
                                                           ggCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATc                                                                                        <  1:481134/149‑1 (MQ=255)
                                                            gCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCa                                                                                       >  1:1787601/1‑149 (MQ=255)
                                                            gCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCa                                                                                       <  1:167946/149‑1 (MQ=255)
                                                               cGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCCGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGAt                                                                                    <  1:4106623/149‑1 (MQ=255)
                                                                gTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATc                                                                                   >  2:1195230/1‑149 (MQ=255)
                                                                 tGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCa                                                                                  <  2:3298942/149‑1 (MQ=255)
                                                                   gCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAg                                                                                <  2:2938409/149‑1 (MQ=255)
                                                                   gCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAg                                                                                >  2:3891358/1‑149 (MQ=255)
                                                                    cTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGt                                                                               <  1:983735/149‑1 (MQ=255)
                                                                    cTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGt                                                                               <  2:178260/149‑1 (MQ=255)
                                                                    cTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGt                                                                               <  2:3314411/149‑1 (MQ=255)
                                                                     ttGCGCTCGGCAGCATGATAAACCTGGGGGATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGtt                                                                              <  2:345460/149‑1 (MQ=255)
                                                                         gcTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGt                                                                          <  2:2500793/149‑1 (MQ=255)
                                                                          cTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGtt                                                                         >  2:281418/1‑149 (MQ=255)
                                                                          cTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGtt                                                                         >  2:1958653/1‑149 (MQ=255)
                                                                           tCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGttt                                                                        >  1:2488859/1‑149 (MQ=255)
                                                                           tCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGttt                                                                        >  2:4064137/1‑149 (MQ=255)
                                                                            cGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTc                                                                       >  1:4076869/1‑149 (MQ=255)
                                                                             ggCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTcc                                                                      <  2:1505949/149‑1 (MQ=255)
                                                                                agcaTGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTccagc                                                                   <  2:175967/149‑1 (MQ=255)
                                                                                 gcaTGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTccagcc                                                                  >  1:2215848/1‑149 (MQ=255)
                                                                                 gcaTGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTccagcc                                                                  >  1:3472410/1‑149 (MQ=255)
                                                                                  caTGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTccagcca                                                                 <  2:2190896/149‑1 (MQ=255)
                                                                                    tGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGt                                                               >  1:1334548/1‑149 (MQ=255)
                                                                                     gATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTc                                                              <  1:1805328/149‑1 (MQ=255)
                                                                                     gATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTc                                                              <  1:1093953/149‑1 (MQ=255)
                                                                                     gATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTc                                                              <  2:495955/149‑1 (MQ=255)
                                                                                      aTAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCg                                                             <  1:1439624/149‑1 (MQ=255)
                                                                                       tacaCCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGa                                                            <  2:2639661/146‑1 (MQ=255)
                                                                                       tAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGa                                                            >  2:3861027/1‑149 (MQ=255)
                                                                                       tAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGa                                                            <  1:2788355/149‑1 (MQ=255)
                                                                                          aCCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGatat                                                         <  1:3897844/149‑1 (MQ=255)
                                                                                            cTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATgg                                                       >  2:3247026/1‑149 (MQ=255)
                                                                                             tGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGc                                                      >  2:3889093/1‑149 (MQ=255)
                                                                                              ggggTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGcc                                                     >  1:2481002/1‑149 (MQ=255)
                                                                                                ggTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCctt                                                   <  1:3119623/149‑1 (MQ=255)
                                                                                                  tATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCcttct                                                 <  2:2566702/149‑1 (MQ=255)
                                                                                                    tCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCttcttc                                               >  2:169596/1‑149 (MQ=255)
                                                                                                      tGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTtcat                                             >  2:411870/1‑149 (MQ=255)
                                                                                                       ggCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTtcatc                                            >  2:2057231/1‑149 (MQ=255)
                                                                                                       ggCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTtcatc                                            >  2:1143650/1‑149 (MQ=255)
                                                                                                         ccTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCAc                                          <  2:1581218/149‑1 (MQ=255)
                                                                                                         ccTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCAc                                          >  1:572266/1‑149 (MQ=255)
                                                                                                            ggATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCa                                       >  2:1660950/1‑149 (MQ=255)
                                                                                                                ctctTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAAt                                   <  1:3360219/149‑1 (MQ=255)
                                                                                                                   ttattaACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTc                                <  2:2447839/149‑1 (MQ=255)
                                                                                                                    tattaACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCt                               >  1:2725502/1‑149 (MQ=255)
                                                                                                                      ttaACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTat                             >  2:2057886/1‑149 (MQ=255)
                                                                                                                        aaCCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTatca                           >  2:519897/1‑149 (MQ=255)
                                                                                                                        aaCCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTatca                           >  1:1055892/1‑149 (MQ=255)
                                                                                                                         aCCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTatcat                          >  1:930178/1‑149 (MQ=255)
                                                                                                                               aGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATc                    >  2:3884157/1‑149 (MQ=255)
                                                                                                                                  gTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATcgcg                 >  1:1216738/1‑149 (MQ=255)
                                                                                                                                      tccctcccGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTg             >  2:1375710/1‑149 (MQ=255)
                                                                                                                                       ccctcccGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGa            >  1:2383034/1‑149 (MQ=255)
                                                                                                                                          tcccGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTATGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCATTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGt         >  1:493287/1‑149 (MQ=255)
                                                                                                                                            ccGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTaa       <  2:3927368/149‑1 (MQ=255)
                                                                                                                                             cGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAAt      <  1:4060174/149‑1 (MQ=255)
                                                                                                                                             cGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAAt      <  1:1362362/149‑1 (MQ=255)
                                                                                                                                               acacGGCTCAACCTTCGTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCa    <  2:2271670/149‑1 (MQ=255)
                                                                                                                                                 acGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATg  >  2:3536836/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TTGGTTACGCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATG  >  CP009273/3459458‑3459751

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: