Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 271,650 G→T 10.9% P205P (CCG→CCT mmuP → CP4‑6 prophage; putative S‑methylmethionine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273271,6500GT10.9% 598.1 / 30.7 179P205P (CCG→CCTmmuPCP4‑6 prophage; putative S‑methylmethionine transporter
Reads supporting (aligned to +/- strand):  ref base G (76/82);  new base T (7/13);  total (83/96)
Fisher's exact test for biased strand distribution p-value = 3.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 179 total aligned reads displayed.

TCGCTGGTCAAAGTGGTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGATTATTTTCTTTATCG  >  CP009273/271504‑271795
                                                                                                                                                  |                                                                                                                                                 
tcGCTGGTCAAAGTGGTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGa                                                                                                                                                  <  2:1706332/148‑1 (MQ=255)
tcGCTGGTCAAAGTGGTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGa                                                                                                                                                  <  1:1326487/148‑1 (MQ=255)
tcGCTGGTCAAAGTGGTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGa                                                                                                                                                  <  2:3599169/148‑1 (MQ=255)
     ggTCAAAGTGGTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTtg                                                                                                                                             <  1:2920191/148‑1 (MQ=255)
       tCAAAGTGGTCACTATCATCGCCTTTCTCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTtgat                                                                                                                                           <  2:1646337/148‑1 (MQ=255)
       tCAAAGTGGTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGCTGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTtgat                                                                                                                                           >  2:978138/1‑148 (MQ=255)
       tCAAAGTGGTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTtgat                                                                                                                                           >  2:1833271/1‑148 (MQ=255)
        cAAAGTGGTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCTATTTtgatg                                                                                                                                          >  2:1911984/1‑148 (MQ=255)
         aaaGTGGTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTtgatga                                                                                                                                         <  1:1605439/148‑1 (MQ=255)
         aaaGTGGTCACTATCATCGCCTTTATCATCCTCGCTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTtgatga                                                                                                                                         <  1:2854296/148‑1 (MQ=255)
           aGTGGTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACt                                                                                                                                       >  2:3101731/1‑148 (MQ=255)
           aGTGGTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACt                                                                                                                                       >  1:2405191/1‑148 (MQ=255)
           aGTGGTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACt                                                                                                                                       >  2:143462/1‑148 (MQ=255)
            gTGGTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTa                                                                                                                                      <  2:2991302/148‑1 (MQ=255)
               gTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTAtgg                                                                                                                                   <  2:1876501/148‑1 (MQ=255)
               gTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTAtgg                                                                                                                                   <  2:2169657/148‑1 (MQ=255)
                   cTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGc                                                                                                                               >  1:3488698/1‑148 (MQ=255)
                       catcGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATTCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTACCGCACGGTGGCTTACCGATTTTTATGACTATGGTTGCAGtt                                                                                                                           >  1:2612485/1‑147 (MQ=255)
                       catcGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTg                                                                                                                           >  1:3040278/1‑148 (MQ=255)
                         tcGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGaa                                                                                                                         <  2:628267/148‑1 (MQ=255)
                         tcGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGaa                                                                                                                         >  1:293950/1‑148 (MQ=255)
                          cGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACGATGGTGGCAGTGAAc                                                                                                                        >  1:20974/1‑148 (MQ=255)
                           gCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACt                                                                                                                       >  1:187860/1‑148 (MQ=255)
                           gCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACt                                                                                                                       <  2:3673407/148‑1 (MQ=255)
                            ccTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACtt                                                                                                                      >  2:1504202/1‑148 (MQ=255)
                            ccTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACtt                                                                                                                      >  2:2513130/1‑148 (MQ=255)
                            ccTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACtt                                                                                                                      >  2:2513245/1‑148 (MQ=255)
                              tttATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTg                                                                                                                    <  1:720304/148‑1 (MQ=255)
                                tATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCt                                                                                                                  <  1:962571/148‑1 (MQ=255)
                                 atcatcCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCtt                                                                                                                 <  2:2550673/148‑1 (MQ=255)
                                  tcatcCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCttt                                                                                                                >  1:1188956/1‑148 (MQ=255)
                                  tcatcCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCttt                                                                                                                <  1:1447229/148‑1 (MQ=255)
                                   catcCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCtttt                                                                                                               <  1:3638933/148‑1 (MQ=255)
                                     tcCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTc                                                                                                             >  1:3281446/1‑148 (MQ=255)
                                      ccTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCTATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCg                                                                                                            >  2:1713270/1‑148 (MQ=255)
                                      ccTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCg                                                                                                            >  1:3770771/1‑148 (MQ=255)
                                         cGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGt                                                                                                         >  2:1686759/1‑148 (MQ=255)
                                                gcggcgATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGc                                                                                                  <  1:3719350/148‑1 (MQ=255)
                                                  ggcgATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCTATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCtt                                                                                                <  2:2680171/148‑1 (MQ=255)
                                                  ggcgATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCtt                                                                                                <  1:2088242/148‑1 (MQ=255)
                                                  ggcgATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCtt                                                                                                <  1:470825/148‑1 (MQ=255)
                                                  ggcgATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCtt                                                                                                <  1:3527092/148‑1 (MQ=255)
                                                   gcgATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTa                                                                                               <  2:1188956/148‑1 (MQ=255)
                                                    cgATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTAt                                                                                              <  2:582028/148‑1 (MQ=255)
                                                      aTTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCg                                                                                            >  1:3385623/1‑148 (MQ=255)
                                                      aTTTTCGGCTTTATTCCGATGAAGGATGGCTCGCCAGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCg                                                                                            >  2:3042057/1‑148 (MQ=255)
                                                       ttttCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCgg                                                                                           >  2:2697986/1‑148 (MQ=255)
                                                         ttCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCa                                                                                         >  2:1012763/1‑148 (MQ=255)
                                                            ggCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTg                                                                                      <  2:3548305/148‑1 (MQ=255)
                                                              cTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTgcc                                                                                    >  1:1366521/1‑148 (MQ=255)
                                                               tttATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTgccg                                                                                   <  1:235766/148‑1 (MQ=255)
                                                               tttAGTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTgccg                                                                                   <  1:307504/148‑1 (MQ=255)
                                                                    tCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGt                                                                              >  1:404151/1‑148 (MQ=255)
                                                                       gATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCTATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGaa                                                                           >  1:2153943/1‑148 (MQ=255)
                                                                        aTGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGaaa                                                                          <  1:3215612/148‑1 (MQ=255)
                                                                             ggATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGaa                                                                     >  1:1584020/1‑148 (MQ=255)
                                                                              gATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGaaa                                                                    >  2:436783/1‑148 (MQ=255)
                                                                                         cgcgCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCaaa                                                         <  2:624952/148‑1 (MQ=255)
                                                                                           cgcCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGt                                                       <  1:495584/148‑1 (MQ=255)
                                                                                            gcCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGtt                                                      <  1:2113184/148‑1 (MQ=255)
                                                                                             ccGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTa                                                     <  2:3281446/148‑1 (MQ=255)
                                                                                                gggCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATcc                                                  <  2:3406631/148‑1 (MQ=255)
                                                                                                 ggCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATccc                                                 >  2:2389788/1‑148 (MQ=255)
                                                                                                      aGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAg                                            >  1:915719/1‑148 (MQ=255)
                                                                                                      aGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAg                                            >  1:11549/1‑148 (MQ=255)
                                                                                                      aGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAg                                            >  1:30173/1‑148 (MQ=255)
                                                                                                           tatCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTAGCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGAtt                                       <  1:1843794/148‑1 (MQ=255)
                                                                                                           tatCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGagt                                       >  1:928739/1‑146 (MQ=255)
                                                                                                           tatCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGAtt                                       <  1:3345515/148‑1 (MQ=255)
                                                                                                            atCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTc                                      >  2:3815108/1‑148 (MQ=255)
                                                                                                             tCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCTATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCg                                     <  2:2339241/148‑1 (MQ=255)
                                                                                                             tCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCTATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCg                                     <  2:417528/148‑1 (MQ=255)
                                                                                                                    aGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTTTGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTAcc                              <  1:1646839/148‑1 (MQ=255)
                                                                                                                    aGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTAcc                              <  1:3557202/148‑1 (MQ=255)
                                                                                                                    aGAAGGCTGGGTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTAcc                              <  1:2049231/148‑1 (MQ=255)
                                                                                                                      aaGGCTGGTTCCCGCACGGTGGCTTACCTATTTTGATGAATATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCAt                            >  1:3626582/1‑148 (MQ=255)
                                                                                                                      aaGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCAt                            >  2:2845996/1‑148 (MQ=255)
                                                                                                                         gCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCAATGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATcgc                         >  1:2791220/1‑148 (MQ=255)
                                                                                                                            ggTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATcgcgcg                      >  2:1538701/1‑148 (MQ=255)
                                                                                                                            ggTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATcgcgcg                      >  1:2752833/1‑148 (MQ=255)
                                                                                                                            ggTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATcgcgcg                      >  2:2615961/1‑148 (MQ=255)
                                                                                                                             gTTCCCGCACGGTGGCTTACCTATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGa                     <  1:2650902/148‑1 (MQ=255)
                                                                                                                             gTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGa                     <  1:1035878/148‑1 (MQ=255)
                                                                                                                              gtCCCGCACGGTGGCTTACCGATTGTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATcgcgcgtc                    <  1:1096631/147‑3 (MQ=255)
                                                                                                                                  cGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGat                <  2:738843/148‑1 (MQ=255)
                                                                                                                                   gCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGatt               <  1:2167320/148‑1 (MQ=255)
                                                                                                                                   gCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGatt               <  1:2167252/148‑1 (MQ=255)
                                                                                                                                    cACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGatta              >  2:3227236/1‑148 (MQ=255)
                                                                                                                                    cACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGatta              <  1:2513245/148‑1 (MQ=255)
                                                                                                                                    cACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGatta              <  1:1504202/148‑1 (MQ=255)
                                                                                                                                     aCGGTGGCTTACCTATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGattat             <  1:9792/148‑1 (MQ=255)
                                                                                                                                         tGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGATTATTTTc         >  1:2813636/1‑148 (MQ=255)
                                                                                                                                         tGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGATTATTTTc         >  2:2628139/1‑148 (MQ=255)
                                                                                                                                          ggCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGATTATTTTCt        <  2:2612485/148‑1 (MQ=255)
                                                                                                                                           gCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGATTATTTTCtt       >  2:3213703/1‑148 (MQ=255)
                                                                                                                                            cTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGATTATTTTCttt      <  1:1761097/148‑1 (MQ=255)
                                                                                                                                             ttACCTATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGTAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGATTATTTTCTTTa     <  2:2144030/148‑1 (MQ=255)
                                                                                                                                             ttACCTATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGATTATTTTCTTTa     <  2:2715955/148‑1 (MQ=255)
                                                                                                                                             ttACCTATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGATTATTTTCTTTa     <  2:2144025/148‑1 (MQ=255)
                                                                                                                                                ccGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGATTATTTTCTTTATCg  >  2:3720785/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TCGCTGGTCAAAGTGGTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGATTATTTTCTTTATCG  >  CP009273/271504‑271795

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: