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% 576.7 / 28.4 173P205P (CCG→CCTmmuPCP4‑6 prophage; putative S‑methylmethionine transporter
Reads supporting (aligned to +/- strand):  ref base G (80/73);  new base T (9/10);  total (90/83)
Fisher's exact test for biased strand distribution p-value = 8.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 180 total aligned reads displayed.

CTCGCTGGTCAAAGTGGTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGATTATTTTCTTTATCGGC  >  CP009273/271503‑271797
                                                                                                                                                   |                                                                                                                                                   
  cGCTGGTCAAAGTGGTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCTAt                                                                                                                                                   >  2:3472773/1‑148 (MQ=255)
   gCTGGTCAAAGTGGTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGAtt                                                                                                                                                  <  2:2130654/148‑1 (MQ=255)
    cTGGTCAAAGTGGTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGAttt                                                                                                                                                 >  1:1934215/1‑148 (MQ=255)
      ggTCAAAGTGGTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATGTTCGGCTTTATTCCGATGCCCGATGGCTCGCCCGCGCCGGGGCTGAGTAATATAACGGCAGAAGGCTGGTTCCCGCACGGTGGCGTACCGATTTtg                                                                                                                                               >  1:711362/1‑148 (MQ=255)
          aaaGTGGTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTtgatga                                                                                                                                           >  2:1697763/1‑148 (MQ=255)
          aaaGTGGTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTtgatga                                                                                                                                           <  1:3349419/148‑1 (MQ=255)
           aaGTGGTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGAc                                                                                                                                          >  2:2407012/1‑148 (MQ=255)
             gTGGTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTa                                                                                                                                        >  1:1134783/1‑148 (MQ=255)
              tGGTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTAt                                                                                                                                       >  2:3090510/1‑148 (MQ=255)
               ggTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTAtg                                                                                                                                      <  2:2205392/148‑1 (MQ=255)
               ggTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTAtg                                                                                                                                      <  1:3051401/148‑1 (MQ=255)
                gTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGTAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTAtgg                                                                                                                                     <  2:2458831/148‑1 (MQ=255)
                gTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTAtgg                                                                                                                                     <  2:1851672/148‑1 (MQ=255)
                gTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTAtgg                                                                                                                                     <  2:47599/148‑1 (MQ=255)
                  cACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGGATTATTACGATTCAGGATGGCTCGCCCCCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTACCGCACGGTGGCTTAACAAATTTTATGAATAttgtt                                                                                                                                   >  1:3388950/1‑147 (MQ=255)
                  cACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTAtggtg                                                                                                                                   <  1:2274427/148‑1 (MQ=255)
                    cTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGc                                                                                                                                 >  1:2077426/1‑148 (MQ=255)
                     tATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCTATTTTGATGACTATGGTGGCa                                                                                                                                <  2:3187813/148‑1 (MQ=255)
                     tATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCa                                                                                                                                >  2:603617/1‑148 (MQ=255)
                     tATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCa                                                                                                                                >  1:3297498/1‑148 (MQ=255)
                      atcatcGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAg                                                                                                                               <  1:222208/148‑1 (MQ=255)
                           cGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAAc                                                                                                                          <  2:826326/148‑1 (MQ=255)
                           cGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAAc                                                                                                                          >  2:167178/1‑148 (MQ=255)
                            gCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACt                                                                                                                         <  1:603547/148‑1 (MQ=255)
                                ttATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGc                                                                                                                     >  1:3044990/1‑148 (MQ=255)
                                 tATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCt                                                                                                                    <  1:1097955/148‑1 (MQ=255)
                                 tATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCt                                                                                                                    <  2:1820747/148‑1 (MQ=255)
                                  atcatcCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCtt                                                                                                                   >  1:2221298/1‑148 (MQ=255)
                                   tcatcCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCTATTTTGATGACTATGGTGGCAGTGAACTTTGCttt                                                                                                                  <  2:1983415/148‑1 (MQ=255)
                                      tcCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTc                                                                                                               >  2:195286/1‑148 (MQ=255)
                                      tcCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTc                                                                                                               >  2:1594150/1‑148 (MQ=255)
                                      tcCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTc                                                                                                               >  2:1593087/1‑148 (MQ=255)
                                        cTCGGTGGGGCGGCGATTTTCGGCTTTAATCCGATGCAGGATGGCACGCCCGCGCCGGGGCTGAGTAATATAACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCgg                                                                                                             >  1:2084105/1‑148 (MQ=255)
                                         tCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCggg                                                                                                            >  1:1963160/1‑148 (MQ=255)
                                         tCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCggg                                                                                                            >  1:3431633/1‑148 (MQ=255)
                                             tGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTAcc                                                                                                        <  1:2850838/148‑1 (MQ=255)
                                             tGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTAcc                                                                                                        >  1:2511327/1‑148 (MQ=255)
                                              ggggCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCg                                                                                                       <  2:2630157/148‑1 (MQ=255)
                                               gggCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGa                                                                                                      >  2:235177/1‑148 (MQ=255)
                                                ggcggcGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAg                                                                                                     <  2:2634507/148‑1 (MQ=255)
                                                  cggcgATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCt                                                                                                   <  1:2530368/148‑1 (MQ=255)
                                                    gcgATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTa                                                                                                 >  2:1808681/1‑148 (MQ=255)
                                                     cgATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTAt                                                                                                >  1:455975/1‑148 (MQ=255)
                                                     cgATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTAt                                                                                                >  2:2697810/1‑148 (MQ=255)
                                                       aTTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATct                                                                                              >  1:3366455/1‑147 (MQ=255)
                                                       aTTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATATTGATCACTATGGTGGCAGTGAACTTTGCATTTTCGGGTACCGAGCTAATCg                                                                                              >  2:721359/1‑148 (MQ=255)
                                                        ttttCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCgg                                                                                             <  1:352296/148‑1 (MQ=255)
                                                         tttCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGc                                                                                            <  2:2853862/148‑1 (MQ=255)
                                                          ttCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCTATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCa                                                                                           <  1:1757904/148‑1 (MQ=255)
                                                          ttCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCa                                                                                           <  2:3296892/148‑1 (MQ=255)
                                                          ttCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCa                                                                                           <  1:2754529/148‑1 (MQ=255)
                                                          ttCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCa                                                                                           <  2:2985711/148‑1 (MQ=255)
                                                          ttCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCa                                                                                           <  2:3297287/148‑1 (MQ=255)
                                                           tCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCAt                                                                                          <  2:1232281/148‑1 (MQ=255)
                                                              gCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTgc                                                                                       >  2:1001156/1‑148 (MQ=255)
                                                              gCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTgc                                                                                       >  1:451695/1‑148 (MQ=255)
                                                               cTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACNGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTgcc                                                                                      >  1:1805724/1‑148 (MQ=255)
                                                                tttATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTgccg                                                                                     >  1:1224750/1‑148 (MQ=255)
                                                                  tATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTgccgcc                                                                                   >  2:3457952/1‑148 (MQ=255)
                                                                     tCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGt                                                                                <  1:714222/148‑1 (MQ=255)
                                                                       cGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCTATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGa                                                                              >  2:603181/1‑148 (MQ=255)
                                                                          tGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAAc                                                                           <  2:3047613/148‑1 (MQ=255)
                                                                            cAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACgg                                                                         <  2:1934215/148‑1 (MQ=255)
                                                                              ggATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGaa                                                                       >  2:1015431/1‑148 (MQ=255)
                                                                               gATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGaaa                                                                      >  1:1954622/1‑148 (MQ=255)
                                                                                aTGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGaaaa                                                                     >  1:2369931/1‑148 (MQ=255)
                                                                                  ggCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAAcc                                                                   >  2:1667073/1‑148 (MQ=255)
                                                                                     tCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCcgc                                                                <  2:1281935/148‑1 (MQ=255)
                                                                                         ccGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCaa                                                            <  1:2597338/148‑1 (MQ=255)
                                                                                            cgcCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGt                                                         >  1:1617285/1‑148 (MQ=255)
                                                                                              ccGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTa                                                       <  1:138208/148‑1 (MQ=255)
                                                                                              ccGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTa                                                       <  2:3016447/148‑1 (MQ=255)
                                                                                                ggggCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATc                                                     <  2:208929/148‑1 (MQ=255)
                                                                                                 gggCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGTAAACCCGCGCAAAGTTATcc                                                    <  2:3560535/148‑1 (MQ=255)
                                                                                                         tAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCg                                            >  1:1518279/1‑148 (MQ=255)
                                                                                                           atatCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGTCAGTGAACTTTGCTTTTTCGGGTACCGAGATTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCTGTAGCGAt                                          <  2:2221264/148‑1 (MQ=255)
                                                                                                           atatCACGGCAGAAGGCTGGGTCCCGCACGGTGGCTTACCTATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGAt                                          <  1:2360660/148‑1 (MQ=255)
                                                                                                               cACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGt                                      >  2:3379430/1‑148 (MQ=255)
                                                                                                                 cGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGtac                                    >  2:627321/1‑148 (MQ=255)
                                                                                                                   gCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGtacta                                  <  1:3090510/148‑1 (MQ=255)
                                                                                                                     aGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTAcc                                <  1:1455215/148‑1 (MQ=255)
                                                                                                                       aaGGCTGGTTCCCGCACGGTGGCTTACCTATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCAt                              >  1:6341/1‑148 (MQ=255)
                                                                                                                        aGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATc                             >  1:2119920/1‑148 (MQ=255)
                                                                                                                           cTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATcgcg                          >  2:2286564/1‑148 (MQ=255)
                                                                                                                              gTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCTATTCGTACTACCATCGCGCGa                       <  2:3388950/148‑1 (MQ=255)
                                                                                                                              gTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCAGTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGa                       <  2:1813776/148‑1 (MQ=255)
                                                                                                                                 cccGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTg                    >  1:557874/1‑148 (MQ=255)
                                                                                                                                 cccGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTg                    >  1:1590713/1‑148 (MQ=255)
                                                                                                                                     cACGGTGGCTTACCTATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGatta                >  1:1109265/1‑148 (MQ=255)
                                                                                                                                     cACGGTGGCTTACCTATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGatta                >  1:1888135/1‑148 (MQ=255)
                                                                                                                                      aCGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGattat               <  1:1624784/148‑1 (MQ=255)
                                                                                                                                      aCGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGattat               <  2:1028220/148‑1 (MQ=255)
                                                                                                                                         gTGGCTTACCTATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGATTAtttt            <  1:2002162/148‑1 (MQ=255)
                                                                                                                                            gCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGATTATTTTCtt         >  2:2758714/1‑148 (MQ=255)
                                                                                                                                               tACCTATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGATTATTTTCTTTAt      >  2:643901/1‑148 (MQ=255)
                                                                                                                                                aCCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGATTATTTTCTTTATc     >  1:891690/1‑148 (MQ=255)
                                                                                                                                                 ccGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGATTATTTTCTTTATCg    >  2:1406899/1‑148 (MQ=255)
                                                                                                                                                   gATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGATTATTTTCTTTATCGGc  <  1:507334/148‑1 (MQ=255)
                                                                                                                                                   gATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGATTATTTTCTTTATCGGc  <  2:711362/148‑1 (MQ=255)
                                                                                                                                                   gATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGATTATTTTCTTTATCGGc  <  1:2697810/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CTCGCTGGTCAAAGTGGTCACTATCATCGCCTTTATCATCCTCGGTGGGGCGGCGATTTTCGGCTTTATTCCGATGCAGGATGGCTCGCCCGCGCCGGGGCTGAGTAATATCACGGCAGAAGGCTGGTTCCCGCACGGTGGCTTACCGATTTTGATGACTATGGTGGCAGTGAACTTTGCTTTTTCGGGTACCGAGCTTATCGGCATTGCCGCCGGTGAAACGGAAAACCCGCGCAAAGTTATCCCGGTAGCGATTCGTACTACCATCGCGCGACTGATTATTTTCTTTATCGGC  >  CP009273/271503‑271797

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