Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,459,603 A→T 63.1% 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,6030AT63.1% 217.7 / 241.8 214intergenic (+24/+5)gspO/bfrbifunctional prepilin leader peptidase/ methylase/bacterioferritin, iron storage and detoxification protein
Reads supporting (aligned to +/- strand):  ref base A (42/37);  new base T (66/69);  total (108/106)
Fisher's exact test for biased strand distribution p-value = 5.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 215 total aligned reads displayed.

TTGGTTACGCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATC  >  CP009273/3459458‑3459748
                                                                                                                                                 |                                                                                                                                                 
ttGGTTACGCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCg                                                                                                                                                 >  1:918673/1‑148 (MQ=255)
 tGGTTACGCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCgg                                                                                                                                                <  2:2256213/148‑1 (MQ=255)
     tACGCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCa                                                                                                                                            >  2:2872395/1‑148 (MQ=255)
      aCGCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCaa                                                                                                                                           >  1:2874287/1‑148 (MQ=255)
      aCGCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCaa                                                                                                                                           >  1:1528504/1‑148 (MQ=255)
       cGCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAAc                                                                                                                                          <  1:2793278/148‑1 (MQ=255)
        gCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAAcc                                                                                                                                         <  2:1612973/148‑1 (MQ=255)
        gCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAAcc                                                                                                                                         <  2:2112515/148‑1 (MQ=255)
          cATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACctt                                                                                                                                       <  1:1675093/148‑1 (MQ=255)
           aTCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACcttc                                                                                                                                      >  2:1803820/1‑148 (MQ=255)
             cGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACcttctt                                                                                                                                    <  2:3368058/148‑1 (MQ=255)
             cGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACcttctt                                                                                                                                    >  2:2270798/1‑148 (MQ=255)
                ttCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTcgc                                                                                                                                 >  1:894891/1‑148 (MQ=255)
                   acaACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGa                                                                                                                              <  1:2866644/148‑1 (MQ=255)
                    caACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGAt                                                                                                                             <  1:544665/148‑1 (MQ=255)
                    caACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGAt                                                                                                                             <  2:350474/148‑1 (MQ=255)
                    caACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGAt                                                                                                                             <  2:2663377/148‑1 (MQ=255)
                     aaCTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGTATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATc                                                                                                                            <  1:3350523/148‑1 (MQ=255)
                      aCTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCt                                                                                                                           >  2:1160347/1‑148 (MQ=255)
                       cTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCtg                                                                                                                          >  1:1879588/1‑148 (MQ=255)
                       cTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCtg                                                                                                                          <  2:805475/148‑1 (MQ=255)
                          cTGAAACGCCTGACTATAACCAATATCGATTTCGGCCGGTGGGTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGTTCTCTTTTTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGc                                                                                                                       <  1:1696232/148‑1 (MQ=255)
                             cAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTg                                                                                                                    <  1:1484446/148‑1 (MQ=255)
                             cAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTg                                                                                                                    <  1:733110/148‑1 (MQ=255)
                             cAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTg                                                                                                                    <  1:1484438/148‑1 (MQ=255)
                             cAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTg                                                                                                                    >  1:119337/1‑148 (MQ=255)
                             cAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTg                                                                                                                    >  2:691815/1‑148 (MQ=255)
                             cAACGCCGGACTATAACCACTATCGCATTCGGCCCGAGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTg                                                                                                                    >  1:3317277/1‑148 (MQ=255)
                              aaCGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGc                                                                                                                   <  2:2354728/148‑1 (MQ=255)
                                cGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAg                                                                                                                 >  2:2939521/1‑148 (MQ=255)
                                cGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAg                                                                                                                 >  1:1045325/1‑148 (MQ=255)
                                cGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAg                                                                                                                 >  1:2612833/1‑148 (MQ=255)
                                cGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAg                                                                                                                 <  2:567693/148‑1 (MQ=255)
                                    ggACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATaa                                                                                                             <  1:3252865/148‑1 (MQ=255)
                                      aCTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAAtt                                                                                                           <  1:1414477/148‑1 (MQ=255)
                                           aaCCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCa                                                                                                      <  1:1028305/148‑1 (MQ=255)
                                            aCCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAg                                                                                                     >  1:1992804/1‑148 (MQ=255)
                                            aCCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAg                                                                                                     >  2:216120/1‑148 (MQ=255)
                                               aCTCTCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGcc                                                                                                  <  2:3016509/148‑1 (MQ=255)
                                                cTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGccc                                                                                                 >  2:145065/1‑148 (MQ=255)
                                                 tATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCa                                                                                                <  2:173220/148‑1 (MQ=255)
                                                  aTCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCAt                                                                                               >  1:301404/1‑148 (MQ=255)
                                                   tCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCAtc                                                                                              <  1:463840/148‑1 (MQ=255)
                                                    cGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCAtct                                                                                             <  2:1705931/148‑1 (MQ=255)
                                                    cGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCAtct                                                                                             >  1:417945/1‑148 (MQ=255)
                                                    cGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCAtct                                                                                             >  1:1953094/1‑148 (MQ=255)
                                                     gCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCAtctt                                                                                            >  1:2982199/1‑148 (MQ=255)
                                                      cATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCAtcttc                                                                                           >  2:2715182/1‑148 (MQ=255)
                                                       aTTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCAtcttct                                                                                          >  1:641463/1‑148 (MQ=255)
                                                       aTTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACAGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCAGATCTGTGCTTGCAGATAATTTTGCAGGCCCAtcttct                                                                                          <  2:311392/148‑1 (MQ=255)
                                                        ttCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCAATCTTCTg                                                                                         >  1:3452700/1‑148 (MQ=255)
                                                         tCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGa                                                                                        >  1:1925032/1‑148 (MQ=255)
                                                          cGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGaa                                                                                       <  2:2172634/148‑1 (MQ=255)
                                                          cGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGaa                                                                                       >  2:360931/1‑148 (MQ=255)
                                                           ggCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAAt                                                                                      <  1:2963932/148‑1 (MQ=255)
                                                                  ggCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCa                                                                               >  2:468401/1‑148 (MQ=255)
                                                                   gCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCaa                                                                              <  1:2609979/148‑1 (MQ=255)
                                                                    cTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAg                                                                             <  1:694392/148‑1 (MQ=255)
                                                                    cTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAg                                                                             <  1:307773/148‑1 (MQ=255)
                                                                           tCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGtt                                                                      <  2:1511962/148‑1 (MQ=255)
                                                                           tCGGCACCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGtt                                                                      <  2:3107584/148‑1 (MQ=255)
                                                                               cagcaTGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTACGTTTcca                                                                  >  1:1839752/1‑148 (MQ=255)
                                                                                 gcaTGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTccagc                                                                >  2:914795/1‑148 (MQ=255)
                                                                                 gcaTGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTccagc                                                                >  2:1416082/1‑148 (MQ=255)
                                                                                     gATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGt                                                            <  1:1112733/148‑1 (MQ=255)
                                                                                      aTAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTc                                                           >  1:2441422/1‑148 (MQ=255)
                                                                                       tAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCg                                                          <  2:1766048/148‑1 (MQ=255)
                                                                                        aaaCCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCAGTTTCCAGCCAGTCGa                                                         >  2:3341660/1‑148 (MQ=255)
                                                                                         aaCCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGat                                                        <  2:1229052/148‑1 (MQ=255)
                                                                                         aaCCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGat                                                        <  1:3390666/148‑1 (MQ=255)
                                                                                           ccTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGatat                                                      <  1:402880/148‑1 (MQ=255)
                                                                                           ccTGGGATATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGatat                                                      <  1:1210139/148‑1 (MQ=255)
                                                                                             tGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATgg                                                    >  2:2759402/1‑148 (MQ=255)
                                                                                             tGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATgg                                                    <  1:2270798/148‑1 (MQ=255)
                                                                                              ggggTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGc                                                   <  2:818540/148‑1 (MQ=255)
                                                                                               gggTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGcc                                                  >  2:1514572/1‑148 (MQ=255)
                                                                                               gggTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGcc                                                  >  2:246682/1‑148 (MQ=255)
                                                                                                ggTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCct                                                 <  2:2818533/148‑1 (MQ=255)
                                                                                                  tATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCcttc                                               <  1:1160347/148‑1 (MQ=255)
                                                                                                    tCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCcttctt                                             >  1:133791/1‑148 (MQ=255)
                                                                                                         ccTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTtcatca                                        >  2:154453/1‑148 (MQ=255)
                                                                                                         ccTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTtcatca                                        >  2:3240467/1‑148 (MQ=255)
                                                                                                          cTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCAc                                       >  2:3404550/1‑148 (MQ=255)
                                                                                                              aTCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCaa                                   >  1:1777652/1‑148 (MQ=255)
                                                                                                               tctcTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCaaa                                  >  1:3103613/1‑148 (MQ=255)
                                                                                                                 tctTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAAt                                >  1:1983708/1‑148 (MQ=255)
                                                                                                                    tattaACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTc                             <  1:575232/148‑1 (MQ=255)
                                                                                                                      ttaACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTa                           <  2:417900/148‑1 (MQ=255)
                                                                                                                        aaCCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTtctctc                         >  2:2441508/1‑148 (MQ=255)
                                                                                                                        aaCCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTatc                         <  2:1103968/148‑1 (MQ=255)
                                                                                                                        aaCCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTGCTGAATCAGATCAAGTTCCGTTTCCAGCCAGGCGATATGGCCTGCTTCATCACGCCACCTTTCTatc                         >  2:1898874/1‑148 (MQ=255)
                                                                                                                         aCCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTatca                        <  2:3180285/148‑1 (MQ=255)
                                                                                                                          ccGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTatcat                       <  1:2369313/148‑1 (MQ=255)
                                                                                                                          ccGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGAAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTatcat                       >  1:655779/1‑148 (MQ=255)
                                                                                                                             ggAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATcatcat                    <  1:3422425/148‑1 (MQ=255)
                                                                                                                                          tcccGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACg       >  2:426545/1‑148 (MQ=255)
                                                                                                                                           cccGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGt      <  2:1968353/148‑1 (MQ=255)
                                                                                                                                            ccGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTa     <  2:850816/148‑1 (MQ=255)
                                                                                                                                             cGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTaa    >  1:1442089/1‑148 (MQ=255)
                                                                                                                                              gACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAAt   >  2:2684028/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TTGGTTACGCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATC  >  CP009273/3459458‑3459748

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: