Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,222,656 A→C 30.0% intergenic (+112/+261) yohP → / ← dusC uncharacterized protein/tRNA‑dihydrouridine synthase C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,222,6560AC30.0% 62.2 / 5.4 69intergenic (+112/+261)yohP/dusCuncharacterized protein/tRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand):  ref base A (3/28);  new base C (4/9);  total (13/56)
Fisher's exact test for biased strand distribution p-value = 1.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCC  >  CP009273/2222509‑2222802
                                                                                                                                                    |                                                                                                                                                    
cTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCa                                                                                                                                                      <  2:1645275/148‑1 (MQ=255)
    tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAgggg                                                                                                                                                  <  1:3072047/148‑1 (MQ=255)
     ggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAggggg                                                                                                                                                 <  1:171843/148‑1 (MQ=255)
        gACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggcgg                                                                                                                                              <  1:2127139/148‑1 (MQ=255)
         aCTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggcggc                                                                                                                                             <  2:796675/148‑1 (MQ=255)
           tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCa                                                                                                                                           <  1:1774754/148‑1 (MQ=255)
            ggCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCaa                                                                                                                                          <  1:143122/148‑1 (MQ=255)
            ggCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCaa                                                                                                                                          <  1:2558651/148‑1 (MQ=255)
             gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGgggggggggggggcg                                                                                                                                         >  1:376907/1‑143 (MQ=255)
             gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGgggggggggggggca                                                                                                                                         >  1:2173661/1‑143 (MQ=255)
             gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGggggggggggggccg                                                                                                                                         >  1:1305148/1‑143 (MQ=255)
             gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGgggggggggggcccg                                                                                                                                         >  2:1509586/1‑143 (MQ=255)
               gTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATcg                                                                                                                                      <  2:789400/148‑1 (MQ=255)
                tATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGGGGGCCaacgc                                                                                                                                       >  2:3085761/1‑145 (MQ=255)
                tATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGGCGGGGGGGccccaccc                                                                                                                                      >  2:3085999/1‑142 (MQ=255)
                tATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATcgc                                                                                                                                     <  1:638204/148‑1 (MQ=255)
                 aTTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGGGGGcccacgcc                                                                                                                                      >  2:1076352/1‑142 (MQ=255)
                   ttAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCa                                                                                                                                           <  1:887022/140‑1 (MQ=255)
                   ttAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCa                                                                                                                                           >  2:887022/1‑140 (MQ=255)
                         tgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGcccc                                                                                                                            <  1:3156169/148‑1 (MQ=255)
                           atatTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCTGGGGGGGGCCAGAGGCGCCGcccccc                                                                                                                          >  2:2616421/1‑147 (MQ=255)
                              ttttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGCGGGCCAACGCCCCCCCCCCCCGcc                                                                                                                       >  2:3510468/1‑147 (MQ=255)
                              ttttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCa                                                                                                                       <  2:96169/148‑1 (MQ=255)
                              ttttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCa                                                                                                                       <  2:1507388/148‑1 (MQ=255)
                                    aattaattAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGGGGGGGGGGGCGCCCCCCGCccccccccccccccccc                                                                                                                 >  1:769841/1‑147 (MQ=255)
                                       taattcatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGGGGGCCAACCCCCCCCCCCCCCGCAACCCCCGGc                                                                                                              >  2:586804/1‑148 (MQ=37)
                                       taattaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCGTGTTGCGCTCCTGTCATGTGCATTGCTTAATATAATCAGTGGCGAATGGCGCGCGCTGGGGGGGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGc                                                                                                              <  1:1762710/148‑1 (MQ=37)
                                       taattaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCCGGGGGGGCCCAAAACCCGCCCCCCCCCCAcccccccgc                                                                                                              >  2:586936/1‑144 (MQ=255)
                                       taattaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGAGGGGGGGGGCCAACGCCGCCCCCCCCAGCACCCCCCGGc                                                                                                              >  2:2377490/1‑148 (MQ=37)
                                        aattaatGTCATCAGGTCAGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCTTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGct                                                                                                             <  2:447818/148‑1 (MQ=255)
                                          ttaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCCGGGGGGGCCCACCGCCCCCCCCCCCCCCCCCCCCCGGccca                                                                                                           >  2:3194148/1‑145 (MQ=255)
                                           taatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGGGGGCCCAACCCGCCCCCCCCCCCCCCCCCCGGcccgg                                                                                                          >  2:1181408/1‑144 (MQ=255)
                                                           gAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGGGGCCCAAAGCCCCCCCCCCCGGCACCCCCCGGCGCTGGCCGAACaaaaacccg                                                                                          >  2:2269311/1‑143 (MQ=25)
                                                            aaaaTAACGAGAATTTTTCTGTCTCTCATCGTGTTGCGGTCGTGTCTTGTGCTTTGCTTCATATAATCCGGGGCGCATGGAGCGCGCTGGGGGCGGCAAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCGGGCGAACAAAATcgccgc                                                                                        <  2:2438892/148‑1 (MQ=18)
                                                                tAACGATAATATTTCTGTCTCTCAGCGTGTTGCGCTCCGTTCTTGTGCATTGCTTCATATAATCACGGGCGCATGGAGCGCGCTGGGGGCGGCAAATCGCCGCCGCCCCCGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTgcg                                                                                    <  1:2943635/148‑1 (MQ=18)
                                                                   ctagaTTATTTCAGTCTCTCATCGTTTTGCGCTCCTGTCATGTGAATTGCTTCATATAATCACGGGCGCATGGAGCGCGCTGGGGGCGGCCAACCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCgg                                                                                 <  1:550048/146‑1 (MQ=18)
                                                                        ttatTTCAGTCTCTCATCGTGTTGCGCTCGTGACTTGTGCATTGCTACAAATACTCACTGGCGCACGGAGCGCGCCGGGGGCGGCAAAGCGCCGCCCCCCCGGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGccc                                                                            <  2:1231759/147‑1 (MQ=17)
                                                                        atttttCGGTCTCTCATCCGGTGGCGCCCCTGTCTTGTGATTGGCTTCATATAACCCCGGGCGCCGGGCGCGCGCGGGGGGCGCCCACGCGCCGCCGCCCCCTCCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGccc                                                                            <  1:2880536/145‑1 (MQ=12)
                                                                             tCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGc                                                                       <  2:2498367/148‑1 (MQ=255)
                                                                                 tctctcTTCCTGTTGCGCTCCTGTCTTGTGCTTTGCTTAATATAATCACGGGCGCACGGCGCGCGCGGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTAt                                                                   <  1:2297444/148‑1 (MQ=18)
                                                                                 tctctcATCCTGTGGCGCTCCTGTAATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTAt                                                                   <  2:2157372/148‑1 (MQ=255)
                                                                                          ttGTTGCGCTCCTGTCTTTTGATTTTATTAATATAAACCCTGGCGAATGGAGCGCGCTGGGGGGGGCAAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACgg                                                          <  1:1098902/147‑1 (MQ=14)
                                                                                          gtGTTGCGCTCGTGTATGGTGCTTTGCTTCATATACTCCGTGGCGCATGGAGCGCGCGGGGGGCGGCAAATCGCCGCCGCCCCCGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACgg                                                          <  1:1101035/147‑1 (MQ=14)
                                                                                          cTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACGGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACgg                                                          <  2:2976614/148‑1 (MQ=255)
                                                                                                   tCCTTTCATGTGATTTGTTTAATATAATCACGGGCCCATGGCGGGCGCCGGGGGCGCCAATTCGCCGCCGCCCCCGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGt                                                 <  1:1089472/148‑1 (MQ=14)
                                                                                                   tCCTGTCATGTGCATTGCTTCATATAATCCCGGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGt                                                 <  2:1127335/148‑1 (MQ=255)
                                                                                                   gccTGTCTTGTGATTTGCTACATATAATCCCGGGCGCACGGAGCGCGCGGGGGGCGCCCAAGCGCCGCCGCCCCGGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGt                                                 <  1:3126102/147‑1 (MQ=14)
                                                                                                          aTGTGATTTGCTTCATATAATCACTGGCGCAAGGCGCGCGCCGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGgcgc                                          <  1:2061552/148‑1 (MQ=255)
                                                                                                            gtgCTTTGCTTCATATACCCCCGGGCGCACGGCGCGGGCGGGGGGCGCCCAAGCCCCCCCCCCCCCGCCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGa                                        <  1:344157/148‑1 (MQ=11)
                                                                                                             tgattggCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAg                                       <  2:913748/142‑1 (MQ=255)
                                                                                                             tgCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAg                                       <  2:1752164/148‑1 (MQ=255)
                                                                                                             tgCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATAGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAg                                       <  2:208604/148‑1 (MQ=255)
                                                                                                             tgCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATAGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAg                                       <  2:208479/148‑1 (MQ=255)
                                                                                                             tgCATTGCTTCATATAAGCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATAGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAg                                       <  2:1575436/148‑1 (MQ=255)
                                                                                                              gctttGTTTCATATAACCACGGGCGCAAGGCGCGCGCGGGGGGCGGCCAAGCGCCGCCGCCCCCGGCACCCCCGGGCTCGGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAgg                                      <  1:2539247/145‑1 (MQ=11)
                                                                                                               cATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGt                                     <  2:3033949/148‑1 (MQ=255)
                                                                                                               cATTGCTACATATAATCCCGGGCGCACGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGt                                     <  1:573539/148‑1 (MQ=14)
                                                                                                                 ttGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTaa                                   <  2:376907/148‑1 (MQ=255)
                                                                                                                 ttGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTaa                                   <  2:2173661/148‑1 (MQ=255)
                                                                                                                 ttGCTTAATAAAACCACGGGCGCACGGCGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTCCACCCCCGGGCTCGGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTaa                                   <  1:2498014/148‑1 (MQ=11)
                                                                                                                       cATATACTCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATccc                             <  2:3197813/148‑1 (MQ=25)
                                                                                                                             atccggggCGCAGGGCGCGGGCGGGGGGGGCCAACGCGCCGCGGCCCCGTCCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTaaa                       <  1:2984798/142‑1 (MQ=11)
                                                                                                                                tcTGGCGCAGGGTGGGCGCGGGGGGCGGCCACGCGCCGCCGCCCCCGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcg                    <  1:934588/147‑1 (MQ=11)
                                                                                                                                  gggCGCATGGCGCGCGCGGGGGGGGCCAAAGGGCCGCCCCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgcg                  <  1:2992508/147‑1 (MQ=11)
                                                                                                                                    gcgcAAGGAGCGCGCCGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGcc                <  2:3068007/148‑1 (MQ=17)
                                                                                                                                     cgcATGGAGCGCGCAGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGccc               <  2:1231974/148‑1 (MQ=17)
                                                                                                                                     cgcAAGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGccc               <  2:1231828/148‑1 (MQ=17)
                                                                                                                                      gcAAGGAGCGCGCCGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCt              <  2:1559211/148‑1 (MQ=17)
                                                                                                                                      gcAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCt              <  2:1559105/148‑1 (MQ=32)
                                                                                                                                         aGGCGCGCGCCGGGGGCGCCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATAGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAg           <  1:325848/148‑1 (MQ=11)
                                                                                                                                         aGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAg           <  1:1477968/148‑1 (MQ=17)
                                                                                                                                          ggAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGc          <  1:3526890/148‑1 (MQ=17)
                                                                                                                                                cgcCGGGGGGGGCCAAGCGCCGCCGCCCCCGCCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACAt    <  1:1277387/148‑1 (MQ=11)
                                                                                                                                                   cAGGGGGCGGCCACTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATcc  <  2:68837/148‑1 (MQ=17)
                                                                                                                                                    |                                                                                                                                                    
CTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCC  >  CP009273/2222509‑2222802

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