Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913646,7470TG42.4% 145.6 / 25.5 85T255P (ACC→CCC) citG2‑(5''‑triphosphoribosyl)‑3'‑dephosphocoenzyme‑A synthase
Reads supporting (aligned to +/- strand):  ref base T (5/44);  new base G (36/0);  total (41/44)
Fisher's exact test for biased strand distribution p-value = 6.02e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCC  >  NC_000913/646599‑646897
                                                                                                                                                     |                                                                                                                                                         
aaTATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGtt                                                                                                                                                          <  1:360864/150‑1 (MQ=255)
   atTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑Ga                                                                                                                                                      <  1:140158/150‑1 (MQ=255)
     ttCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑t                                                                                                                                                   <  2:1264950/150‑1 (MQ=255)
     ttCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑t                                                                                                                                                   <  2:618633/150‑1 (MQ=255)
      tCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑Tg                                                                                                                                                  <  2:1242707/150‑1 (MQ=255)
       ccAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCACGTAAGGATCAGTAGGCCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTc                                                                                                                                                         >  2:933682/1‑144 (MQ=255)
           aTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCt                                                                                                                                           <  2:1348080/151‑1 (MQ=255)
              tATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCttt                                                                                                                                         <  2:87924/150‑1 (MQ=255)
                tCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑GGA‑TG‑CCCCCcttttt                                                                                                                                      >  2:918959/1‑150 (MQ=255)
                   aGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCTTTTTGCaa                                                                                                                                   <  1:949868/151‑1 (MQ=255)
                    gTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAA‑TG‑CCCCCctttt                                                                                                                                       >  1:405855/1‑146 (MQ=255)
                     tGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTG‑GAA‑GG‑CCCCCttttt                                                                                                                                       >  2:1062351/1‑146 (MQ=255)
                     tGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GGA‑TG‑CCCCCCTTTTGTaat                                                                                                                                  >  1:1164595/1‑150 (MQ=255)
                      gCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCGCTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCTTTTTGCaata                                                                                                                                 <  2:309003/150‑1 (MQ=255)
                      gCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑GAA‑TGTCCCCCttttttaata                                                                                                                                  >  1:767722/1‑145 (MQ=255)
                         tAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAA‑TT‑CCCCCcttttttcatta                                                                                                                                >  1:1439006/1‑141 (MQ=255)
                          aaaTAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCTTTTTTCAATAATGtt                                                                                                                            >  2:1623948/1‑151 (MQ=255)
                          aaaTAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCTTTTTGCAATAATGt                                                                                                                             <  1:426189/150‑1 (MQ=255)
                          aaaTAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑GAC‑TG‑ACCCCCTTTTTCAAATATGt                                                                                                                             >  1:1016814/1‑150 (MQ=255)
                           aataatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCTTTTTGCAATAATGttt                                                                                                                           <  2:1003619/151‑1 (MQ=255)
                            ataatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GGA‑TG‑CCCCCCTTTTTTAATAATtttt                                                                                                                           >  1:769738/1‑150 (MQ=255)
                              aatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCTTTTTGCAATAATGTTTgcg                                                                                                                        <  1:1340018/151‑1 (MQ=255)
                              aatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTACCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GGA‑TG‑CCCCCCTTTTTCAATAAAGttt                                                                                                                           >  1:277956/1‑148 (MQ=255)
                               atTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCTTTTTGCAATAATGTTTgcgc                                                                                                                       <  1:567998/151‑1 (MQ=255)
                                   aaTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCTTTTTGCAATAATGTTTGCGCCTcg                                                                                                                   <  1:1378214/151‑1 (MQ=255)
                                      cTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑GAA‑AG‑CCCCCCTTTTTTAATAATGtt                                                                                                                            >  1:1615420/1‑139 (MQ=255)
                                      cTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAA‑TG‑GCCCCCTTTTTCAATAATTTTTTCGcccc                                                                                                                    >  2:1260302/1‑145 (MQ=255)
                                       tgtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑GAA‑TG‑CCCCCCTTTTTCAATAAAGTTTTCGcccccc                                                                                                                  >  2:1062887/1‑144 (MQ=255)
                                        gtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGCGATAATCGAGATCGGC‑GGGGGGGC‑AAA‑AG‑ACCCCCTTTTGTTATTATGTTTActcctctccc                                                                                                                >  1:614680/1‑145 (MQ=255)
                                        gtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCTTTTTGCAATAATGTTTGCGCCTcgc                                                                                                                  <  2:1086048/147‑1 (MQ=255)
                                        gtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑GAA‑TG‑CCCCCCTTTTTCAATAATTTTTTCGCCCcgcg                                                                                                                 >  2:1070934/1‑148 (MQ=255)
                                            tAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAA‑TG‑CCCCCCTTTTTCAATAATGTTTTCGCCCCCCGCCGTAGc                                                                                                          >  2:1151214/1‑151 (MQ=255)
                                             aaaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGcc                                                                                                         <  1:862215/151‑1 (MQ=255)
                                             aaaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGc                                                                                                          <  2:262890/150‑1 (MQ=255)
                                              aaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCa                                                                                                        <  2:1537718/151‑1 (MQ=255)
                                              aaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCa                                                                                                        <  1:862000/151‑1 (MQ=255)
                                              aaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAA‑TG‑CCCCCTTTTTTCAATAATGTTTTCGCCCCGCGCTGTAGcc                                                                                                         >  2:532058/1‑150 (MQ=255)
                                              aaaaCCAGGTAAGGATCAGAAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑CG‑CCCCCTTTTTCCAATAATGTTTGCACCTCGCTCTTTAGcta                                                                                                        <  1:686037/151‑3 (MQ=255)
                                                    aGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAg                                                                                                  <  1:421749/151‑1 (MQ=255)
                                                       tAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑GAA‑TG‑CCCCCCTTTTTTCATAATGTTTTCGCCTCCCGCCGTAGCCAGCGCCGGcc                                                                                               >  2:389947/1‑151 (MQ=255)
                                                           gATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGGAA‑GG‑CCCCCTTTTTTC‑ATAATATTTTTGCCCCCCGCCGTTGCCCGCGCCGGgccccc                                                                                           >  1:990445/1‑151 (MQ=255)
                                                             tCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCC‑GAA‑GG‑CCCCCTTTTT‑TAATAATGTTTTGGCCCCGCCCCGGAGCCCGCGCCGGgccccccc                                                                                         >  1:995952/1‑149 (MQ=255)
                                                              cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTc                                                                                         <  2:1545492/150‑1 (MQ=255)
                                                              cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTc                                                                                         <  2:1416962/150‑1 (MQ=255)
                                                               aGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑GCA‑TG‑CCCCCCTTTTTCAATAATGTTTTTGCCCCCCGCCGTAGCCAGCGCCGGcccccccccc                                                                                       >  1:1458371/1‑147 (MQ=255)
                                                                  aGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCAc                                                                                    <  2:1662692/151‑1 (MQ=255)
                                                                  aGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAA‑TG‑CCCCCCTTTTTTAATAATGTTTTCGCCCCCCGCTGTAGCCAGCGCAGGccccccccccccc                                                                                    >  2:610012/1‑144 (MQ=255)
                                                                    gTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACcg                                                                                  <  1:222418/151‑1 (MQ=255)
                                                                     tCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAA‑TG‑CCCCCCTTTTTCAATAATTTTTTCGCCCCCCGCTGTAGCCCGCGCCGGgccccccccccccccc                                                                                 >  2:1571852/1‑141 (MQ=255)
                                                                      cAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACcgcg                                                                                <  2:31960/151‑1 (MQ=255)
                                                                      cAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACcgc                                                                                 <  1:1015415/150‑1 (MQ=255)
                                                                      cAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAA‑TG‑CCCCCCTTTTTCAATAATGTTTTTGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCCCGCcccccc                                                                                 >  2:1180936/1‑148 (MQ=255)
                                                                        gCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑GCA‑TG‑TCCCCCTTTTTTAATAATATTTTCGCCCCCCGCTGTAGCCAGCGCCGGcccccccccccccccc                                                                                 >  2:504646/1‑138 (MQ=255)
                                                                          aCTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GGG‑AG‑GCCCCTTTTTTCAATAATGTTTTCGCCTCACGCGGGAGGCAACGCGGGCCCCCCCCGccgcaccc                                                                                >  2:467274/1‑141 (MQ=255)
                                                                            tGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAA‑TG‑CCCCCTTTTTTCAATAATGTTTGCGCCTCCCCCTGTAGCCAGCGCCGGCCCCCCCCCCCACCCCGCGCTg                                                                           >  2:156476/1‑150 (MQ=255)
                                                                               cgccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GCA‑TG‑CCCCCTCTTTGCAATAAATTTTTTTCCTCCCGCCGGAGCCAGCGGCGGCCCCCCCCCCCAcccccccacgcca                                                                        >  1:1643478/1‑139 (MQ=255)
                                                                               cgccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAA‑TT‑CCCCCCTTTTTCAATAATGTTTGCGGCTCGCCCCGTAGCCCGCCCAGGCCCCCCCCCCCCCcgcccccagcca                                                                        >  2:613332/1‑141 (MQ=255)
                                                                                    ggACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTgg                                                                  <  1:211967/151‑1 (MQ=255)
                                                                                    ggACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTgg                                                                  <  2:1451747/151‑1 (MQ=255)
                                                                                       cTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GGA‑AG‑CCCCCCTTTTTCAATAATATTTTGGGCCCCCGCTGGAGGCAGCGCAGGGCccccccgcccccccgc                                                                               >  1:517626/1‑135 (MQ=255)
                                                                                            aTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCg                                                          <  2:292180/151‑1 (MQ=255)
                                                                                               tcgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGtt                                                        <  2:820265/150‑1 (MQ=255)
                                                                                               tcgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑GAA‑TG‑CCCCCCTTTTTCAAAAATGTTTTCGCCTCCCGCCGTAGCCAGCGCAGGCCCCCCCCCCCACCCCCCCCAGCAACGTTGGGATCGcccctg                                                       >  2:885826/1‑147 (MQ=255)
                                                                                                cgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTg                                                       <  1:1118528/150‑1 (MQ=255)
                                                                                                 gttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGAt                                                     <  2:294550/151‑1 (MQ=255)
                                                                                                 gttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAATTC‑CCCCCTTTTTTCAAAATTATTTTC‑CCCCGCCGTGCAGCCCGAGCAGGCCCCCCCCCCCCCCGCGCCCTGCAAcactgtgatccccg                                                          >  1:682106/1‑133 (MQ=255)
                                                                                                       tACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑GAA‑AG‑CCCCCCTTTTTCAATAATGTTTTCGCCCCCCCCTGTAGCCAGCGCCGGCCCCCCCCCCCCCCCCGCGATGCAACGCTGGGATCGCCCTTGATCGcc                                                 >  1:606499/1‑149 (MQ=255)
                                                                                                        acacTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCAt                                               <  1:1184303/150‑1 (MQ=255)
                                                                                                        acacTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑GGA‑TG‑CCCCCTTTTTTCCATAATATTTTTGCCCCCCCCTGTAGCCAGCGCAGGGCCCCCCCCCCACCCCCCGGTGCAACGTTGGGGTTGCCGCTGATCGCCCTc                                              >  1:1646288/1‑151 (MQ=255)
                                                                                                         cacTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCa                                             <  1:1685561/151‑1 (MQ=255)
                                                                                                          acTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATc                                              <  1:798121/149‑1 (MQ=255)
                                                                                                                 gTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAA‑TG‑CCCCCCTTTTTCAATAATGTTTTTGCCCCGCGCTGTAGCCAGCGCAGGCCCCCCCCCCCCCCCCGCGGAGCAACGTTGGGATCGCCGCTGATTGCCCTCAGTa                                          >  2:52997/1‑146 (MQ=255)
                                                                                                                     aaCTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAgg                                 <  1:871512/151‑1 (MQ=255)
                                                                                                                     aaCTGCCGGAGATAATCGAGATCGGC‑GGGGGGTG‑GAA‑AG‑GCCCCTTTTTTCAAAAAAGTTTTTGCCCCGCCCCGTAGGAAAAGGAGGGCCCCCCCGCCAACCGGGGAGGGAAAGGTGGGAAAGGCGGTGATAGCCATCAGAAGGAGAAAgg                                 >  2:713592/1‑151 (MQ=255)
                                                                                                                      aCTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCTTTTTGCAATAATGTGTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGGCATCAGTAGGAGCAAGGt                                <  1:1613522/151‑1 (MQ=255)
                                                                                                                         gCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATc                             <  1:838983/151‑1 (MQ=255)
                                                                                                                         gCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAA‑TG‑CCCCCCTTTTTCAATAAAGTTTGTGCCCCCCGCCGTTGCCCGGGCCGGCCCCCCCCGCCCCCCCGCGCCGCGACGTTTGTATCGGCCTTGGTCGCCCTCCGTAGGAGCCAggggac                             >  1:25692/1‑147 (MQ=255)
                                                                                                                           cGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGa                           <  1:1595313/151‑1 (MQ=255)
                                                                                                                           cGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGa                           <  2:780303/151‑1 (MQ=255)
                                                                                                                            ggAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCg                            <  2:1696811/149‑1 (MQ=255)
                                                                                                                              agaTAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCa                        <  2:707174/151‑1 (MQ=255)
                                                                                                                                   aTCGAGATCGGC‑GGGGGGTC‑GAA‑TG‑CCCCCTTTTTTTAATAAATTTTTTGCCTCGCGCTGTAGCCAGCGCAGGGCCCCCCCCCCACCCCGCGCTGCAAAGGTGGGTTCGCCCGTGGTCGCCCCCCGGAGGAGGAAGGGATCggggca                       >  1:209595/1‑142 (MQ=255)
                                                                                                                                             gc‑ggGGGTTC‑GAA‑TG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGAt         <  2:405870/151‑1 (MQ=255)
                                                                                                                                             gc‑ggGGGTTC‑GAA‑TG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGAt         <  2:1084255/151‑1 (MQ=255)
                                                                                                                                                 ggggTTC‑GAA‑TG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCt       >  2:631156/1‑150 (MQ=255)
                                                                                                                                                    gTTC‑GAA‑TG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAAcc   >  2:721389/1‑151 (MQ=255)
                                                                                                                                                     ttC‑GAA‑TG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAAccc  >  2:843892/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                         
AATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAA‑TG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCC  >  NC_000913/646599‑646897

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