Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,643,955 C→A 43.8% A30A (GCC→GCA dicB → Qin prophage; cell division inhibition protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,643,9550CA43.8% 254.0 / inf 306A30A (GCC→GCAdicBQin prophage; cell division inhibition protein
Reads supporting (aligned to +/- strand):  ref base C (90/82);  new base A (68/66);  total (158/148)
Fisher's exact test for biased strand distribution p-value = 8.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.15e-01

Only 100 of 309 total aligned reads displayed.

AACGTCACCTGGAGGCACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGC  >  CP009273/1643808‑1644102
                                                                                                                                                   |                                                                                                                                                   
aaCGTCACCTGGAGGCACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGcc                                                                                                                                                     <  2:311322/148‑1 (MQ=255)
    tCACCTGGAGGCACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTa                                                                                                                                                 <  2:1550745/148‑1 (MQ=255)
     cACCTGGAGGCACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTat                                                                                                                                                >  2:3428143/1‑147 (MQ=255)
     cACCTGGAGGCACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTaa                                                                                                                                                >  2:871145/1‑148 (MQ=255)
     cACCTGGAGGCACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTaa                                                                                                                                                >  2:1694396/1‑148 (MQ=255)
     cACCTGGAGGCACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTaa                                                                                                                                                >  2:1412248/1‑148 (MQ=255)
       ccTGGAGGCACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGTTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACa                                                                                                                                              <  2:1420218/148‑1 (MQ=255)
        cTGGAGGCACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACaa                                                                                                                                             <  2:2945911/148‑1 (MQ=255)
         tGGAGGCACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAag                                                                                                                                            >  1:2846125/1‑148 (MQ=255)
         tGGAGGCACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAag                                                                                                                                            <  1:2243161/148‑1 (MQ=255)
           gAGGCACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAagag                                                                                                                                          >  2:2259105/1‑148 (MQ=255)
            aGGCACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAgaga                                                                                                                                         <  1:1694396/148‑1 (MQ=255)
             ggCACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGaa                                                                                                                                        >  1:218376/1‑148 (MQ=255)
               cACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGaaaa                                                                                                                                      <  1:1378354/148‑1 (MQ=255)
               cACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGaaaa                                                                                                                                      >  2:735025/1‑148 (MQ=255)
                      cACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAAc                                                                                                                               <  2:3783045/148‑1 (MQ=255)
                      cACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAAc                                                                                                                               >  1:513146/1‑148 (MQ=255)
                         tGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACggg                                                                                                                            <  1:1666117/148‑1 (MQ=255)
                         tGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACggg                                                                                                                            <  1:3284336/148‑1 (MQ=255)
                         tGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACggg                                                                                                                            >  2:3600598/1‑148 (MQ=255)
                         tGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACggg                                                                                                                            <  1:3505995/148‑1 (MQ=255)
                          gCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGa                                                                                                                           >  2:1098996/1‑148 (MQ=255)
                             tCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCt                                                                                                                        <  2:2418399/148‑1 (MQ=255)
                                  aaaTTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATtaa                                                                                                                   >  2:1158942/1‑148 (MQ=255)
                                   aaTTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATtaaa                                                                                                                  >  2:838212/1‑148 (MQ=255)
                                   aaTTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATtaaa                                                                                                                  <  2:1228912/148‑1 (MQ=255)
                                     ttCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATtaaata                                                                                                                >  2:1277777/1‑148 (MQ=255)
                                        aTTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTaaataaaa                                                                                                             <  1:1847082/148‑1 (MQ=255)
                                         ttgttgAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTaaataaaat                                                                                                            <  2:1038435/148‑1 (MQ=255)
                                         ttgttgAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTaaataaaat                                                                                                            <  1:2092289/148‑1 (MQ=255)
                                          tgttgAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTaaataaaata                                                                                                           <  2:1001424/148‑1 (MQ=255)
                                          tgttgAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTaaataaaata                                                                                                           >  2:223462/1‑148 (MQ=255)
                                              gAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCa                                                                                                       <  1:38761/148‑1 (MQ=255)
                                               aGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCAt                                                                                                      <  2:1813117/148‑1 (MQ=255)
                                               aGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCAt                                                                                                      <  2:1718801/148‑1 (MQ=255)
                                                    gcgATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGttt                                                                                                 <  2:1308910/148‑1 (MQ=255)
                                                       aTAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCaa                                                                                              <  1:2041705/148‑1 (MQ=255)
                                                       aTAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCaa                                                                                              >  1:579549/1‑148 (MQ=255)
                                                        tAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAAt                                                                                             >  2:1924164/1‑148 (MQ=255)
                                                        tAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAAt                                                                                             >  1:654392/1‑148 (MQ=255)
                                                        tAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAAt                                                                                             >  2:28539/1‑148 (MQ=255)
                                                         aaTGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATg                                                                                            >  1:3633749/1‑148 (MQ=255)
                                                         aaTGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATg                                                                                            >  2:3235411/1‑148 (MQ=255)
                                                            gAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTg                                                                                         >  2:930053/1‑148 (MQ=255)
                                                              aaaCGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGc                                                                                       >  1:3389289/1‑148 (MQ=255)
                                                                aCGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTc                                                                                     >  1:3801233/1‑148 (MQ=255)
                                                                 cGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCg                                                                                    <  2:2765206/148‑1 (MQ=255)
                                                                  gTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGt                                                                                   <  2:2105102/148‑1 (MQ=255)
                                                                    tattaCCAAACGTTTATAAGTCTGTAGTTTGTTTTGATATTGGTGTCACTATCAGTAACCAAGTTTTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGttt                                                                                 <  2:3542088/148‑1 (MQ=255)
                                                                       taCCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACg                                                                              >  2:3012867/1‑148 (MQ=255)
                                                                           aaaCGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTg                                                                          <  1:310453/148‑1 (MQ=255)
                                                                           aaaCGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTg                                                                          >  1:3594888/1‑148 (MQ=255)
                                                                            aaCGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGa                                                                         >  2:1704432/1‑148 (MQ=255)
                                                                                ttAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGcc                                                                     >  1:1549137/1‑148 (MQ=255)
                                                                                ttAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGcc                                                                     >  1:3060328/1‑148 (MQ=255)
                                                                                 tAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCa                                                                    >  1:304365/1‑148 (MQ=255)
                                                                                   aTACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCaaa                                                                  >  2:3364979/1‑148 (MQ=255)
                                                                                   aTACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCaaa                                                                  <  1:735025/148‑1 (MQ=255)
                                                                                     aCGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAg                                                                >  2:2792685/1‑148 (MQ=255)
                                                                                     aCGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAg                                                                >  2:3709947/1‑148 (MQ=255)
                                                                                         cTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGAtg                                                            >  1:474313/1‑148 (MQ=255)
                                                                                         cTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGAtg                                                            >  1:2364307/1‑148 (MQ=255)
                                                                                         cTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGAtg                                                            >  1:946092/1‑148 (MQ=255)
                                                                                           gAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGAtgtg                                                          <  2:366742/148‑1 (MQ=255)
                                                                                            aaGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGc                                                         <  1:838212/148‑1 (MQ=255)
                                                                                               gttgttTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTg                                                                          <  2:631601/128‑1 (MQ=255)
                                                                                               gttgttTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGcaca                                                      >  1:781371/1‑148 (MQ=255)
                                                                                               gttgttTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGcaca                                                      <  1:1587958/148‑1 (MQ=255)
                                                                                                ttgttTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACaa                                                     <  1:2943510/148‑1 (MQ=255)
                                                                                                ttgttTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACaa                                                     >  1:454908/1‑148 (MQ=255)
                                                                                                   ttttGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGa                                                  <  1:199202/148‑1 (MQ=255)
                                                                                                       gAAATTTGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGAAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTc                                              <  2:3164773/148‑1 (MQ=255)
                                                                                                       gAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTc                                              <  1:2557203/148‑1 (MQ=255)
                                                                                                       gAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTc                                              <  1:2533742/148‑1 (MQ=255)
                                                                                                         aaTTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAg                                            <  2:147251/148‑1 (MQ=255)
                                                                                                         aaTTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAg                                            <  2:3048844/148‑1 (MQ=255)
                                                                                                          aTTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGc                                           >  1:420376/1‑148 (MQ=255)
                                                                                                           ttGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACACTGAATTCAGCa                                          >  2:3457624/1‑148 (MQ=255)
                                                                                                            tGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTGTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCAt                                         >  2:1978673/1‑148 (MQ=255)
                                                                                                                 tCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTtgtg                                    >  2:1799187/1‑148 (MQ=255)
                                                                                                                       tCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGtt                              >  1:1277382/1‑148 (MQ=255)
                                                                                                                          gTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTg                           <  2:2730888/148‑1 (MQ=255)
                                                                                                                          gTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTg                           <  1:3374692/148‑1 (MQ=255)
                                                                                                                          gTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTg                           <  1:1474524/148‑1 (MQ=255)
                                                                                                                             aCCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACa                        >  1:1105475/1‑148 (MQ=255)
                                                                                                                             aCCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACa                        >  2:2181054/1‑148 (MQ=255)
                                                                                                                             aCCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACa                        >  1:1320044/1‑148 (MQ=255)
                                                                                                                             aCCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACa                        >  1:3779724/1‑148 (MQ=255)
                                                                                                                              cccAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGCTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAg                       >  2:1847767/1‑148 (MQ=255)
                                                                                                                               ccAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGt                      <  2:555259/148‑1 (MQ=255)
                                                                                                                                 aGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGttt                    <  2:2546407/148‑1 (MQ=255)
                                                                                                                                   tATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGttttt                  <  2:2702279/148‑1 (MQ=255)
                                                                                                                                   tATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGttttt                  >  1:1410385/1‑148 (MQ=255)
                                                                                                                                    aTTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTc                 <  1:1037201/148‑1 (MQ=255)
                                                                                                                                     tttACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCt                <  1:1158942/148‑1 (MQ=255)
                                                                                                                                     tttACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAACTCAGCATTTGTGCTTGTTCTGACAGTTTTTCt                <  1:2798430/148‑1 (MQ=255)
                                                                                                                                     tttACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCt                <  1:2645709/148‑1 (MQ=255)
                                                                                                                                      ttACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCtt               >  2:2837077/1‑148 (MQ=255)
                                                                                                                                          tGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGttt           <  1:3209780/148‑1 (MQ=255)
                                                                                                                                                tGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGgag     >  1:2318586/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AACGTCACCTGGAGGCACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGC  >  CP009273/1643808‑1644102

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: