Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,222,656 A→C 23.8% 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,6560AC23.8% 92.7 / 2.0 72intergenic (+112/+261)yohP/dusCuncharacterized protein/tRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand):  ref base A (4/38);  new base C (4/9);  total (15/57)
Fisher's exact test for biased strand distribution p-value = 7.91e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCA  >  CP009273/2222509‑2222803
                                                                                                                                                     |                                                                                                                                                    
cTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCa                                                                                                                                                      <  2:1106038/148‑1 (MQ=255)
cTGGTGGTGACTGGCGTATCTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCa                                                                                                                                                      <  2:1103005/148‑1 (MQ=255)
  ggtggtGACTGGCGTATTTAAGATGATATTTTAAAATTAATTACTGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAgg                                                                                                                                                    <  2:505567/148‑1 (MQ=255)
   gtggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCATGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGAAggg                                                                                                                                                   >  2:2453205/1‑148 (MQ=255)
    tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAgggg                                                                                                                                                  <  1:2387852/148‑1 (MQ=255)
      gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggc                                                                                                                                                <  2:2085578/148‑1 (MQ=255)
      gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggc                                                                                                                                                <  2:2085430/148‑1 (MQ=255)
         aCTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggcggc                                                                                                                                            <  1:2923076/148‑1 (MQ=255)
         aCTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggcggc                                                                                                                                            <  1:134011/148‑1 (MQ=255)
         aCTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggcggc                                                                                                                                            <  1:1640749/148‑1 (MQ=255)
            ggCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCgggggggggccc                                                                                                                                          >  2:1944760/1‑144 (MQ=255)
            ggCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAAGCACTGGCGCAAGGAGCGCGgggggggggggccc                                                                                                                                          >  1:278953/1‑144 (MQ=255)
             gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAt                                                                                                                                        <  2:1530579/148‑1 (MQ=255)
             gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCGTCATATAAGCACTGGCGCAAGGAGCGCGGCCGGGGGGGGcccc                                                                                                                                         >  2:1905444/1‑146 (MQ=255)
             gCGTATTTAAGATGATATTTTAAAAATAATTAAAGTCATCAGGTCCGAAAATAACTAGAATATTTCAGTCTCTCATCCTGTTGAGATCCTGTCATGTGCATTGCTTCATATAATCACTGGCTCAAGTAGCGCGCAGGGGGCGGCCAAt                                                                                                                                        >  1:2596383/1‑148 (MQ=255)
              cGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGGGGGCCaaa                                                                                                                                        >  2:1941081/1‑147 (MQ=255)
              cGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGgggggggggggggcgg                                                                                                                                        >  1:630893/1‑142 (MQ=255)
              cGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGggggggggggggcgcc                                                                                                                                        >  2:1208981/1‑142 (MQ=255)
              cGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGggggggggggggccac                                                                                                                                        >  2:2664390/1‑142 (MQ=255)
              cGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCGGGGGGGGGcccagc                                                                                                                                       >  2:1941198/1‑144 (MQ=255)
                    tAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATcgccgcc                                                                                                                                 <  1:1093013/148‑1 (MQ=255)
                     aaGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATcgccgccg                                                                                                                                <  2:3147517/148‑1 (MQ=255)
                     aaGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATcgccgccg                                                                                                                                <  2:3147535/148‑1 (MQ=255)
                        atgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGGGGGCCAACGCgcccgccc                                                                                                                             >  2:2087802/1‑148 (MQ=255)
                        atgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGGGGGGGGGGGCGAAAGCCCCCGccc                                                                                                                             >  2:1386597/1‑148 (MQ=255)
                         tgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGcccc                                                                                                                            <  2:3470362/148‑1 (MQ=255)
                          gatATTTTAAAATTAATTAAGGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGccccc                                                                                                                           <  2:2596383/148‑1 (MQ=255)
                               tttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCAGGGGGCGGCCACACccccccccccccccac                                                                                                                       >  2:26081/1‑144 (MQ=38)
                                       taattaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGGGGGCCCAAGCCCCCCCCCCCCCCCCCCCCCCgc                                                                                                                >  2:158647/1‑147 (MQ=255)
                                       taattaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGGGGGGGGGGGCCCCCGcccccccccccccccccccccccc                                                                                                              >  1:1163202/1‑144 (MQ=255)
                                           taatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCGGGGGGGGGGCCAAAGCCGCCCCCCCCAGCCCCCCCCGGCGCTg                                                                                                          >  2:2303928/1‑148 (MQ=37)
                                           taatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGGGGGCCAAAGCCGCGCCCCCCAGCACCCCCCGGcgcag                                                                                                          >  2:555409/1‑144 (MQ=37)
                                                tcatcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGGGGCCCCCAGCCCCCGCCCCCCGCCCCCCCCCGCCCTGGCGaa                                                                                                     >  2:2727731/1‑148 (MQ=255)
                                                     aGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGAAAGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACaaaa                                                                                                <  2:1556884/148‑1 (MQ=255)
                                                        tCCTAAAATAACGTTAATTTTTCTGTCTCTCATCCTGTTGCGCGCGTGTCTTGTGCATTGCTTCATATAATCAGTGGCGCATGGCGCGCGCTGGGGGCGGCAAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATcg                                                                                             <  1:1801910/148‑1 (MQ=21)
                                                        tCCTAAAATAACGAGAATATTTCTGTCTCTCATCGTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACGGGCGCAAGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATcg                                                                                             <  1:3420224/148‑1 (MQ=21)
                                                          cGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATcgcc                                                                                           <  2:576107/148‑1 (MQ=255)
                                                                   cGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCTTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCgg                                                                                  <  2:2915663/148‑1 (MQ=255)
                                                                    gataaTATTTCTGTCTCTCATCGTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCACGGAGCGCGCCGGGGGCGGCCAATCCCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGt                                                                                 <  2:309765/145‑1 (MQ=18)
                                                                    gataaTATTTCTGTCTCTCATCGTGTTGCGCTCCTGACATGTGCTTTGCTTCATTTAATCCCTGGCGCAAGGAGGGCGCTGGGGGCGGCAACTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGt                                                                                 <  2:309787/145‑1 (MQ=18)
                                                                       aaTATTTCAGTCTCCCAGCCTGTTGCGCTCCTGTCATGTGCATTGCTACATATAAGCACTGGCGCACGGAGCGCGCAGGGGGCGGCAAATCCCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGcc                                                                              <  2:1538922/148‑1 (MQ=18)
                                                                         tttttCATTCTCTCATCGTTTTGCGCTCCTGTATTGTGATTTGCTTAATATAATCACGGGCCAATGGAGCGCGCCGGGGGCGCCAAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCt                                                                            <  1:3466005/146‑1 (MQ=17)
                                                                              cAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCACGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCt                                                                       <  1:2347551/148‑1 (MQ=255)
                                                                              cAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTCGCGCAAGGAGCGCGCCAGGGGCGGCCAAACGCCCCCGCCCCCAGCACCCCCCCCCACAGGCGAAAAAAAAAACCCCAGCGCGGGGCCCACGGcc                                                                       >  2:591159/1‑147 (MQ=255)
                                                                              cAGTCTCTCATCCTGTGGCGCTCCTGTCATGTGCTTTGCTTCATATAATCCCGGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCt                                                                       <  1:2347303/148‑1 (MQ=255)
                                                                                ggctctcATCCTGTTGCGCTCCTGTAATGTGATTTGCTTAATATAATCACTGGCGCAAGGAGCGCGCCGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGATGCGCGGTGCCCTCGGCTTa                                                                     <  1:2330540/146‑1 (MQ=17)
                                                                                 tctctcATCCTGTTGGGCTCCTGTCATGTGCATTGCTTCATATAATCAGTGGCGAAAGGAGCGCGCTGGGGGGGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTAt                                                                    <  2:1842605/148‑1 (MQ=18)
                                                                                 tctctcATCCTGTTGCGCTCCGGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTAt                                                                    <  2:468521/148‑1 (MQ=255)
                                                                                  ctctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATc                                                                   <  2:589199/148‑1 (MQ=255)
                                                                                  ctctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATc                                                                   <  2:590382/148‑1 (MQ=255)
                                                                                      cATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCACGGAGCGCGCCGGGGGCGGCAAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCtt                                                               <  1:1463150/148‑1 (MQ=255)
                                                                                       aTCGTGTTGCGCTCTTTTATTGTGATTTTTTTAATATAATCAGGGGCGAATGGCGGGCGCGGGGGGCGCCAAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTa                                                              <  1:1325395/148‑1 (MQ=2)
                                                                                           tGTTGCGCTCCTGTCTTTTGATTTTCTTAATATAACCCCTGGCGCACGGCGCGCGCTGGGGGCGCCCAACCCCCCCCCCCCCCGCCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGc                                                          <  1:141599/148‑1 (MQ=14)
                                                                                           tGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCACGGAGCGCGCCGGGGGCGGCCAAGCGCCGCCCCCCCGTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGc                                                          <  1:2553271/148‑1 (MQ=17)
                                                                                           tGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGc                                                          <  2:2630752/148‑1 (MQ=255)
                                                                                           tGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGc                                                          <  2:2627470/148‑1 (MQ=255)
                                                                                           tGTTGCGCTCCTGTAATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGc                                                          <  2:2630757/148‑1 (MQ=255)
                                                                                                  cTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCggg                                                   <  2:2356418/148‑1 (MQ=255)
                                                                                                     cTGTCATGTGCATTGCTTAATATAATCCCTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCg                                                <  2:2572359/148‑1 (MQ=255)
                                                                                                       gTCTTGTGCATTGATTCATATAATCCCTGGCGCAAGGAGCGCGCAGGGGGCGCCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCggg                                              <  2:2364256/148‑1 (MQ=255)
                                                                                                           tgtgCATTGATAAATATAATCACTGGCGCAAGGAGCGCGCCGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGcgcg                                          <  1:1386597/148‑1 (MQ=255)
                                                                                                                atggCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTa                                     <  2:699481/145‑1 (MQ=255)
                                                                                                                  tgataaatatAAGCCCGGGCGCACGGTGCGCGCCGGGGGCGCCCAAGCGCCGCCCCCCCCGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAAc                                   <  1:2990541/142‑1 (MQ=11)
                                                                                                                  tGCTTAATATAACCCGTGGCGCATGGAGCGCGCGGGGGGCGCCCAACCGCCGCCCCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAAc                                   <  1:2992462/148‑1 (MQ=11)
                                                                                                                    cTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACAt                                 <  2:430919/148‑1 (MQ=40)
                                                                                                                     ttaatatAACCAGTGGCGCACGGAGGGCGCGGGGGGCGCCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATc                                <  1:3376167/145‑1 (MQ=11)
                                                                                                                         tataATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTg                            <  2:395660/148‑1 (MQ=37)
                                                                                                                             aTCACTGGCGCATGGTGCGCGCTGGGGGCGGCCACTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTaaa                        <  1:821254/148‑1 (MQ=255)
                                                                                                                             aTCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTaaa                        <  2:1847453/148‑1 (MQ=33)
                                                                                                                                aCTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcg                     <  2:3496484/148‑1 (MQ=32)
                                                                                                                                aCTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcg                     <  2:3496929/148‑1 (MQ=32)
                                                                                                                                aCTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcg                     <  2:3500566/148‑1 (MQ=32)
                                                                                                                                aCTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcg                     <  2:3500980/148‑1 (MQ=32)
                                                                                                                                   ggCGCACGGAGGGCGCCGGGGGCGCCCAATCGCGGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                  <  2:2090793/148‑1 (MQ=11)
                                                                                                                                   ggCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                  <  2:321805/148‑1 (MQ=32)
                                                                                                                                      gaaaGGAGGGCGGCGGGGGGGGCAACTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTTCCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCt               <  1:232842/146‑1 (MQ=11)
                                                                                                                                                  gcAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATc    <  2:3509434/148‑1 (MQ=32)
                                                                                                                                                    cAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATcc   <  2:1163207/148‑1 (MQ=32)
                                                                                                                                                     aGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCa  <  2:1358952/148‑1 (MQ=32)
                                                                                                                                                     |                                                                                                                                                    
CTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCA  >  CP009273/2222509‑2222803

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: