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

AATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCC  >  NC_000913/646599‑646897
                                                                                                                                                     |                                                                                                                                                         
aaTATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGtt                                                                                                                                                          <  1:360864/150‑1 (MQ=255)
   atTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGa                                                                                                                                                      <  1:140158/150‑1 (MQ=255)
     ttCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAt                                                                                                                                                   <  2:1264950/150‑1 (MQ=255)
     ttCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAt                                                                                                                                                   <  2:618633/150‑1 (MQ=255)
      tCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATg                                                                                                                                                  <  2:1242707/150‑1 (MQ=255)
       ccAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCACGTAAGGATCAGTAGGCCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTc                                                                                                                                                         >  2:933682/1‑144 (MQ=255)
           aTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCt                                                                                                                                           <  2:1348080/151‑1 (MQ=255)
              tATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCttt                                                                                                                                         <  2:87924/150‑1 (MQ=255)
                tCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGATGCCCCCcttttt                                                                                                                                      >  2:918959/1‑150 (MQ=255)
                   aGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCaa                                                                                                                                   <  1:949868/151‑1 (MQ=255)
                    gTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCctttt                                                                                                                                       >  1:405855/1‑146 (MQ=255)
                     tGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGAAGGCCCCCttttt                                                                                                                                       >  2:1062351/1‑146 (MQ=255)
                     tGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGCCCCCCTTTTGTaat                                                                                                                                  >  1:1164595/1‑150 (MQ=255)
                      gCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCGCTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCaata                                                                                                                                 <  2:309003/150‑1 (MQ=255)
                      gCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAATGTCCCCCttttttaata                                                                                                                                  >  1:767722/1‑145 (MQ=255)
                         tAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATTCCCCCcttttttcatta                                                                                                                                >  1:1439006/1‑141 (MQ=255)
                          aaaTAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTTCAATAATGtt                                                                                                                            >  2:1623948/1‑151 (MQ=255)
                          aaaTAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGt                                                                                                                             <  1:426189/150‑1 (MQ=255)
                          aaaTAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGACTGACCCCCTTTTTCAAATATGt                                                                                                                             >  1:1016814/1‑150 (MQ=255)
                           aataatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGttt                                                                                                                           <  2:1003619/151‑1 (MQ=255)
                            ataatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGCCCCCCTTTTTTAATAATtttt                                                                                                                           >  1:769738/1‑150 (MQ=255)
                              aatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTgcg                                                                                                                        <  1:1340018/151‑1 (MQ=255)
                              aatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTACCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGCCCCCCTTTTTCAATAAAGttt                                                                                                                           >  1:277956/1‑148 (MQ=255)
                               atTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTgcgc                                                                                                                       <  1:567998/151‑1 (MQ=255)
                                   aaTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTcg                                                                                                                   <  1:1378214/151‑1 (MQ=255)
                                      cTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAAAGCCCCCCTTTTTTAATAATGtt                                                                                                                            >  1:1615420/1‑139 (MQ=255)
                                      cTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGGCCCCCTTTTTCAATAATTTTTTCGcccc                                                                                                                    >  2:1260302/1‑145 (MQ=255)
                                       tgtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAATGCCCCCCTTTTTCAATAAAGTTTTCGcccccc                                                                                                                  >  2:1062887/1‑144 (MQ=255)
                                        gtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGCGATAATCGAGATCGGCGGGGGGGCAAAAGACCCCCTTTTGTTATTATGTTTActcctctccc                                                                                                                >  1:614680/1‑145 (MQ=255)
                                        gtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTcgc                                                                                                                  <  2:1086048/147‑1 (MQ=255)
                                        gtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAATGCCCCCCTTTTTCAATAATTTTTTCGCCCcgcg                                                                                                                 >  2:1070934/1‑148 (MQ=255)
                                            tAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAATGTTTTCGCCCCCCGCCGTAGc                                                                                                          >  2:1151214/1‑151 (MQ=255)
                                             aaaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGcc                                                                                                         <  1:862215/151‑1 (MQ=255)
                                             aaaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGc                                                                                                          <  2:262890/150‑1 (MQ=255)
                                              aaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCa                                                                                                        <  2:1537718/151‑1 (MQ=255)
                                              aaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCa                                                                                                        <  1:862000/151‑1 (MQ=255)
                                              aaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCTTTTTTCAATAATGTTTTCGCCCCGCGCTGTAGcc                                                                                                         >  2:532058/1‑150 (MQ=255)
                                              aaaaCCAGGTAAGGATCAGAAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAACGCCCCCTTTTTCCAATAATGTTTGCACCTCGCTCTTTAGcta                                                                                                        <  1:686037/151‑3 (MQ=255)
                                                    aGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAg                                                                                                  <  1:421749/151‑1 (MQ=255)
                                                       tAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAATGCCCCCCTTTTTTCATAATGTTTTCGCCTCCCGCCGTAGCCAGCGCCGGcc                                                                                               >  2:389947/1‑151 (MQ=255)
                                                           gATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGAAGGCCCCCTTTTTTCATAATATTTTTGCCCCCCGCCGTTGCCCGCGCCGGgccccc                                                                                           >  1:990445/1‑151 (MQ=255)
                                                             tCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCCGAAGGCCCCCTTTTTTAATAATGTTTTGGCCCCGCCCCGGAGCCCGCGCCGGgccccccc                                                                                         >  1:995952/1‑149 (MQ=255)
                                                              cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTc                                                                                         <  2:1545492/150‑1 (MQ=255)
                                                              cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTc                                                                                         <  2:1416962/150‑1 (MQ=255)
                                                               aGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGCATGCCCCCCTTTTTCAATAATGTTTTTGCCCCCCGCCGTAGCCAGCGCCGGcccccccccc                                                                                       >  1:1458371/1‑147 (MQ=255)
                                                                  aGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCAc                                                                                    <  2:1662692/151‑1 (MQ=255)
                                                                  aGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTTAATAATGTTTTCGCCCCCCGCTGTAGCCAGCGCAGGccccccccccccc                                                                                    >  2:610012/1‑144 (MQ=255)
                                                                    gTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACcg                                                                                  <  1:222418/151‑1 (MQ=255)
                                                                     tCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAATTTTTTCGCCCCCCGCTGTAGCCCGCGCCGGgccccccccccccccc                                                                                 >  2:1571852/1‑141 (MQ=255)
                                                                      cAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACcgcg                                                                                <  2:31960/151‑1 (MQ=255)
                                                                      cAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACcgc                                                                                 <  1:1015415/150‑1 (MQ=255)
                                                                      cAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAATGTTTTTGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCCCGCcccccc                                                                                 >  2:1180936/1‑148 (MQ=255)
                                                                        gCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGCATGTCCCCCTTTTTTAATAATATTTTCGCCCCCCGCTGTAGCCAGCGCCGGcccccccccccccccc                                                                                 >  2:504646/1‑138 (MQ=255)
                                                                          aCTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGGAGGCCCCTTTTTTCAATAATGTTTTCGCCTCACGCGGGAGGCAACGCGGGCCCCCCCCGccgcaccc                                                                                >  2:467274/1‑141 (MQ=255)
                                                                            tGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCTTTTTTCAATAATGTTTGCGCCTCCCCCTGTAGCCAGCGCCGGCCCCCCCCCCCACCCCGCGCTg                                                                           >  2:156476/1‑150 (MQ=255)
                                                                               cgccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGCATGCCCCCTCTTTGCAATAAATTTTTTTCCTCCCGCCGGAGCCAGCGGCGGCCCCCCCCCCCAcccccccacgcca                                                                        >  1:1643478/1‑139 (MQ=255)
                                                                               cgccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATTCCCCCCTTTTTCAATAATGTTTGCGGCTCGCCCCGTAGCCCGCCCAGGCCCCCCCCCCCCCcgcccccagcca                                                                        >  2:613332/1‑141 (MQ=255)
                                                                                    ggACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTgg                                                                  <  1:211967/151‑1 (MQ=255)
                                                                                    ggACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTgg                                                                  <  2:1451747/151‑1 (MQ=255)
                                                                                       cTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGAAGCCCCCCTTTTTCAATAATATTTTGGGCCCCCGCTGGAGGCAGCGCAGGGCccccccgcccccccgc                                                                               >  1:517626/1‑135 (MQ=255)
                                                                                            aTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCg                                                          <  2:292180/151‑1 (MQ=255)
                                                                                               tcgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGtt                                                        <  2:820265/150‑1 (MQ=255)
                                                                                               tcgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAATGCCCCCCTTTTTCAAAAATGTTTTCGCCTCCCGCCGTAGCCAGCGCAGGCCCCCCCCCCCACCCCCCCCAGCAACGTTGGGATCGcccctg                                                       >  2:885826/1‑147 (MQ=255)
                                                                                                cgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTg                                                       <  1:1118528/150‑1 (MQ=255)
                                                                                                 gttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGAt                                                     <  2:294550/151‑1 (MQ=255)
                                                                                                 gttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATTCCCCCCTTTTTTCAAAATTATTTTCCCCCGCCGTGCAGCCCGAGCAGGCCCCCCCCCCCCCCGCGCCCTGCAAcactgtgatccccg                                                          >  1:682106/1‑133 (MQ=255)
                                                                                                       tACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAAAGCCCCCCTTTTTCAATAATGTTTTCGCCCCCCCCTGTAGCCAGCGCCGGCCCCCCCCCCCCCCCCGCGATGCAACGCTGGGATCGCCCTTGATCGcc                                                 >  1:606499/1‑149 (MQ=255)
                                                                                                        acacTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCAt                                               <  1:1184303/150‑1 (MQ=255)
                                                                                                        acacTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGATGCCCCCTTTTTTCCATAATATTTTTGCCCCCCCCTGTAGCCAGCGCAGGGCCCCCCCCCCACCCCCCGGTGCAACGTTGGGGTTGCCGCTGATCGCCCTc                                              >  1:1646288/1‑151 (MQ=255)
                                                                                                         cacTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCa                                             <  1:1685561/151‑1 (MQ=255)
                                                                                                          acTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATc                                              <  1:798121/149‑1 (MQ=255)
                                                                                                                 gTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAATGTTTTTGCCCCGCGCTGTAGCCAGCGCAGGCCCCCCCCCCCCCCCCGCGGAGCAACGTTGGGATCGCCGCTGATTGCCCTCAGTa                                          >  2:52997/1‑146 (MQ=255)
                                                                                                                     aaCTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAgg                                 <  1:871512/151‑1 (MQ=255)
                                                                                                                     aaCTGCCGGAGATAATCGAGATCGGCGGGGGGTGGAAAGGCCCCTTTTTTCAAAAAAGTTTTTGCCCCGCCCCGTAGGAAAAGGAGGGCCCCCCCGCCAACCGGGGAGGGAAAGGTGGGAAAGGCGGTGATAGCCATCAGAAGGAGAAAgg                                 >  2:713592/1‑151 (MQ=255)
                                                                                                                      aCTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTGTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGGCATCAGTAGGAGCAAGGt                                <  1:1613522/151‑1 (MQ=255)
                                                                                                                         gCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATc                             <  1:838983/151‑1 (MQ=255)
                                                                                                                         gCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAAAGTTTGTGCCCCCCGCCGTTGCCCGGGCCGGCCCCCCCCGCCCCCCCGCGCCGCGACGTTTGTATCGGCCTTGGTCGCCCTCCGTAGGAGCCAggggac                             >  1:25692/1‑147 (MQ=255)
                                                                                                                           cGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGa                           <  1:1595313/151‑1 (MQ=255)
                                                                                                                           cGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGa                           <  2:780303/151‑1 (MQ=255)
                                                                                                                            ggAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCg                            <  2:1696811/149‑1 (MQ=255)
                                                                                                                              agaTAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCa                        <  2:707174/151‑1 (MQ=255)
                                                                                                                                   aTCGAGATCGGCGGGGGGTCGAATGCCCCCTTTTTTTAATAAATTTTTTGCCTCGCGCTGTAGCCAGCGCAGGGCCCCCCCCCCACCCCGCGCTGCAAAGGTGGGTTCGCCCGTGGTCGCCCCCCGGAGGAGGAAGGGATCggggca                       >  1:209595/1‑142 (MQ=255)
                                                                                                                                             gcggGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGAt         <  2:405870/151‑1 (MQ=255)
                                                                                                                                             gcggGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGAt         <  2:1084255/151‑1 (MQ=255)
                                                                                                                                                 ggggTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCt       >  2:631156/1‑150 (MQ=255)
                                                                                                                                                    gTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAAcc   >  2:721389/1‑151 (MQ=255)
                                                                                                                                                     ttCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAAccc  >  2:843892/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                         
AATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCC  >  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: