Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,107,8780TG42.4% 82.6 / 21.7 61G109G (GGT→GGG) pfkA6‑phosphofructokinase I
Reads supporting (aligned to +/- strand):  ref base T (6/28);  new base G (25/0);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 4.54e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACATGATCAACCGTGGCGGTACGTTCCTCGGTTCTGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTA‑CA‑T‑‑GGGTGCAATGCGTC‑TGA‑CC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCT  >  NC_000913/4107730‑4108028
                                                                                                                                                        |                                                                                                                                                          
aCATGATCAACCGTGGCGGTACGTTCCTCGGTTCTGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTA‑CA‑T‑‑GGGTGc                                                                                                                                                          <  2:754145/151‑1 (MQ=255)
               gCGGTACGTTCCTCGGTTCTGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCGGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTAGCGGCGTTGACGGTTCGTA‑CA‑T‑‑GGGTGCAATGCGtgt                                                                                                                                                 <  2:467789/145‑3 (MQ=255)
                    aCGTTCCTCGGTTCTGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGGGACGGGTCCTT‑CA‑T‑‑GGGGGCAATGCGTC‑TGG‑CC‑G‑CAAAgg                                                                                                                                  >  1:633897/1‑151 (MQ=255)
                     cGTTCCTCGGTTCTGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTA‑CA‑T‑‑GGGTGCAATGCGTC‑TGA‑CC‑G‑AAATggg                                                                                                                                 <  1:1065649/151‑1 (MQ=255)
                            cGGTTCTGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGGGGCGGGTCCTACCA‑T‑‑TGG‑GGAAAGCGGC‑TGA‑CC‑G‑GAAAGGGGGTcccc                                                                                                                          >  1:329915/1‑150 (MQ=255)
                                   gcgcGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTA‑CA‑T‑‑GGGTGCAATGCGTC‑TGA‑CC‑G‑AAATGGGCTTCCCGTGCATCg                                                                                                                   <  2:666745/151‑1 (MQ=255)
                                     gcgTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGACAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGGGAAGGTTCCTA‑TC‑T‑‑GGGGGCAATGCGTC‑TTA‑CC‑G‑AAAAGGGAGTccc                                                                                                                           >  2:812544/1‑141 (MQ=255)
                                       gTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGGCGGGTTCTT‑CA‑T‑‑TGGGGGAATGCGGCTTGC‑CC‑G‑AAAGGGGCTcccctgcctccg                                                                                                                   >  1:473320/1‑140 (MQ=255)
                                       gTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTA‑CA‑T‑‑GGGTGCAATGCGTC‑TGA‑CC‑G‑AAATGGGCTTCCCGTGCATCGGTCt                                                                                                               <  1:301967/151‑1 (MQ=255)
                                           cccGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTA‑CA‑T‑‑GGGTGCAATGCGTC‑TGA‑CC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGcc                                                                                                            <  1:385942/150‑1 (MQ=255)
                                            ccGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTA‑CA‑T‑‑GGGTGCAATGCGTC‑TGA‑CC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCgg                                                                                                          <  2:1032824/151‑1 (MQ=255)
                                            ccGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGCGACGGTTTCTC‑CA‑T‑‑TGGTGGGATGCGGC‑TCG‑CC‑G‑CAAAGGGGTTgcccccccgccgcgt                                                                                                               >  2:318833/1‑134 (MQ=255)
                                                aaTTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTA‑CA‑T‑‑GGGTGCAATGCGTC‑TGA‑CC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCGGGCAc                                                                                                      <  1:783942/151‑1 (MQ=255)
                                                    ccGCGACGAGAGCATCCGCGCGGTGCCTATGGAACACCTGAAAAAACGTTGTATCGTCGCGCTGGTGGTTATCAGCGGTGACGGTTCCTA‑CA‑T‑‑GGCTGCAATGCGTC‑TGA‑GC‑G‑AAATGGTCTTCCCGTGCATCGGTCTGCCGTGCACTATc                                                                                                  <  2:174850/151‑1 (MQ=255)
                                                    ccGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTTACGGGTCCTC‑CA‑TC‑GGGGGCAA‑GTGTC‑TGACCC‑G‑AAAGGGGC‑TCCCCTGCCTCgggtttccgg                                                                                                          >  1:129455/1‑135 (MQ=255)
                                                    ccGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGGTTCTT‑CC‑T‑‑GGGGGCCATGCGTC‑TGA‑CC‑G‑CAAAGGGGTTCCCCTGCCTCGGTCTGCCGGGccccat                                                                                                   >  1:1066789/1‑145 (MQ=255)
                                                      gcgACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTA‑CA‑T‑‑GGGTGCAATGCGTC‑TGA‑CC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGa                                                                                                <  2:871350/151‑1 (MQ=255)
                                                       cgacgaGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTA‑CA‑T‑‑GGGTGCAATGCGTC‑TGA‑CC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGAc                                                                                               <  2:1048079/151‑1 (MQ=255)
                                                        gacgaGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGGCGGGTCCTA‑CC‑T‑‑TGGGGCCATGCGTC‑CGG‑CCGG‑AAAGGGGGTTCCC‑CGCCTCGGTCTGCCggggcctct                                                                                                   >  2:734465/1‑140 (MQ=255)
                                                        gacgaGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGGTCCTA‑CC‑A‑‑GGGGGGCATGGGGC‑TGG‑CC‑G‑ACAAGGGGGGCCCGTGCATCGGTGTGCCcgggcccat                                                                                                   >  2:418613/1‑140 (MQ=255)
                                                              aaCATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGGCGGATGCTA‑CC‑T‑‑TGGGGCCAAGCGGC‑TGG‑GC‑GAAAAGGGGCTTCCCGGGC‑TCCGTCTGGCCGGCCCGATccatca                                                                                             >  2:326059/1‑141 (MQ=255)
                                                                  tacgcgCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGGGGGGTGCTA‑AA‑TTGGGGGGAAACGTGTC‑TG‑‑CG‑G‑AAAGGGGGTTCCCGGGCTCCGGGCTG‑CGGGGGCTATAGACAAcaccaac                                                                                      >  2:112479/3‑143 (MQ=255)
                                                                     gcgcCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGAAGGGTTCTA‑CA‑T‑‑GGGGGCAAAGGGTC‑TGG‑CC‑G‑AAAAGGGGGTTCCCTGCCTCCGTCTGGCGGGCACTCTCGACAAAGACATCaaa                                                                                   >  2:1023288/1‑149 (MQ=255)
                                                                        ccGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGCGGCGGATTCTA‑CA‑T‑‑AGGGGCCATGCGTC‑TGC‑CC‑GAAAAGGGGCTTCCGGT‑CATTGGTCTTCCGGGGGATCTCGCCCACGACATCAAAg                                                                                  >  2:77593/1‑147 (MQ=255)
                                                                           tGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGAAGGGTCCCC‑CA‑A‑‑GGGGGGAAAGCGGG‑TGG‑CC‑C‑AAATGGGGGTTCCCCGGATCGGTCTGCCGGGCACCATCGAAAACCACAACaaa                                                                                   >  2:634662/1‑143 (MQ=255)
                                                                            ggCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTA‑CA‑T‑‑GGGTGCAATGCGTC‑TGA‑CC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGa                                                                           <  1:856267/150‑1 (MQ=255)
                                                                              cTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTT‑CA‑T‑‑GGGGGCAATGCGGC‑TGA‑CC‑G‑AAAAAGGGTTCCCCCGCCTCGGTCTGCCCGGGCCTATCGACAAAAACAACaaa                                                                                   >  1:541294/1‑140 (MQ=255)
                                                                               tATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGCGGAGGGTTCTC‑CA‑T‑‑GGGGGGGAAGCGTC‑CGA‑CC‑G‑GAAAGGGGTTCCCCGGCCTCGGGCTGCCCGGCACTATCGGCAACGACATCCAAAGccccgcatac                                                                       >  2:9631/1‑142 (MQ=255)
                                                                                aTCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTA‑CA‑T‑‑GGGTGCAATGCGTC‑TGA‑CC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTa                                                                        <  1:556051/149‑1 (MQ=255)
                                                                                aTCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTA‑CA‑T‑‑GGGTGCAATGCGTC‑TGA‑CC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTa                                                                        <  1:917782/149‑1 (MQ=255)
                                                                                 tCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGGTTCCA‑CA‑A‑‑GGGGGGAACGCGGC‑TGA‑CC‑G‑GAATGGGGTTCCCGCGCCTCGGGCTGCCGGGGGCAATCTACAACGACCTCAAAAGCAcc                                                                             >  2:75041/1‑142 (MQ=255)
                                                                                   gAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGGTCCTA‑AA‑T‑‑TGGGGCAATGCGGCTTAA‑CC‑C‑AAAGGGGC‑TCCCCTGCCTCGGGGCGCCCGGCGCCATCGACAACGACAACAAAgggaacgccatcac                                                                     >  2:799517/1‑137 (MQ=255)
                                                                                    aaaaCCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGGCGGGTTCTT‑CATT‑‑GGGGGCAA‑GCGTC‑TGA‑CC‑C‑AAAAGGGGTTTCCCTGCATCGGGCTGCCGGGGACTATCGACAACGACATCAAAAGCACCGACTaca                                                                      >  1:542242/1‑147 (MQ=255)
                                                                                     aaaCCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTT‑CA‑T‑‑GGGTGCAATGCGTC‑TGA‑CC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCTGGCACTATCGACAACGACATCAAAGGCACTGATTACACTa                                                                   <  2:879665/149‑1 (MQ=255)
                                                                                          tGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTA‑CA‑T‑‑GGGGGGAAAGCGTC‑TGA‑CC‑G‑GAAAAGGCTTCCCCTGCATCGGGCTTCCGGGGGCTATCGACAACGACATCAAAAGCACTGACTCCACTATCGGttt                                                            >  1:200716/1‑151 (MQ=255)
                                                                                           gAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTA‑CA‑T‑‑GGGTGCAATGCGTC‑TGA‑CC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGtt                                                             <  1:411859/149‑1 (MQ=255)
                                                                                           gAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTA‑CA‑T‑‑GGGTGCAATGCGTC‑TGA‑CC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTttc                                                           <  1:1056309/151‑1 (MQ=255)
                                                                                           gAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGGTTCCC‑CA‑T‑‑GGGGGGAAAGCGTC‑TTA‑CC‑G‑AAAAGGGGCTCCCCTGCCTCGGGCTGCCCGGGGCTATCCACAAAAAAATCAAAAGCACCGACTTCACTATCGGgtt                                                            >  1:393175/1‑150 (MQ=255)
                                                                                            aaaaaaCGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGGTCCTT‑CA‑T‑‑GGGGGCGATGCGGC‑TGA‑CC‑C‑CAAAAGGGTTCCCCCGCCTCCGGGTGGCGGGGGCCATCGACAAACACAACAAAAGCACTGACTACACCATACGttttt                                                          >  1:869673/1‑149 (MQ=255)
                                                                                             aaaaaCGTCGTAATGACTCGCTGTTGGTTGTCCGGGGTGGCCGTTCCTA‑CC‑T‑‑GGCTGCAATCCTCC‑TGA‑CC‑T‑AAATGGGCTTCCCGTGCATCAGTCTGCCGGGCACTATCGACAACTTCATCAAAGTCACTGACTACACTATCGGTttct                                                          <  2:793930/150‑1 (MQ=255)
                                                                                              aaaaCGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGGTCCTA‑CA‑T‑‑GGGGGGAAAGCGTC‑TGA‑CC‑G‑AAATGGGGTTCCCCCGCATCGGTCTGCCGGGCACTATCGCCAACGACACCCAAAGCACTGACCACACCATCGGTttctt                                                         >  1:111765/1‑150 (MQ=255)
                                                                                               aaaCGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGAAGGGTTCTACCA‑T‑‑GGG‑GGAAAGGGTC‑TGA‑CC‑G‑CAAAGGGCTTCCCGTGCCTCGGGCTGCCGGGGACTATCGACAACAACAACCAAAGCACTGACCACACTATCGGTttctt                                                         >  2:296428/1‑149 (MQ=255)
                                                                                                 aCGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGGTCCTA‑CA‑T‑‑TGGGGGAAAGCGGC‑TGG‑CC‑GAAAAGGGGGTTCCCG‑GCCTCCGTCTGGCCGGGACTATCGACCACCACATCAAAAGCACAGCCTACACCATCGGGttttt                                                         >  2:1022762/1‑147 (MQ=255)
                                                                                                 aCGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGGTCCTA‑CA‑T‑‑GGGGGGAATGCGTC‑TGA‑CC‑G‑CAAAGGGCTTCCCCTGCATCGGTCTGCCGGGGACTATCGACAACGACATCAAAAGCACTGACTACACTATCGGTTTCTTTAc                                                      >  1:485329/1‑150 (MQ=255)
                                                                                                        aTCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTA‑CA‑T‑‑GGGTGCAATGCGTC‑TGA‑CC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGa                                              <  1:254480/151‑1 (MQ=255)
                                                                                                          cGACGCGCTGGTGGTTATCGGCGGTGACGGGTCCTC‑CA‑T‑‑GGGGGCAAAGCGGC‑TGA‑CC‑G‑GAAAGGGCTTCCCCCGCCTCGGGCTGCCGGGGCCCATCCACAACGACATCAAAAGGACCGACTACACCATCGGGTTTTTTCCTCCGCTgg                                              >  1:455680/1‑148 (MQ=255)
                                                                                                          cGACGCGCTGGTGGTTATCGGCGGAGCGCGTTTCCC‑CA‑T‑‑CGGGGCAATGGGAG‑CGA‑CC‑G‑GAAAGGGGTTTCCCTGCCTCCGGCCTCCCGGCCCCCTCGACAACGAAATCAAAAGCCCCTACCACACTATCGGTTTCTTCCCtcccc                                                 >  1:698245/1‑142 (MQ=255)
                                                                                                              gcgcTGGTGGTTATCGGCGGTGACGGGTCCTT‑CA‑T‑‑GGGTGGAATGCGTC‑TGG‑CC‑G‑AAAAGGGGGTTCCGTGCATCCGTGTGCCGGGCCCTATTGGCAACCAAAACAAAAGCACTGACTACACTATCGGTTTTTTCTCTTCGCTGGGCa                                           >  2:865730/1‑148 (MQ=255)
                                                                                                                   ggtggtTATCGGCGGTGACGGTTCCTA‑CA‑T‑‑GGGTGCAATGCGTC‑TGA‑CC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGt                                    <  1:492881/150‑1 (MQ=255)
                                                                                                                      ggtTATCGGCGGTGACGGTTCCTA‑CA‑T‑‑GGGTGCAATGCGTC‑TGA‑CC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGa                                 <  1:208730/150‑1 (MQ=255)
                                                                                                                            cggcggTGACGGTTCCTA‑CA‑T‑‑GGGTGCAATGCGTC‑TGA‑CC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGAt                           <  1:831813/150‑1 (MQ=255)
                                                                                                                            cggcggTGACGGTTCCTA‑CA‑T‑‑GGGTGCAATGCGTC‑TGA‑CC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATc                          >  1:72719/1‑151 (MQ=255)
                                                                                                                              gcggTGACGGTTCCTA‑CA‑T‑‑GGGTGCAATGCGTC‑TGA‑CC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGa                        <  1:326258/151‑1 (MQ=255)
                                                                                                                                ggTGACGGTTCCTA‑CA‑T‑‑GGGTGCAATGCGTC‑TGA‑CC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGAcc                      <  1:263948/151‑1 (MQ=255)
                                                                                                                                ggTGACGGTTCCTA‑CA‑T‑‑GGGTGCAATGCGTC‑TGA‑CC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGAcc                      >  2:123074/1‑151 (MQ=255)
                                                                                                                                    aCGGTTCCTA‑CA‑T‑‑GGGTGCAATGCGTC‑TGA‑CC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGt                    <  2:457133/149‑1 (MQ=255)
                                                                                                                                     cGGTTCCTA‑CA‑T‑‑GGGTGCAATGCGTC‑TGA‑CC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTg                 <  1:288223/151‑1 (MQ=255)
                                                                                                                                       gTTCCTA‑CA‑T‑‑GGGTGCAATGCGTC‑TGA‑CC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCg               >  1:600763/1‑151 (MQ=255)
                                                                                                                                               cA‑T‑‑GGGTGCAATGCGTC‑TGA‑CC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACAcc        >  2:801492/1‑151 (MQ=255)
                                                                                                                                                     gggTGCAATGCGTC‑TGA‑CC‑G‑AAATTGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGTCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTATAAGCGATCGACCGTCTTCGTAACACctc      <  2:100054/150‑1 (MQ=255)
                                                                                                                                                     gggTGCAATGCGTC‑TGA‑CC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACctct     <  1:734757/151‑1 (MQ=255)
                                                                                                                                                        tGCAATGCGTC‑TGA‑CC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCtcttct  >  1:459767/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                          
ACATGATCAACCGTGGCGGTACGTTCCTCGGTTCTGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTA‑CA‑T‑‑GGGTGCAATGCGTC‑TGA‑CC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCT  >  NC_000913/4107730‑4108028

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