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

ATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACC  >  NC_000913/646602‑646896
                                                                                                                                                 |                                                                                                                                                       
atTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGa                                                                                                                                                    <  1:598043/150‑1 (MQ=255)
  ttCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAt                                                                                                                                                  <  1:535062/150‑1 (MQ=255)
      aaaTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGGTGAATGcccc                                                                                                                                            >  1:871612/1‑151 (MQ=255)
       aaTTTATCAAGTGCTTAAATAATTATATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTTAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCCAATAgcccc                                                                                                                                           >  1:465335/1‑151 (MQ=255)
       aaTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGccccc                                                                                                                                           >  2:291728/1‑151 (MQ=255)
        aTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAATGccccc                                                                                                                                           >  2:877386/1‑150 (MQ=255)
        aTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAAAGccccc                                                                                                                                           >  2:321398/1‑150 (MQ=255)
         tttATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGAAGccccc                                                                                                                                           >  2:446684/1‑149 (MQ=255)
           tATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGGAAGACCCCcttt                                                                                                                                       >  2:734752/1‑151 (MQ=255)
           tATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCttt                                                                                                                                        <  2:531674/150‑1 (MQ=255)
           tATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCttt                                                                                                                                        <  1:862558/150‑1 (MQ=255)
           tATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCtttt                                                                                                                                       >  2:177685/1‑151 (MQ=255)
            aTCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGCATGGCCCCttttt                                                                                                                                      >  2:84237/1‑151 (MQ=255)
            aTCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAGTGCCCCCttttt                                                                                                                                      >  1:918979/1‑151 (MQ=255)
            aTCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCctttt                                                                                                                                      >  1:270112/1‑151 (MQ=255)
              cacgTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGATGTCCCCCtttttt                                                                                                                                     >  1:244375/4‑150 (MQ=255)
               aaGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGACTGCCCCCctttt                                                                                                                                      >  2:505838/1‑148 (MQ=255)
                      tAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATg                                                                                                                             <  1:576485/150‑1 (MQ=255)
                       aaaTAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAATGtt                                                                                                                           >  2:62957/1‑151 (MQ=255)
                        aataatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGttt                                                                                                                          <  2:123452/151‑1 (MQ=255)
                          taatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGCCCCCCTTTTTCAATAATGttt                                                                                                                          >  2:847152/1‑149 (MQ=255)
                          taatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGGCGGATTCCCCCCTTTTTCaataaatttt                                                                                                                          >  1:300301/1‑144 (MQ=255)
                           aatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTgcg                                                                                                                       <  1:564731/151‑1 (MQ=255)
                           aatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTgc                                                                                                                        <  1:371180/150‑1 (MQ=255)
                              tAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCTGATGTACCCCTTTTGCAATAATAttttgtcc                                                                                                                     >  1:233166/1‑145 (MQ=255)
                                aaTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTcg                                                                                                                  <  2:64571/151‑1 (MQ=255)
                                aaTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGATGCCCTCTTTTTTCcataaatatttttgctc                                                                                                                   >  2:857835/1‑138 (MQ=255)
                                   cTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAATGTTTTCGcccc                                                                                                                   >  2:751955/1‑145 (MQ=255)
                                        cTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCTTTTTTTTATAATGtttttgcccccc                                                                                                                 >  2:422467/1‑135 (MQ=255)
                                             aaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAgc                                                                                                     <  1:753742/151‑1 (MQ=255)
                                               ccAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTGGAATGCCCCCCTTTTTGAATAATTTTTTGGCCTCCCCCCGTAGCCAgcg                                                                                                    >  1:331479/1‑150 (MQ=255)
                                               ccAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAATGTTTTCGCCCCGCGCTGTAGcc                                                                                                        >  2:932291/1‑146 (MQ=255)
                                                      aGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTTTAACACCCCCCTTTTTCCCTAATTTTTTTGCCTCGCGCCGGAGCCCGCCCACGcccc                                                                                            >  1:581233/1‑151 (MQ=255)
                                                        gATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTGAATAATGTTTTCGGCTCGCGCGGTTGGCAGCGCCGGgccccc                                                                                          >  2:69486/1‑151 (MQ=255)
                                                           cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGTATGCCCCCCTTTTTCAATAATGTTTTCGCCCCCCGCCGGTGGCCGCGCAGGcccccccc                                                                                        >  2:449524/1‑148 (MQ=255)
                                                           cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGATGCCCCCCTTTTTCAATAATATTTTTGCCCCCCGCTGTAGCCCGCGCAGGcccccccc                                                                                        >  2:305491/1‑148 (MQ=255)
                                                              tAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGCTCGGATGCCCCCCTTTTTGACAAATGTTTTGGCCCCCCGCCGGTGGCAACGCAGGGCCCCCTCGcc                                                                                     >  1:710931/1‑150 (MQ=255)
                                                               aGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCa                                                                                    <  1:259611/150‑1 (MQ=255)
                                                                    aGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACcgcgc                                                                              <  2:780427/151‑1 (MQ=255)
                                                                              ccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGt                                                                    <  2:433894/151‑1 (MQ=255)
                                                                              ccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGt                                                                    <  1:119111/151‑1 (MQ=255)
                                                                                 ggACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTg                                                                  <  1:593020/150‑1 (MQ=255)
                                                                                      gagaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATc                                                             <  1:544211/150‑1 (MQ=255)
                                                                                             cgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTg                                                      <  1:776086/150‑1 (MQ=255)
                                                                                             cgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAAAGCCCCCCTTTTTCAATAATGTTTTCGCCCCGCGCCGTAGCCAGCGCAGGCCCCCCCCCCCCCCCCCCGATGCAACGTTGGTATCGcc                                                          >  2:140385/1‑146 (MQ=255)
                                                                                                     acacTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTGCGGATGACCCCCTTTTTCCATAAAGATTTCGCCTCCCCCAGTAGGCCGCGCAGGCACCCCCCCCCCCCCCCCCAGGCCGCCGAGGGTTCGCCGGTGATCGCCATc                                             >  1:15313/1‑151 (MQ=255)
                                                                                                      cacTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGATGCCCCCCTTTTTTCATAATGTTTTTGCCCCGCGCTGTAGCCCGCGCAGGCCCCCCCCCCCCCCCCCCGGTGCAACGATGGTATCGCCGGTGGTCGCCATc                                             >  2:824976/1‑150 (MQ=255)
                                                                                                          cccTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTa                                         <  2:908284/150‑1 (MQ=255)
                                                                                                          cccTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCTGATTGCCCCTTTTTTGAAAAATGTTTTTTGCTCCCGCTTGAGCCAGCCCAAGCCCCCCCCCCCACCCCGCGCTGCCACGTTGGTGTCGCCGTTGGTCg                                                  >  1:631590/1‑141 (MQ=255)
                                                                                                          cccTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGTATGCCCCCCTTTTTCTATAATGTTTTCGCCCCGCGCTTGTGCCAGGGCCAGGCCCCCTCCCCCCCCCGCGGTGCAAAGGTGGTATCCCCCCTTTTCGCCACCAGt                                          >  2:529061/1‑149 (MQ=255)
                                                                                                            cTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATTCCCCCTTTTTTCAATAATGTTTTCGCCCCGCGCTGTAGGCCGCGCCGGGCCCCCCCGCCCCCCCCCGGTTCAAAGGTGGTATCGCCGCTGGTCGCCATCACTAgg                                       >  2:905672/1‑150 (MQ=255)
                                                                                                                cGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCaa                                  <  2:460477/151‑1 (MQ=255)
                                                                                                                cGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGACCCCCTTTTGCAATAATGTTTGTGCCTCCCGCTGGAGCCCGCGCCGGCCCCCCTCGCCCCCCCGCGCTGCAACGTTGGTATCGCCGTCGGTCGCCATCCGAAGGAGc                                    >  1:755034/1‑149 (MQ=255)
                                                                                                                     tGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAgga                               <  1:277218/149‑2 (MQ=255)
                                                                                                                      gCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTa                              <  2:726746/149‑1 (MQ=255)
                                                                                                                        cGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGa                          <  1:692576/151‑1 (MQ=255)
                                                                                                                          gagaTAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGc                        <  1:739053/151‑1 (MQ=255)
                                                                                                                                    agaTCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGtt              <  2:7546/151‑1 (MQ=255)
                                                                                                                                          gcggGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGAt        <  2:932397/151‑1 (MQ=255)
                                                                                                                                           cggGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATc       <  1:234682/151‑1 (MQ=255)
                                                                                                                                                gTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCGGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAAcc  <  2:830091/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                       
ATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACC  >  NC_000913/646602‑646896

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