Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,626,583 G→T 7.3% intergenic (‑41/‑746) ydfJ ← / → ydfK pseudogene, MFS transporter family; interrupted by Qin prophage;Phage or Prophage Related; putative transport protein/cold shock protein, function unknown, Qin prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,626,5830GT 7.3% 911.9 / 25.7 254intergenic (‑41/‑746)ydfJ/ydfKpseudogene, MFS transporter family; interrupted by Qin prophage;Phage or Prophage Related; putative transport protein/cold shock protein, function unknown, Qin prophage
Reads supporting (aligned to +/- strand):  ref base G (111/124);  new base T (9/10);  total (120/134)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 256 total aligned reads displayed.

CCGTAGGTTCCCATTGCCAGAATTAACGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATGGATTTATGC  >  CP009273/1626436‑1626726
                                                                                                                                                   |                                                                                                                                               
 cGTAGGTTCCCATTGCCAGAATTAACGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCgg                                                                                                                                                >  2:2323250/1‑147 (MQ=255)
 cGTAGGTTCCCATTGCCAGAATTAACGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGt                                                                                                                                                <  2:623057/148‑1 (MQ=255)
   tAGGTTCCCATTGCCAGAATTAACGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGttt                                                                                                                                              >  1:2224675/1‑148 (MQ=255)
    aGGTTCCCATTGCCAGAATTAACGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGtttt                                                                                                                                             <  1:3216572/148‑1 (MQ=255)
    aGGTTCCCATTGCCAGAATTAACGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGtttt                                                                                                                                             >  1:792076/1‑148 (MQ=255)
    aGGTTCCCATTGCCAGAATTAACGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGtttt                                                                                                                                             <  2:1373621/148‑1 (MQ=255)
     ggTTCCCATTGCCAGAATTAACGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTc                                                                                                                                            >  1:3163547/1‑148 (MQ=255)
     ggTTCCCATTGCCAGAATTAACGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTc                                                                                                                                            >  1:1337565/1‑148 (MQ=255)
     ggTTCCCATTGCCAGAATTAACGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTc                                                                                                                                            >  1:26529/1‑148 (MQ=255)
      gTTCCCATTGCCAGAATTAACGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATACTTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCt                                                                                                                                           <  2:2197990/148‑1 (MQ=255)
      gTTCCCATTGCCAGAATTAACGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAAATCACAGGTGTTTTTTCCCATCCTGTGGGTTCCTTGGCGTTTTCt                                                                                                                                           >  2:3559860/1‑148 (MQ=255)
         cccATTGCCAGAATTAACGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAgg                                                                                                                                        <  2:2239355/148‑1 (MQ=255)
           cagtGCCAGAATTAACGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGtt                                                                                                                                      <  1:2940432/145‑1 (MQ=255)
            aTTGCCAGAATTAACGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCTAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGttt                                                                                                                                     <  1:672240/148‑1 (MQ=255)
            aTTGCCAGAATTAACGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGttt                                                                                                                                     >  2:3277941/1‑148 (MQ=255)
             ttGCCAGAATTAACGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGtttt                                                                                                                                    <  1:3329219/148‑1 (MQ=255)
              tGCCAGAATTAACGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGttttt                                                                                                                                   >  1:96514/1‑148 (MQ=255)
                ccAGAATTAACGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTc                                                                                                                                 <  1:3305655/148‑1 (MQ=255)
                 cAGAATTAACGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCa                                                                                                                                >  1:1397939/1‑148 (MQ=255)
                  aGAATTAACGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTGTAGGTTTTTTCAg                                                                                                                               >  1:1733049/1‑148 (MQ=255)
                     aTTAACGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATa                                                                                                                            <  2:3701666/148‑1 (MQ=255)
                         aCGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATATTTg                                                                                                                        >  1:544226/1‑148 (MQ=255)
                          cGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTGAGATAGtttg                                                                                                                       >  1:1946338/1‑146 (MQ=255)
                          cGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGc                                                                                                                       >  1:2189287/1‑148 (MQ=255)
                           gCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCa                                                                                                                      <  2:3456447/148‑1 (MQ=255)
                             cATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTGTAGGTTTTTTCAGATAGTTGCAtt                                                                                                                    >  1:2763822/1‑148 (MQ=255)
                             cATTGCAGTTGATGATTCATGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCAtt                                                                                                                    <  2:3009585/148‑1 (MQ=255)
                              aTTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCAttt                                                                                                                   >  2:1386813/1‑148 (MQ=255)
                              aTTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCAttt                                                                                                                   <  1:2291229/148‑1 (MQ=255)
                               ttGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCAtttt                                                                                                                  <  2:1593875/148‑1 (MQ=255)
                                   cGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGAGAATATAACTGTAAATCCATAATTCACATGTGTTTTTTCCCATCCTTTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTa                                                                                                              <  1:3663016/148‑1 (MQ=255)
                                    gTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTaa                                                                                                             >  1:831741/1‑148 (MQ=255)
                                    gTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTaa                                                                                                             <  1:1182907/148‑1 (MQ=255)
                                      tgatgaTTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTaaaa                                                                                                           <  1:3594294/148‑1 (MQ=255)
                                        atgatTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAg                                                                                                         >  2:675651/1‑148 (MQ=255)
                                        atgatTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAg                                                                                                         >  2:169698/1‑148 (MQ=255)
                                             tCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATcc                                                                                                    <  2:663121/148‑1 (MQ=255)
                                                ggAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTaa                                                                                                 <  1:856074/148‑1 (MQ=255)
                                                 gAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAg                                                                                                <  1:2715520/148‑1 (MQ=255)
                                                  aaaaaaTATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGt                                                                                               >  2:1542951/1‑148 (MQ=255)
                                                  aaaaaaTATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGt                                                                                               <  1:619881/148‑1 (MQ=255)
                                                   aaaaaTATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGtt                                                                                              <  2:1474992/148‑1 (MQ=255)
                                                       atatTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCTTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGAt                                                                                          >  2:2873898/1‑148 (MQ=255)
                                                         atTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTGTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATACTAAGTTAGAtct                                                                                        >  1:1650495/1‑148 (MQ=255)
                                                               tGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAgtgt                                                                                  <  1:290346/148‑1 (MQ=255)
                                                               tGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAgtgt                                                                                  <  1:1200473/148‑1 (MQ=255)
                                                                gAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTc                                                                                 <  1:3783824/148‑1 (MQ=255)
                                                                  aaCACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCCCAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTa                                                                               <  2:2714443/148‑1 (MQ=255)
                                                                   acacTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTAt                                                                              >  2:2799899/1‑148 (MQ=255)
                                                                   acacTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTAt                                                                              <  1:987696/148‑1 (MQ=255)
                                                                    cacTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCCGTGTCTATc                                                                             <  1:1849199/148‑1 (MQ=255)
                                                                    cacTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATc                                                                             <  1:2795170/148‑1 (MQ=255)
                                                                      cTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTg                                                                           <  2:1161601/148‑1 (MQ=255)
                                                                      cTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTg                                                                           <  2:1215757/148‑1 (MQ=255)
                                                                      cTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTg                                                                           >  2:3159959/1‑148 (MQ=255)
                                                                      cTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTg                                                                           <  1:2456056/148‑1 (MQ=255)
                                                                           gCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCCGTGTCTATCTGGGGcc                                                                      <  1:1702308/148‑1 (MQ=255)
                                                                            cTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCTTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCt                                                                     >  1:2299374/1‑148 (MQ=255)
                                                                             tGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTa                                                                    <  2:789512/148‑1 (MQ=255)
                                                                               cgcCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCCTCCTTAGTTCGTTCTCAGTGTCTATCTGGGGCCTAtt                                                                  >  1:3033123/1‑148 (MQ=255)
                                                                                gcCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCTTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTAttt                                                                 <  1:2049894/148‑1 (MQ=255)
                                                                                      cGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTcc                                                           <  2:1145103/148‑1 (MQ=255)
                                                                                        aaTATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCTTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCa                                                         <  2:763582/148‑1 (MQ=255)
                                                                                         atatAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCat                                                        >  2:1342697/1‑148 (MQ=255)
                                                                                          tataACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCata                                                       >  1:1521833/1‑148 (MQ=255)
                                                                                              aCTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATg                                                   >  2:3177164/1‑148 (MQ=255)
                                                                                               cTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGc                                                  <  2:1971347/148‑1 (MQ=255)
                                                                                                   aaaTCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCa                                              >  2:3007524/1‑148 (MQ=255)
                                                                                                     aTCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCATGGCGATTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCAACCAAAGTTCGATATCAGTGTCTATCTGGGGCCTATTTCTGTCCCAAAAATGCCCCaaa                                            >  2:3408335/1‑148 (MQ=255)
                                                                                                      tCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCTTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCaaaa                                           >  2:2548631/1‑148 (MQ=255)
                                                                                                      tCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAACTTAGATCTAAGTGTCTATCTGGGGCCTAATTCTGTCCCAAATATGCCCCaaaa                                           >  1:2039738/1‑148 (MQ=255)
                                                                                                       ccATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCaaaaa                                          >  2:952044/1‑148 (MQ=255)
                                                                                                       ccATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCaaaaa                                          >  2:591529/1‑148 (MQ=255)
                                                                                                         aTAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCaaaaaaa                                        >  2:2783034/1‑148 (MQ=255)
                                                                                                              tCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTccc                                   <  1:2363708/148‑1 (MQ=255)
                                                                                                                     tgtTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGa                            <  2:1733307/148‑1 (MQ=255)
                                                                                                                     tgtTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGa                            <  2:3513984/148‑1 (MQ=255)
                                                                                                                       ttttttCCCATCCTGTGGTTTCCTTGGCGTTTTCTATGTTTTTTAAGATAGTTGCATTTTTTTTAAAATCATCCTAAGTTAGATCTCAGTGTCTATATGGGGGCTATTTCTGTCCCATAAATTACCCAAAAAAAATACCCAACACATa                          >  1:1702200/1‑148 (MQ=255)
                                                                                                                         ttttCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATATTTGCATTTTTTTAAAAAGCATCCTAAGTTAGATCTCAGTGTCTATCTGGGGCCTATTTCTTTCCCATATATGCCCCAAAAAAACTACCCAACAGATAag                        >  1:1903754/1‑148 (MQ=255)
                                                                                                                         ttttCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCCGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAag                        <  1:1717982/148‑1 (MQ=255)
                                                                                                                           ttCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAagta                      <  2:2280451/148‑1 (MQ=255)
                                                                                                                              ccATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGtt                   >  2:1850117/1‑148 (MQ=255)
                                                                                                                              ccATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGtt                   <  1:187425/148‑1 (MQ=255)
                                                                                                                              ccATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGtt                   <  2:443173/148‑1 (MQ=255)
                                                                                                                               cATCCTGTGGTTTCCTTGGCGTTTTGTAGGTTTTTTCAGATAGTTTCATTTTTTTAAAAAACAACCAAAGTTAGATCTCAGTGTCTATAATGGGCCTATTTATGTCCCAAATATTCCCCAAAAAAAAACCCCCACACAAAATTAGttt                  >  1:1953696/1‑148 (MQ=255)
                                                                                                                                 tCCTGTGGTTTCCTTGGCTTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGttttt                >  2:2452493/1‑148 (MQ=255)
                                                                                                                                 tCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGttttt                <  1:1542951/148‑1 (MQ=255)
                                                                                                                                  ccTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATATTTTCATTTTTTTAAAAATCATCCTAAGTTCGATATCAGTGTCTATATGGGGCCTATTTATGTCCCAAATATTCCCCAAAAAAACTACCCAACAGATAAGTAGtttttt               >  1:2357254/1‑148 (MQ=255)
                                                                                                                                  ccTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTAAGATAGTTTCATTTTTTTAAAAAGCATCCTAAGTTCGATATCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGtttttt               >  1:3121315/1‑148 (MQ=255)
                                                                                                                                   cTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTc              <  1:405796/148‑1 (MQ=255)
                                                                                                                                   cTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATACTAAGTTCGATCTCAGTGTCTATATTGGGCCTATTTATGTCCCATATATTCCCCAAAAAAACACCCCAACACATAAGTAGTTTTTTc              >  1:3670645/1‑148 (MQ=255)
                                                                                                                                    tgtgGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCa             >  1:1273505/1‑148 (MQ=255)
                                                                                                                                    tgtgGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCa             <  2:290813/148‑1 (MQ=255)
                                                                                                                                      tgGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCAtt           >  2:3325530/1‑147 (MQ=255)
                                                                                                                                       ggTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATgg          <  1:1940736/148‑1 (MQ=255)
                                                                                                                                       ggTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATgg          <  1:1980487/148‑1 (MQ=255)
                                                                                                                                       ggTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATgg          <  2:1106699/148‑1 (MQ=255)
                                                                                                                                          ttCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCAAACAGATAAGTAGTTTTTTCATGGAtt       >  2:1965071/1‑148 (MQ=255)
                                                                                                                                              ttGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATGGATTTATg   <  1:3256361/148‑1 (MQ=255)
                                                                                                                                              ttGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACAACCCAACAGATAAGTAGTTTTTTAATGGATTTATg   >  1:1013991/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
CCGTAGGTTCCCATTGCCAGAATTAACGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATGGATTTATGC  >  CP009273/1626436‑1626726

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