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

ATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACC  >  NC_000913/646602‑646896
                                                                                                                                                 |                                                                                                                                                       
atTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑Ga                                                                                                                                                    <  1:598043/150‑1 (MQ=255)
  ttCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAAt                                                                                                                                                  <  1:535062/150‑1 (MQ=255)
      aaaTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGGT‑GAATG‑cccc                                                                                                                                            >  1:871612/1‑151 (MQ=255)
       aaTTTATCAAGTGCTTAAATAATTATATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTTAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTC‑CAATA‑gcccc                                                                                                                                           >  1:465335/1‑151 (MQ=255)
       aaTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTC‑GAATG‑ccccc                                                                                                                                           >  2:291728/1‑151 (MQ=255)
        aTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTT‑GAATG‑ccccc                                                                                                                                           >  2:877386/1‑150 (MQ=255)
        aTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTC‑GAAAG‑ccccc                                                                                                                                           >  2:321398/1‑150 (MQ=255)
         tttATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTC‑GGAAG‑ccccc                                                                                                                                           >  2:446684/1‑149 (MQ=255)
           tATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GGAAG‑ACCCCcttt                                                                                                                                       >  2:734752/1‑151 (MQ=255)
           tATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATG‑CCCCCttt                                                                                                                                        <  2:531674/150‑1 (MQ=255)
           tATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATG‑CCCCCttt                                                                                                                                        <  1:862558/150‑1 (MQ=255)
           tATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTC‑GAATG‑CCCCCtttt                                                                                                                                       >  2:177685/1‑151 (MQ=255)
            aTCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GCATG‑GCCCCttttt                                                                                                                                      >  2:84237/1‑151 (MQ=255)
            aTCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTC‑GAGTG‑CCCCCttttt                                                                                                                                      >  1:918979/1‑151 (MQ=255)
            aTCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTC‑GAATG‑CCCCCctttt                                                                                                                                      >  1:270112/1‑151 (MQ=255)
              cacgTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTT‑GGATGTCCCCCtttttt                                                                                                                                     >  1:244375/4‑150 (MQ=255)
               aaGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTC‑GACTG‑CCCCCctttt                                                                                                                                      >  2:505838/1‑148 (MQ=255)
                      tAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATg                                                                                                                             <  1:576485/150‑1 (MQ=255)
                       aaaTAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTC‑GAATG‑CCCCCCTTTTTCAATAATGtt                                                                                                                           >  2:62957/1‑151 (MQ=255)
                        aataatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGttt                                                                                                                          <  2:123452/151‑1 (MQ=255)
                          taatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTC‑GGATG‑CCCCCCTTTTTCAATAATGttt                                                                                                                          >  2:847152/1‑149 (MQ=255)
                          taatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGGC‑GGATT‑CCCCCCTTTTTCaataaatttt                                                                                                                          >  1:300301/1‑144 (MQ=255)
                           aatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTgcg                                                                                                                       <  1:564731/151‑1 (MQ=255)
                           aatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTgc                                                                                                                        <  1:371180/150‑1 (MQ=255)
                              tAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTC‑TGATG‑TACCCCTTTTGCAATAATAttttgtcc                                                                                                                     >  1:233166/1‑145 (MQ=255)
                                aaTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTcg                                                                                                                  <  2:64571/151‑1 (MQ=255)
                                aaTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTT‑GGATG‑CCCTCTTTTTTCcataaatatttttgctc                                                                                                                   >  2:857835/1‑138 (MQ=255)
                                   cTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTC‑GAATG‑CCCCCCTTTTTCAATAATGTTTTCGcccc                                                                                                                   >  2:751955/1‑145 (MQ=255)
                                        cTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTC‑GAATG‑CCCCCTTTTTTTTATAATGtttttgcccccc                                                                                                                 >  2:422467/1‑135 (MQ=255)
                                             aaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAgc                                                                                                     <  1:753742/151‑1 (MQ=255)
                                               ccAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTG‑GAATG‑CCCCCCTTTTTGAATAATTTTTTGGCCTCCCCCCGTAGCCAgcg                                                                                                    >  1:331479/1‑150 (MQ=255)
                                               ccAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTC‑GAATG‑CCCCCCTTTTTCAATAATGTTTTCGCCCCGCGCTGTAGcc                                                                                                        >  2:932291/1‑146 (MQ=255)
                                                      aGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTT‑TAACA‑CCCCCCTTTTTCCCTAATTTTTTTGCCTCGCGCCGGAGCCCGCCCACGcccc                                                                                            >  1:581233/1‑151 (MQ=255)
                                                        gATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTC‑GAATG‑CCCCCCTTTTTGAATAATGTTTTCGGCTCGCGCGGTTGGCAGCGCCGGgccccc                                                                                          >  2:69486/1‑151 (MQ=255)
                                                           cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTT‑GTATG‑CCCCCCTTTTTCAATAATGTTTTCGCCCCCCGCCGGTGGCCGCGCAGGcccccccc                                                                                        >  2:449524/1‑148 (MQ=255)
                                                           cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTT‑GGATG‑CCCCCCTTTTTCAATAATATTTTTGCCCCCCGCTGTAGCCCGCGCAGGcccccccc                                                                                        >  2:305491/1‑148 (MQ=255)
                                                              tAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGCTC‑GGATG‑CCCCCCTTTTTGACAAATGTTTTGGCCCCCCGCCGGTGGCAACGCAGGGCCCCCTCGcc                                                                                     >  1:710931/1‑150 (MQ=255)
                                                               aGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCa                                                                                    <  1:259611/150‑1 (MQ=255)
                                                                    aGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACcgcgc                                                                              <  2:780427/151‑1 (MQ=255)
                                                                              ccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGt                                                                    <  2:433894/151‑1 (MQ=255)
                                                                              ccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGt                                                                    <  1:119111/151‑1 (MQ=255)
                                                                                 ggACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTg                                                                  <  1:593020/150‑1 (MQ=255)
                                                                                      gagaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATc                                                             <  1:544211/150‑1 (MQ=255)
                                                                                             cgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTg                                                      <  1:776086/150‑1 (MQ=255)
                                                                                             cgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTC‑GAAAG‑CCCCCCTTTTTCAATAATGTTTTCGCCCCGCGCCGTAGCCAGCGCAGGCCCCCCCCCCCCCCCCCCGATGCAACGTTGGTATCGcc                                                          >  2:140385/1‑146 (MQ=255)
                                                                                                     acacTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTGC‑GGATG‑ACCCCCTTTTTCCATAAAGATTTCGCCTCCCCCAGTAGGCCGCGCAGGCACCCCCCCCCCCCCCCCCAGGCCGCCGAGGGTTCGCCGGTGATCGCCATc                                             >  1:15313/1‑151 (MQ=255)
                                                                                                      cacTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTT‑GGATG‑CCCCCCTTTTTTCATAATGTTTTTGCCCCGCGCTGTAGCCCGCGCAGGCCCCCCCCCCCCCCCCCCGGTGCAACGATGGTATCGCCGGTGGTCGCCATc                                             >  2:824976/1‑150 (MQ=255)
                                                                                                          cccTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTa                                         <  2:908284/150‑1 (MQ=255)
                                                                                                          cccTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCTGATTG‑CCCCTTTTTTG‑AAAAATGTTTTTTGCTCCCGCTTGAGCCAGCCCAAGCCCCCCCCCCCACCCCGCGCTGCCACGTTGGTGTCGCCGTTGGTCg                                                  >  1:631590/1‑141 (MQ=255)
                                                                                                          cccTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTC‑GTATG‑CCCCCCTTTTTCTATAATGTTTTCGCCCCGCGCTTGTGCCAGGGCCAGGCCCCCTCCCCCCCCCGCGGTGCAAAGGTGGTATCCCCCCTTTTCGCCACCAGt                                          >  2:529061/1‑149 (MQ=255)
                                                                                                            cTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTC‑GGATT‑CCCCCTTTTTTCAATAATGTTTTCGCCCCGCGCTGTAGGCCGCGCCGGGCCCCCCCGCCCCCCCCCGGTTCAAAGGTGGTATCGCCGCTGGTCGCCATCACTAgg                                       >  2:905672/1‑150 (MQ=255)
                                                                                                                cGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCaa                                  <  2:460477/151‑1 (MQ=255)
                                                                                                                cGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTC‑GAATG‑ACCCCCTTTTGCAATAATGTTTGTGCCTCCCGCTGGAGCCCGCGCCGGCCCCCCTCGCCCCCCCGCGCTGCAACGTTGGTATCGCCGTCGGTCGCCATCCGAAGGAGc                                    >  1:755034/1‑149 (MQ=255)
                                                                                                                     tGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAgga                               <  1:277218/149‑2 (MQ=255)
                                                                                                                      gCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTa                              <  2:726746/149‑1 (MQ=255)
                                                                                                                        cGGAGATAATCGAGATCGGCGGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGa                          <  1:692576/151‑1 (MQ=255)
                                                                                                                          gagaTAATCGAGATCGGCGGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGc                        <  1:739053/151‑1 (MQ=255)
                                                                                                                                    agaTCGGCGGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGtt              <  2:7546/151‑1 (MQ=255)
                                                                                                                                          gcggGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGAt        <  2:932397/151‑1 (MQ=255)
                                                                                                                                           cggGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATc       <  1:234682/151‑1 (MQ=255)
                                                                                                                                                gTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCGGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAAcc  <  2:830091/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                       
ATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACC  >  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: ‑