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

TATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACC  >  NC_000913/646601‑646896
                                                                                                                                                   |                                                                                                                                                       
tatTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCgg                                                                                                                                                     >  1:246416/1‑150 (MQ=255)
 atTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCg                                                                                                                                                      >  1:246657/1‑149 (MQ=255)
 atTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGa‑a                                                                                                                                                   <  1:495743/151‑1 (MQ=255)
   ttCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑At                                                                                                                                                  <  2:421737/150‑1 (MQ=255)
   ttCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTTGa‑a                                                                                                                                                   >  1:295832/1‑149 (MQ=255)
        aaTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGA‑ATGccccc                                                                                                                                            >  1:895739/1‑151 (MQ=255)
           ttATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCtt                                                                                                                                          <  2:246421/150‑1 (MQ=255)
               cAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTg                                                                                                                                      <  1:445541/150‑1 (MQ=255)
                aaGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCa                                                                                                                                    <  2:884426/151‑1 (MQ=255)
                 aGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTCG‑‑ATGGCCCCttttt                                                                                                                                       >  2:210014/1‑147 (MQ=255)
                     cTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGA‑ATGCCCCCTTTTTTCaataa                                                                                                                                >  2:697367/1‑150 (MQ=255)
                      ttAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTTGG‑ATGCCCCCcttttttgataa                                                                                                                                >  2:62713/1‑142 (MQ=255)
                       tAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATg                                                                                                                              <  2:289167/150‑1 (MQ=255)
                        aaaTAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGt                                                                                                                             <  1:34576/150‑1 (MQ=255)
                        aaaTAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGA‑ATGCCCCCCTTTTTCTATAAtttt                                                                                                                            >  2:220594/1‑151 (MQ=255)
                        aaaTAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGA‑ATGCCCCCCTTTTTCAATAAtttt                                                                                                                            >  2:692715/1‑151 (MQ=255)
                          ataatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCCG‑ATGCCCCCCTTTTGCCATAATAttt                                                                                                                           >  2:178107/1‑150 (MQ=255)
                               tAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGcc                                                                                                                      <  1:343380/150‑1 (MQ=255)
                                aaaTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGGCGG‑ATGCCCCCCTTTTTCAAAAATGtttttgcccc                                                                                                                    >  1:169559/1‑144 (MQ=255)
                                      gtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGA‑ATGCCCCCCTTTTTCAAAAATGTTTTCGccccccg                                                                                                                 >  1:903161/1‑143 (MQ=255)
                                         cTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTTGG‑ATGCCCCCTTTTTTTAATAATATTTTTGCCTCGCGCTGt                                                                                                             >  1:279211/1‑149 (MQ=255)
                                          tAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTTGG‑AAGCCCCCCTTTTTTAAAAAAGTTTTTGCCCCCCGCTGt                                                                                                             >  1:398883/1‑148 (MQ=255)
                                             aaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGG‑ATGCCCCCCTTTTTCAATAATTTTTTCGCCCCCCGCTGTA‑Gcc                                                                                                        >  2:557892/1‑149 (MQ=255)
                                                 cAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCa                                                                                                  <  1:756007/151‑1 (MQ=255)
                                                 cAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGG‑ATGGCCCCCTTTTTTAATAATGTTTTCGCCCCCCGCCGTA‑Gcc                                                                                                        >  2:821337/1‑145 (MQ=255)
                                                   ggTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGG‑ATGCCCCCCTTTTGCAATAATATTTTCGGCTCCCGCCGTA‑GCCCGGGCAgg                                                                                                >  2:92660/1‑151 (MQ=255)
                                                   ggTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGA‑ATGCCCCCCTTTTTTAATAATGTTTTCGCCTCGCGCCGTA‑GCCCGCGCAg                                                                                                 >  1:466924/1‑150 (MQ=255)
                                                       aGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGA‑ATGCCCCCCTTTTTTATTAATGTTTTTGCCCCGCCCTGTA‑GCCCGCGCAGGcccc                                                                                            >  2:867328/1‑151 (MQ=255)
                                                          aTCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCGGGCCCCCCt                                                                                         <  2:28139/151‑1 (MQ=255)
                                                          aTCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCt                                                                                         <  1:560427/151‑1 (MQ=255)
                                                          aTCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTGGA‑ATGCCCCCCTTTTTCTATAATTTTTTTGCCCCCCCCTGTAGGCCACCCAAGGcccccc                                                                                          >  2:501179/1‑151 (MQ=255)
                                                            cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGA‑ATGGCCCCTTTTTGCAATAATGTTTTTGCCTCGCGCTGTA‑GCCCGCGCAGGcccccccc                                                                                        >  2:888061/1‑148 (MQ=255)
                                                                 ggTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCAcc                                                                                  <  1:820665/151‑1 (MQ=255)
                                                                    cAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACcgcg                                                                               <  2:481517/151‑1 (MQ=255)
                                                                    cAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACcgc                                                                                <  1:611998/150‑1 (MQ=255)
                                                                        aCTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGA‑ATGCCCCCCTTTTTTAATAATTTTTTTGCCCCGCGCCGTA‑GCCAGCGCAGGcccccccccccccccccc                                                                              >  2:363606/1‑136 (MQ=255)
                                                                          tGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTGGACAAGCCCCCTTTTT‑CAATAATGTTTGCGCCCCTCGCTGTA‑GCCAGGCCCGGCCCCCCTCTCCACCGCCCGATg                                                                          >  2:813684/1‑150 (MQ=255)
                                                                           gccgccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATg                                                                          <  2:676324/149‑1 (MQ=255)
                                                                           gccgccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCa                                                                        <  1:271353/151‑1 (MQ=255)
                                                                                  ggACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGtt                                                                   <  1:62711/149‑1 (MQ=255)
                                                                                       gagaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATc                                                             <  2:518353/150‑1 (MQ=255)
                                                                                         gaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGG‑CTGCCCCCCTTTTTTTATAATTTTTTTGCCCCCCGCGGTG‑GCCCGCGCGCGCCCCCCCCCCCCCCCCGCGATGCCAccgtggggt                                                              >  1:880868/1‑139 (MQ=255)
                                                                                          aTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGA‑ATGCCCCCCTTTTTCAATAAAGTTTTCGCCTCCCGCTGTA‑GCCAGCGCAGGCCCCCCCCCCCCCCCCGCGGTGCAACGTTTGTATCGcc                                                          >  2:838765/1‑150 (MQ=255)
                                                                                          aTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCCG‑ATGCCCCCCTTTTTCACTAATGTTTTTGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCCCCCCCCCCCGCGACGCAACGTTGGGATCGcc                                                          >  1:221868/1‑150 (MQ=255)
                                                                                            ttcgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGt                                                        <  1:206931/150‑1 (MQ=255)
                                                                                                 tcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATc                                                   <  1:383894/150‑1 (MQ=255)
                                                                                                      acacTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATc                                             <  2:857697/151‑1 (MQ=255)
                                                                                                        acTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGG‑AAGCCCCCCTTTTTCAATAATTTTTTCGCCTCCCGCTGTA‑GCCAGCGCAGGCCCCCCCCCCCACCCCCCGCTGCAACGTTGGGATTGCCGGTTATCGcc                                                >  2:294332/1‑146 (MQ=255)
                                                                                                           cccTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTa                                         <  1:788704/150‑1 (MQ=255)
                                                                                                            ccTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGC‑ATGGCCCCCTTTTTTAATAATGTTTTTGCCCCCCGCTGGA‑GCCACCGCAGGGCCCCCCCCCCCCCCCGCCGTGCAAAGGTGGGGACGCCCGTGATCGCCATCCGCAgg                                       >  1:202272/1‑151 (MQ=255)
                                                                                                               gTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAg                                     <  2:86992/150‑1 (MQ=255)
                                                                                                               gTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGC‑ATGCCCCCCTTTTTCAATAATGTTTTCGCCTCCCGCTGTA‑GCCAGCGCCGGCCCCCCCCCCCCCCCCCCGGTGCCACGGTGGGATCGCCCTTGATCGCCCTCAATAgg                                       >  2:906490/1‑148 (MQ=255)
                                                                                                                 cGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCa                                   <  2:645070/150‑1 (MQ=255)
                                                                                                                   aaCTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAgg                                <  1:491255/151‑1 (MQ=255)
                                                                                                                     cTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGA‑ATGCCCCCCTTTTTCCATAATATTTTTGCCCCCCGCTGGA‑GCCAGCGCCGGGCCCCCCCCCCCCCCCCCGACGCAACGGTGGTATCGCCGGTGGTCGCCATCAGGAGGGGGAAggg                               >  1:685569/1‑149 (MQ=255)
                                                                                                                        ccGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCg                           <  2:886277/151‑1 (MQ=255)
                                                                                                                         cGGAGATAATCGAGATCGGC‑GGGGGGTTGA‑ATGCCCCCCTTTTTCAATAATGTTTTCGCCCCCCGCTGTA‑GGCCGCGCCGGGCCCCCCCCCCCCCCCCCGATGCCACATTGGGATCGCCCGTGATCGCCCTCAGTAGGAGGAAgg                                >  2:371232/1‑145 (MQ=255)
                                                                                                                             gaTAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCaa                      <  2:192470/151‑1 (MQ=255)
                                                                                                                                        tCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGGATCGAGCAATGCCAGTTc             >  1:645390/1‑149 (MQ=255)
                                                                                                                                          ggc‑ggGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCGAGTTCAGGa         <  1:946058/151‑1 (MQ=255)
                                                                                                                                           gc‑ggGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGAt        <  2:895757/151‑1 (MQ=255)
                                                                                                                                              gggggTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATc       <  1:451558/150‑1 (MQ=255)
                                                                                                                                              gggggTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCt      >  2:242309/1‑151 (MQ=255)
                                                                                                                                              gggggTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCt      <  2:237484/151‑1 (MQ=255)
                                                                                                                                               ggggTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCt      >  2:793597/1‑150 (MQ=255)
                                                                                                                                                  gTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAAcc  >  1:893763/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                       
TATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACC  >  NC_000913/646601‑646896

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