Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913646,7470TG43.2% 136.9 / 12.1 81T255P (ACC→CCC) citG2‑(5''‑triphosphoribosyl)‑3'‑dephosphocoenzyme‑A synthase
Reads supporting (aligned to +/- strand):  ref base T (9/37);  new base G (35/0);  total (44/37)
Fisher's exact test for biased strand distribution p-value = 7.57e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGA‑A‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCT‑CGCCACCGCGCGATGCAACGTTGGTATCGCC‑GTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCC  >  NC_000913/646601‑646897
                                                                                                                                                  |                                                                                                                                                          
tatTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCgc                                                                                                                                                        >  1:630271/1‑150 (MQ=255)
 atTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGCGGtttgg                                                                                                                                                        >  1:348648/1‑147 (MQ=255)
   ttCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGG‑A‑Tg                                                                                                                                                   >  1:229529/1‑151 (MQ=255)
     ccAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGA‑A‑TGcc                                                                                                                                                 <  1:189435/151‑1 (MQ=255)
           ttATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGA‑A‑TGCCCCCctt                                                                                                                                           >  2:1258283/1‑151 (MQ=255)
            tATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGA‑A‑TGCCCCCttt                                                                                                                                           <  2:498584/150‑1 (MQ=255)
             aTCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGA‑A‑TGCCCCCttttt                                                                                                                                         <  1:1277594/151‑1 (MQ=255)
              tCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGACA‑GGCCCCCttttg                                                                                                                                         >  2:1126608/1‑150 (MQ=255)
               cAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGA‑A‑TGCCCCCTTTTTgn                                                                                                                                       <  2:269016/151‑2 (MQ=255)
               cAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGA‑A‑TGCCCCCTTTTTGc                                                                                                                                       <  1:131849/151‑1 (MQ=255)
                aaGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGA‑A‑TGCCCCCctttttc                                                                                                                                       >  1:1160777/1‑148 (MQ=255)
                  gTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGA‑A‑TGCCCCCCTTTTTCaat                                                                                                                                    >  1:925657/1‑151 (MQ=255)
                    gCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGA‑A‑TGGCCCCCTTTTTCaata                                                                                                                                   >  1:429443/1‑150 (MQ=255)
                       tAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTTGA‑A‑TGCCCCCTTTTTTCAATAATGt                                                                                                                               >  1:1178766/1‑151 (MQ=255)
                        aaaTAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGA‑A‑TGCCCCCCTTTTTCAATAATGtt                                                                                                                              >  1:1236987/1‑151 (MQ=255)
                          ataatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGA‑A‑TGCCCCCTTTTTGCAATAATGTTTg                                                                                                                            <  1:431347/151‑1 (MQ=255)
                           taatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGA‑A‑ATCCCCCTTTTTGCAATAATGTTTgg                                                                                                                           >  2:1458514/1‑150 (MQ=255)
                               tAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGA‑A‑TGCCCCCCTTTTTCAATAATGTTTGCGCCt                                                                                                                       >  1:622511/1‑151 (MQ=255)
                                   tCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGA‑A‑AGCCCCCCTTTTTCAATAATGTTTTCGccccccg                                                                                                                   >  1:1466598/1‑146 (MQ=255)
                                    cTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGA‑A‑TGCCCCCTTTTTGCAATAATTTTTTCGCCCcgccc                                                                                                                  >  1:1021207/1‑149 (MQ=255)
                                      gtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGA‑A‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTcgc                                                                                                                    <  2:1387014/147‑1 (MQ=255)
                                       tgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGA‑A‑TGCCCCCCTTTTTCAATAAAGTTTGCGCCCCGCGCCGt                                                                                                               >  1:3990/1‑151 (MQ=255)
                                        gCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGA‑A‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGt                                                                                                               <  1:680168/150‑1 (MQ=255)
                                        gCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGA‑A‑TGCCCCCCTTTTTCAATAATGTTTTCGCCCCGCGCTGt                                                                                                               >  1:1184864/1‑150 (MQ=255)
                                             aaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGA‑A‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAg                                                                                                         <  2:1092539/151‑1 (MQ=255)
                                             aaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGA‑A‑TGGCCCCTTTTTTGCATAATGTTTATGCCTCGCCCGGGAGCCac                                                                                                         >  1:1422584/1‑150 (MQ=255)
                                              aaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGA‑A‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAgc                                                                                                        <  1:123177/151‑1 (MQ=255)
                                              aaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGA‑A‑AGCCCCCCTTTTTCAATAATGTTTTTGCCTCCCGCTGGAGCCAgc                                                                                                        >  2:924113/1‑151 (MQ=255)
                                                ccAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGA‑A‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAgcgc                                                                                                      <  1:1372986/151‑1 (MQ=255)
                                                ccAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGA‑A‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAgcgc                                                                                                      <  1:1130980/151‑1 (MQ=255)
                                                  aGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGA‑ATTGCCCCCTTTTT‑CAATAATGTTTTCGCCCCCCCCCGTATCCAacgc                                                                                                      >  2:254426/1‑149 (MQ=255)
                                                  aGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGA‑A‑TGCCCCCTTTTTTCAATAATGTTTTTGCCTCCCGCCGTAGCCAgcgccg                                                                                                    >  2:590825/1‑149 (MQ=255)
                                                   ggTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGA‑A‑TGCCCCCCTTTTTCAATAAAGTTTTCGCCTCGCGCTGTAGCCAGCGCAgg                                                                                                   >  2:568420/1‑151 (MQ=255)
                                                      aaGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGA‑A‑TGGCCCCCTTTTTCCATAATGTTTTCGCCCCGCGCTGTAGCCAGCGCAGGccc                                                                                                >  2:1037916/1‑151 (MQ=255)
                                                      aaGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGCGGGTGGA‑A‑TGCCCCCTTTTTTCAAAAATGGTTTCGCCTCGCCCCGTAGCCAGCCCAGGgcc                                                                                                >  2:1245934/1‑151 (MQ=255)
                                                         gATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGA‑A‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGcccccc                                                                                             <  2:634074/151‑1 (MQ=255)
                                                          aTCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTGGA‑A‑TGCCCCCCTTTTTCAATAATTTATTCGCCCCCCGCTGGAGGCAGCCCAGGccccccc                                                                                            >  2:559300/1‑150 (MQ=255)
                                                          aTCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGA‑A‑TGCCCCCCTTTTTCAATAATGTTTGCGCCCCGCGCCGTAGCCCGCGCAGGgcccccc                                                                                            >  2:398521/1‑150 (MQ=255)
                                                            cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGA‑A‑TGCCCCCCTTTTGCAATAATGTTTTCGCCTCCCGCCGTAGCCCGCGCAGGccccccccc                                                                                          >  1:1024278/1‑148 (MQ=255)
                                                                 ggTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGA‑A‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCT‑CGCCAcc                                                                                    <  1:408842/151‑1 (MQ=255)
                                                                    cAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGA‑A‑TGCCCCCTTTTTGCAATAAATTTTTCGCCCCGCGCCGTAACCCGCGCACGccccccccccccccccc                                                                                  >  2:824377/1‑140 (MQ=255)
                                                                        aCTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGA‑A‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCT‑CGCCACCGCGCGAt                                                                             <  2:885802/151‑1 (MQ=255)
                                                                        aCTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCCA‑C‑CGCCCCCTCTTTTCCATCATATTTTTGCGCCGCGACGCAGCCCACGCCGGCCCCCCT‑CGCCTCcgcccgc                                                                              >  2:1317464/1‑149 (MQ=255)
                                                                           gccgccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGA‑A‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCT‑CGCCACCGCGCGATGCa                                                                          <  2:1403844/151‑1 (MQ=255)
                                                                             cgccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGA‑A‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCT‑CGCCACCGCGCGATGCAAc                                                                        <  2:721060/151‑1 (MQ=255)
                                                                                 tGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGA‑A‑TGCCCCCCTTTTTCAATAATGTTTTCGCCTCGCGCCGTAGCCAGCGCCGGCCCCCCC‑CCCCCCCCCGCGCTGCAACgct                                                                     >  2:339555/1‑148 (MQ=255)
                                                                                  ggACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGA‑A‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCT‑CGCCACCGCGCGATGCAACGTTg                                                                    <  1:981360/150‑1 (MQ=255)
                                                                                       gagaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGA‑A‑AGGCCCCTTTTTTCAATAATGTTTTCGCCCCGCGCCGTAGCCAGCGCAGGCCCCCCC‑CCCCACCCCCCCATGGAACGATGGGATCg                                                              >  1:1159570/1‑151 (MQ=255)
                                                                                         gaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGA‑A‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCT‑CGCCACCGCGCGATGCAACGTTGGTATCGcc                                                            <  2:1352371/151‑1 (MQ=255)
                                                                                          aTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGA‑A‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCT‑CGCCACCGCGCGATGCAACGTTGGTATCGCC‑g                                                          <  1:1125381/151‑1 (MQ=255)
                                                                                          aTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGA‑A‑TGCCCCCCTTTTTCAATAATGTTTTCGCCTCCCGCTGTTGCCCGCGCCGGCCCCCCC‑CCCCACCCCGCGCTGCAACGGTGGGATCGCC‑g                                                          >  2:341357/1‑151 (MQ=255)
                                                                                           tttCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGA‑A‑TGGCCCCTTTTTTCCATAATTTTTTCGCCTCGCGCTGCTGCCAGCGCAGGCCCCCCC‑CGCCCCCGCGCGGTGGCACGGTGGGATccgcgc                                                          >  1:1172630/1‑146 (MQ=255)
                                                                                             tcgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCTC‑A‑TTCCCCCCTTTTTTTATAATATTTGGGCCTCCCGCTGGAGCCAGCGCCGGGCCCCCC‑CCCCCCCCCCCGGTGCAAAGTTTGGATCGcccgtg                                                        >  2:943157/1‑147 (MQ=255)
                                                                                              cgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGA‑A‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCT‑CGCCACCGCGCGATGCAACGTTGGTATCGCC‑GTTGa                                                      <  2:229533/151‑1 (MQ=255)
                                                                                                 tcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGA‑A‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCT‑CGCCACCGCGCGATGCAACGTTGGTATCGCC‑GTTGATc                                                    <  2:203958/150‑1 (MQ=255)
                                                                                                 tcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGA‑A‑AGCCCCCTTTTTTCCGTAAAGTTTTTGCCGCGCGCTCTCGGCAACGCACGCCCCCCC‑CCCCCCCCCGCGATGCAACGTTGGTATCACC‑gccgtt                                                     >  2:244929/1‑144 (MQ=255)
                                                                                                 tcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGA‑A‑AGCCCCCCTTTTTGAAAAATTTTTTCGCCCCCCGCTGGAGCCCGCCCCCGCCCCCCC‑CCCCACCCCCCGGCGCAACGGTTGGATTCCC‑CTTGATcc                                                   >  2:1461826/1‑150 (MQ=255)
                                                                                                  cgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGA‑A‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCT‑CGCCACCGCGCGATGCAACGTTGGTATCGCC‑GTTGATCg                                                   <  1:808606/150‑1 (MQ=255)
                                                                                                  cgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGA‑A‑TGCCCCCCTTTTTCAATAATGTTTTCGCCCCCCGCTGTAGCCAGCGCAGGCCCCCCC‑CCCCCCCCCGCGATGCAACGTTGGTATCGCC‑GTTGATCGc                                                  >  1:1288030/1‑151 (MQ=255)
                                                                                                   gATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGA‑A‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCT‑CGCCACCGCGCGATGCAACGTTGGTATCGCC‑GTTGATCGcc                                                 <  2:258574/151‑1 (MQ=255)
                                                                                                           cccTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGA‑A‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCT‑CGCCACCGCGCGATGCAACGTTGGTATCGCC‑GTTGATCGCCATCAGTa                                          <  1:730059/150‑1 (MQ=255)
                                                                                                            ccTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGA‑A‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCT‑CGCCACCGCGCGATGCAACGTTGGTATCGCC‑GTTGATCGCCATCAGTAgg                                        <  2:9617/151‑1 (MQ=255)
                                                                                                              tGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGGTGA‑A‑TGCCCCCCTTTTTTAAAAATGTTTTCGCCCCGCCCTGGAGCCAGCGCAGGGCCCCCC‑CCCCACCCCGCGGTGCAACGGTGGTATTCCC‑GCTGGTCGCCCACCGTAgg                                        >  2:1056719/1‑149 (MQ=255)
                                                                                                                 cGAGCTGCCGGAGATAATCGAGATCGGCGGGGGTTCGA‑A‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCT‑CGCCACCGCGCGATGCAACGTTGGTATCGCC‑GTTGATCGCCATCAGTAGGAGCa                                    <  1:314348/150‑1 (MQ=255)
                                                                                                                 cGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGA‑A‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCT‑CGCCACCGCGCGATGCAACGTTGGTATCGCC‑GTTGATCGCCATCAGTAGGAGCa                                    <  2:619291/150‑1 (MQ=255)
                                                                                                                   aaCTGCCGGAGATAATCGAGATCGGCGGGGGTTCGA‑A‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCT‑CGCCACCGCGCGATGCAACGTTGGTATCGCC‑GTTGATCGCCATCAGTAGGAGCAAgg                                 <  2:1288265/151‑1 (MQ=255)
                                                                                                                       gCCGGAGATAATCGAGATCGGCGGGGGTTCGA‑A‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCT‑CGCCACCGCGCGATGCAACGTTGGTATCGCC‑GTTGATCGCCATCAGTAGGAGCAAGGTAt                              <  2:1160784/150‑1 (MQ=255)
                                                                                                                       gCCGGAGATAATCGAGATCGGCGGGGGTTCGA‑A‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCT‑CGCCACCGCGCGATGCAACGTTGGTATCGCC‑GTTGATCGCCATCAGTAGGAGCAAGGTATc                             <  1:1436719/151‑1 (MQ=255)
                                                                                                                       gCCGGAGATAATCGAGATCGGCGGGGGTTCGA‑A‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCT‑CGCCACCGCGCGATGCAACGTTGGTATCGCC‑GTTGATCGCCATCAGTAGGAGCAAGGTATc                             <  1:1362748/151‑1 (MQ=255)
                                                                                                                        ccGGAGATAATCGAGATCGGCGGGGGTTCGA‑A‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCT‑CGCCACCGCGCGATGCAACGTTGGTATCGCC‑GTTGATCGCCATCAGTAGGAGCAAGGTATCg                            <  2:1184872/151‑1 (MQ=255)
                                                                                                                        ccGGAGATAATCGAGATCGGCGGGGGTTCGA‑A‑AAACCCCCTTTTGTAATAAATTTTTCTGCCCGCGCCGTAGCCCACGCAAGCCCCCCT‑CACCACCCCCCGGTGCCACGGTGGTATCGCCGGTTGTTC‑CCATCCGCAGGAGCGAGGTATCg                            >  2:456922/1‑151 (MQ=255)
                                                                                                                        ccGGAGATAATCGAGATCGGCGGGGGGTCGA‑A‑AGCCCCCCTTTTTCAATAATTTTTTCGCCCCGCGCTGTAGCCAGCGCAGGGCCCCCC‑CGCCACCCCGCGATGCAACGTTGGGATCGCC‑GTTGATCGCCATCAGGAGGAGGAAGGGATCg                            >  1:954096/1‑151 (MQ=255)
                                                                                                                           gagaTAATCGAGATCGGCGGCGGGTCGA‑A‑AGCCCCCCTTTTTCAATAATATTTTCGCCCCCCGCCGGAGCCAGGGCCCGCCCCCCTCCCCCCCCCCCCCAGCCAACGTTGGTTTGCCC‑GTTGA‑CGCCACCAGTAgggggaaaggatccggc                         >  2:1244872/1‑136 (MQ=255)
                                                                                                                             gaTAATCGAGATCGGCGGGGGTTCGA‑ATTGCCCCCTTTTT‑TAATAAAGTTTTTGCCTCGCGCCGTAGCCAGCGCAGGCCCCCCC‑CCCCCCCCCCCGGTGCAACGGTGGTATCGCC‑GTTGGTCGCCCTCAGGAGGAGCAAGGGATTGAgga                        >  2:42028/1‑148 (MQ=255)
                                                                                                                                 aTCGAGATCGGCGGGGGTTCGA‑A‑TGCCCCCTTTTTTCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCC‑CGCCACCCCGCGATGCAACGTTGGGATCGCC‑GTTGATCGCCATCAGTAGGAGCAAGGGATCGAGCAATGcc                   >  1:1309589/1‑151 (MQ=255)
                                                                                                                                   cGAGATCGGCGGGGGTTCGA‑A‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCT‑CGCCACCGCGCGATGCAACGTTGGTATCGCC‑GTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCa                  <  1:607521/150‑1 (MQ=255)
                                                                                                                                    gagaTCGGCGGGGGGTCGA‑A‑TGCCCCCCTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCT‑CGCCACCCCCCGATGCAACGTTGGTATCGCC‑GTTGATCGCCCTCAGTAGGAGCAAGGGATCGAGCAATGCCAGt                >  2:339022/1‑151 (MQ=255)
                                                                                                                                     ngaTCGGCGGGGGTTCGA‑A‑TGCCCCCCTTTTGCAATAATGTTTTCGCCCCGCGCTGTAGCCCGCGCAGGCCCCCCC‑CGCCCCCGCGCGATGCAAAGTTGGTATCGCC‑GTTTGTCGCCCTCCGGAGGAGCAAGGGATCGAGCAATGCCAGtt               >  2:1259935/2‑151 (MQ=255)
                                                                                                                                      gaTCGGCGGGGGTTCGA‑A‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCT‑CGCCACCGCGCGATGCAACGTTGGTATCGCC‑GTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTc              <  2:699160/151‑1 (MQ=255)
                                                                                                                                          ggcggGGGTTCGA‑A‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCT‑CGCCACCGCGCGATGCAACGTTGGTATCGCC‑GTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAgg           >  2:747052/1‑150 (MQ=255)
                                                                                                                                                ggTTCGA‑A‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCT‑CGCCACCGCGCGATGCAACGTTGGTATCGCC‑GTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAAc    <  2:522362/151‑1 (MQ=255)
                                                                                                                                                 gTTCGA‑A‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCT‑CGCCACCGCGCGATGCAACGTTGGTATCGCC‑GTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAAcc   <  1:1448290/151‑1 (MQ=255)
                                                                                                                                                  ttCGA‑A‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCT‑CGCCACCGCGCGATGCAACGTTGGTATCGCC‑GTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAAccc  >  2:859730/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                          
TATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGA‑A‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCT‑CGCCACCGCGCGATGCAACGTTGGTATCGCC‑GTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCC  >  NC_000913/646601‑646897

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