Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,459,603 A→T 40.3% 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,6030AT40.3% 147.8 / 274.6 211intergenic (+24/+5)gspO/bfrbifunctional prepilin leader peptidase/ methylase/bacterioferritin, iron storage and detoxification protein
Reads supporting (aligned to +/- strand):  ref base A (58/68);  new base T (46/39);  total (104/107)
Fisher's exact test for biased strand distribution p-value = 2.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.29e-01

Only 100 of 213 total aligned reads displayed.

TTGGTTACGCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCA  >  CP009273/3459458‑3459749
                                                                                                                                                 |                                                                                                                                                  
 tGGTTACGCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCgg                                                                                                                                                 >  2:3613538/1‑148 (MQ=255)
 tGGTTACGCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCgg                                                                                                                                                 >  2:2493611/1‑148 (MQ=255)
   gTTACGCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCt                                                                                                                                               <  1:3321613/148‑1 (MQ=255)
    ttACGCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTc                                                                                                                                              >  2:3664743/1‑148 (MQ=255)
     tACGCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCa                                                                                                                                             >  1:3533652/1‑148 (MQ=255)
            tCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACcttct                                                                                                                                      <  1:594227/148‑1 (MQ=255)
            tCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACcttct                                                                                                                                      >  2:525630/1‑148 (MQ=255)
            tCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACcttct                                                                                                                                      >  2:2686226/1‑148 (MQ=255)
             cGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACcttctt                                                                                                                                     >  2:801112/1‑148 (MQ=255)
                ttCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTcgc                                                                                                                                  >  1:301805/1‑148 (MQ=255)
                ttCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTcgc                                                                                                                                  >  2:65736/1‑148 (MQ=255)
                   acaACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGa                                                                                                                               <  2:1463943/148‑1 (MQ=255)
                    caACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGAt                                                                                                                              <  1:512090/148‑1 (MQ=255)
                    caACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGAt                                                                                                                              >  1:1130510/1‑148 (MQ=255)
                    caACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGAt                                                                                                                              <  1:1251110/148‑1 (MQ=255)
                        ttCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCtgt                                                                                                                          >  2:1068368/1‑148 (MQ=255)
                         tCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCtgtg                                                                                                                         <  1:2301390/148‑1 (MQ=255)
                                cGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAg                                                                                                                  >  2:1692180/1‑148 (MQ=255)
                                cGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAg                                                                                                                  >  2:1692941/1‑148 (MQ=255)
                                cGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAg                                                                                                                  <  1:97341/148‑1 (MQ=255)
                                   cGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATa                                                                                                               <  1:335446/148‑1 (MQ=255)
                                   cGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATa                                                                                                               <  2:1559618/148‑1 (MQ=255)
                                     gACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAAt                                                                                                             <  1:547735/148‑1 (MQ=255)
                                      aCTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAAtt                                                                                                            <  1:3577695/148‑1 (MQ=255)
                                         ataACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTg                                                                                                         >  1:2139365/1‑148 (MQ=255)
                                            aCCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAg                                                                                                      >  2:3655889/1‑148 (MQ=255)
                                             ccACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAgg                                                                                                     >  2:3470003/1‑148 (MQ=255)
                                                  aTCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCAt                                                                                                >  2:1524565/1‑148 (MQ=255)
                                                    cGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCAtct                                                                                              >  1:1300866/1‑148 (MQ=255)
                                                     gCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCAtctt                                                                                             <  2:567716/148‑1 (MQ=255)
                                                     gCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCAtctt                                                                                             <  1:2416973/148‑1 (MQ=255)
                                                     gCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCAtctt                                                                                             <  1:3224717/148‑1 (MQ=255)
                                                         tCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGa                                                                                         <  1:200185/148‑1 (MQ=255)
                                                         tCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGa                                                                                         >  2:328924/1‑148 (MQ=255)
                                                         tCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGa                                                                                         <  1:3270599/148‑1 (MQ=255)
                                                          cGGCCCGTGTCTTGCGCTCGGCAGCATGATAACGCTGGGGTATCTGGCGTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATACTTTTGCAGGCCCATCTTCTGaa                                                                                        <  2:3681687/148‑1 (MQ=255)
                                                          cGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGaa                                                                                        <  1:3063248/148‑1 (MQ=255)
                                                            gCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATc                                                                                      >  1:3210038/1‑148 (MQ=255)
                                                            gCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATc                                                                                      >  1:1241603/1‑148 (MQ=255)
                                                               cGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGTCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGa                                                                                   >  1:498476/1‑148 (MQ=255)
                                                               cGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGa                                                                                   <  2:1219768/148‑1 (MQ=255)
                                                               cGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCACGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGa                                                                                   >  2:405259/1‑148 (MQ=255)
                                                                   gCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCaa                                                                               >  2:2319209/1‑148 (MQ=255)
                                                                              gcagcaTGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTcc                                                                    >  2:3251749/1‑148 (MQ=255)
                                                                              gcagcaTGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTcc                                                                    <  2:2916975/148‑1 (MQ=255)
                                                                               cagcaTGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCAGTTTcca                                                                   >  1:2640462/1‑148 (MQ=255)
                                                                                  caTGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTATCCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTccagcc                                                                >  2:1324560/1‑148 (MQ=255)
                                                                                  caTGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTccagcc                                                                >  1:3066794/1‑148 (MQ=255)
                                                                                    tGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTccagccag                                                              >  1:452776/1‑148 (MQ=255)
                                                                                     gATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGt                                                             >  2:338554/1‑148 (MQ=255)
                                                                                      aTAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCACATCTTCTGAATCAGAGCAAGTTCCGTTTCCAGGCAGTc                                                            >  2:111321/1‑148 (MQ=255)
                                                                                         aaCCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGat                                                         <  1:527769/148‑1 (MQ=255)
                                                                                          aCCTGTGGTATCTGTCCTGGATCTCTTATTAACCGGGAGGGTTCTGCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCAATCTTCTGTATCAGATCAAGTTCCGTTTCCAGCCAGTCGtta                                                        <  1:194537/148‑1 (MQ=255)
                                                                                          aCCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGata                                                        >  2:2347312/1‑148 (MQ=255)
                                                                                          aCCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGata                                                        <  1:1067381/148‑1 (MQ=255)
                                                                                          aCCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGata                                                        <  1:1574511/148‑1 (MQ=255)
                                                                                          aCCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGata                                                        >  1:3773793/1‑148 (MQ=255)
                                                                                          aCCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGata                                                        >  1:579559/1‑148 (MQ=255)
                                                                                           ccTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGatat                                                       <  1:2976982/148‑1 (MQ=255)
                                                                                             tGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATgg                                                     >  2:105929/1‑148 (MQ=255)
                                                                                              ggggTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGc                                                    <  1:3664743/148‑1 (MQ=255)
                                                                                               gggTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGcc                                                   <  2:250400/148‑1 (MQ=255)
                                                                                               gggTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTGTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGcc                                                   <  2:255460/148‑1 (MQ=255)
                                                                                                    tCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCcttctt                                              <  2:3189342/148‑1 (MQ=255)
                                                                                                    tCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCcttctt                                              <  2:2833530/148‑1 (MQ=255)
                                                                                                     cTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAATTTCCGTTTCCAGCCAGTCGATATGGCCttcttc                                             >  1:2246615/1‑148 (MQ=255)
                                                                                                     cTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCttcttc                                             >  2:2241882/1‑148 (MQ=255)
                                                                                                       ggCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTtcat                                           <  1:1338330/148‑1 (MQ=255)
                                                                                                        gCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTtcatc                                          >  1:2365335/1‑148 (MQ=255)
                                                                                                         ccTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTtcatca                                         <  1:67478/148‑1 (MQ=255)
                                                                                                         ccTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTtcatca                                         <  2:3520796/148‑1 (MQ=255)
                                                                                                            ggATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGc                                      >  2:1981305/1‑148 (MQ=255)
                                                                                                            ggATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGc                                      >  2:2402340/1‑148 (MQ=255)
                                                                                                              aTCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCaa                                    >  1:2439225/1‑148 (MQ=255)
                                                                                                                ctctTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCaaaa                                  <  2:3724425/148‑1 (MQ=255)
                                                                                                                  ctTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAAtt                                >  2:1683425/1‑148 (MQ=255)
                                                                                                                    tattaACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTc                              <  2:1711384/148‑1 (MQ=255)
                                                                                                                      ttaACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTa                            <  2:3768055/148‑1 (MQ=255)
                                                                                                                      ttaACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTCCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTa                            <  1:886935/148‑1 (MQ=255)
                                                                                                                       taACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTat                           >  1:3395872/1‑148 (MQ=255)
                                                                                                                         aCCGGGAGGGTTCTCCCTCCCGTCTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTatca                         <  1:14594/148‑1 (MQ=255)
                                                                                                                         aCCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTatca                         >  2:700796/1‑148 (MQ=255)
                                                                                                                         aCCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTatca                         >  2:1345533/1‑148 (MQ=255)
                                                                                                                          ccGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTatcat                        >  1:2340074/1‑148 (MQ=255)
                                                                                                                           cGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTatcatc                       <  1:691865/148‑1 (MQ=255)
                                                                                                                            gggAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTAtcatca                      <  2:1616739/148‑1 (MQ=255)
                                                                                                                             ggAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATcatcat                     <  1:1613523/148‑1 (MQ=255)
                                                                                                                             ggAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATcatcat                     <  1:186051/148‑1 (MQ=255)
                                                                                                                             ggAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATcatcat                     >  2:1452094/1‑148 (MQ=255)
                                                                                                                              gAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCata                    <  2:1871556/148‑1 (MQ=255)
                                                                                                                              gAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCata                    <  2:452776/148‑1 (MQ=255)
                                                                                                                                gggTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATc                  >  1:1374134/1‑148 (MQ=255)
                                                                                                                                   ttCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATcgcg               >  1:1805316/1‑148 (MQ=255)
                                                                                                                                     ctccctccCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGc             <  1:1692941/148‑1 (MQ=255)
                                                                                                                                     ctccctccCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGc             <  1:1692079/148‑1 (MQ=255)
                                                                                                                                      tccctcccGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCt            >  2:3077894/1‑148 (MQ=255)
                                                                                                                                         ctcccGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGAc         <  2:1656790/148‑1 (MQ=255)
                                                                                                                                           cccGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGt       >  2:1708414/1‑148 (MQ=255)
                                                                                                                                            ccGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTa      <  2:765597/148‑1 (MQ=255)
                                                                                                                                                cacGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCa  <  1:912214/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TTGGTTACGCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCA  >  CP009273/3459458‑3459749

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: