Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,755,8620CA41.7% 163.4 / 63.0 97S55R (AGC→AGA) pykFpyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base C (41/15);  new base A (0/40);  total (41/56)
Fisher's exact test for biased strand distribution p-value = 7.11e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 101 total aligned reads displayed.

AATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATG  >  NC_000913/1755712‑1756012
                                                                                                                                                        |                                                                                                                                                       
aaTTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAGTGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGc                                                                                                                                                         <  3:11390/151‑1 (MQ=255)
aaTTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGc                                                                                                                                                         >  1:216657/1‑151 (MQ=255)
       tGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACt                                                                                                                                                   <  3:362638/150‑1 (MQ=255)
          aCCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTa                                                                                                                                              <  1:1498763/151‑1 (MQ=255)
          aCCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTa                                                                                                                                              <  4:1361912/151‑1 (MQ=255)
            cATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTaaa                                                                                                                                            <  2:119875/151‑1 (MQ=255)
            cATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTaa                                                                                                                                             >  4:587329/1‑150 (MQ=255)
             aTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTaaa                                                                                                                                            >  2:501284/1‑150 (MQ=255)
               cGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAAcc                                                                                                                                         >  1:703416/1‑151 (MQ=255)
               cGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAAcc                                                                                                                                         <  2:421155/151‑1 (MQ=255)
               cGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAAc                                                                                                                                          <  1:530417/150‑1 (MQ=255)
                ggACCGAAAACCGAATCTGAAGATATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAAccg                                                                                                                                        <  4:1235742/151‑1 (MQ=255)
                ggACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAAccg                                                                                                                                        <  1:1756270/151‑1 (MQ=255)
                ggACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAAccg                                                                                                                                        >  1:1771679/1‑151 (MQ=255)
                ggACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAAccg                                                                                                                                        <  4:777464/151‑1 (MQ=255)
                   ccGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAAccgcc                                                                                                                                      <  4:226216/150‑1 (MQ=255)
                     gAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAACCGCCGCt                                                                                                                                   >  1:1178745/1‑151 (MQ=255)
                     gAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAACCGCCGCt                                                                                                                                   >  1:1139371/1‑151 (MQ=255)
                             aaTCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCt                                                                                                                           >  1:305221/1‑151 (MQ=255)
                             aaTCTAAAGAGTTTTTAGCAAAAATGCGGGCCGCTGGCATAATCTTTATGCGTCTGACCTT‑CTCTCATGGTGATTATGCAAACCGGGGTCAGCCAATTGAAATCCGGCGAAGC‑GTGAGGAGAAAAGCG‑GGAAAAACCGCCGCTATCCTGCt                                                                                                                           <  2:429760/151‑1 (MQ=255)
                                   aaGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATAc                                                                                                                     >  2:46703/1‑151 (MQ=255)
                                      agaTGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCa                                                                                                                   >  4:95436/1‑150 (MQ=255)
                                        aTGTTAGCAAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTAACCTT‑CTCTCAGGGTAACTATGCAAAACACGGCCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGAAAAACG‑GGAAAAACCGCCGCTATCCTGCTTGATACCa                                                                                                                   <  1:185848/148‑1 (MQ=255)
                                         ttttagaaaaagTGCGGGCCGCTGCAATGACCTTTAGGCTTCT‑AACTTCCCCTCGGGGTGACTAGGAAAAACAGGGTCAGCGCATTAAAAATCTGCGAAAC‑GTAAGGACAAAAACT‑GGAAAAACCGCCGCTATCCTGCTTGATACCaaa                                                                                                                 <  4:850703/141‑1 (MQ=255)
                                         tGTTAGCTAAAATGCTGGCCGCTGGCATGAACGTTATGCGTCTAACCTT‑CTCTCATGGTGACTATGCAAAACCCGGTCAGCGCATTCAGAATCGGCGCAAG‑GTGAGGAGAAAAAGT‑GAAAAAACCGCCGCTATCCTGCTTGATACCaa                                                                                                                  <  4:226273/148‑1 (MQ=255)
                                             aGCAAAAATGCGGGACGGTGGCATGAACTTTTGGCGTCTGACCTT‑CTCTCAGGGTGACTAGGAAAAACCCGGTCAGCGCATCAAAAACCTGCGCAGC‑GTGAAGAGAAAAACT‑GAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTc                                                                                                             <  1:1179834/149‑1 (MQ=255)
                                                tAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGa                                                                                                         >  1:175206/1‑150 (MQ=255)
                                                tAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGa                                                                                                         >  4:874594/1‑150 (MQ=255)
                                                 aaagTGCTGGCCGCGGACATAACCTTTATCCTTCTGAACTT‑CTCTAATGGTGCCTATGCAGAACACGGCCGCCGCATTAAGAATCGGCCAAAG‑GTGATGAGCAAAGGT‑GGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGa                                                                                                         <  1:1039806/149‑1 (MQ=255)
                                                 aaaaTGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGaa                                                                                                        >  4:983211/1‑150 (MQ=255)
                                                 aaaaTCCTGGACGCGGGCATAACCGTTGTGCGTCTAAACTC‑CTCTCATGGTGATTAGGCAAACCCCGGTCAGCGCTTCCAGATCCGGCGAAAG‑GTGATGAACAAACGT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGaa                                                                                                        <  1:379648/150‑1 (MQ=255)
                                                   aaTGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATc                                                                                                     >  1:1846816/1‑151 (MQ=255)
                                                   aaTGCGGGACGCTGGCATAACCTTTATGCGTCAGACCTT‑CTCTCATGGTGACTATGAAAAACCCGGTCAGCGTTTTAAAATTCTGCGAAAG‑GTGAGGAAAAAAACT‑GAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATc                                                                                                     <  1:1457258/151‑1 (MQ=255)
                                                       cTGGACGCTGGCATGAACGTTATGCGTCTGAACTT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGc                                                                                                  >  1:1686748/1‑150 (MQ=255)
                                                        tGGACGCTGGCATGAACGTTATGCGTCTGAACTT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCAc                                                                                                >  4:4060/1‑151 (MQ=255)
                                                        tGGACGCTGGCATGAACGTTATGCGTCTGAACTT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCAc                                                                                                >  2:536810/1‑151 (MQ=255)
                                                         ggccgcgGGCATAACCTTTATGCTTCTAACCTT‑CTCTCTTGGTGCCTATGCAAACCCCGGCCACCGCTTTCAGATTCTGCAAACC‑GTGATGAGAAAAACT‑GTAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATcc                                                                                                    <  2:147216/143‑1 (MQ=255)
                                                                 gCATGAACGTTATGCGTCTGAACTT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACt                                                                                       >  4:1594209/1‑151 (MQ=255)
                                                                  cATGAACGTTATGCGTCTGAACTT‑TTCTCATGTTGACTACGCAGAACACGGTAAGCGCATTAAGAATCGCCACAAC‑GTGATGAGCAAAACG‑ATTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACt                                                                                       <  4:1428037/150‑1 (MQ=255)
                                                                      aaCGTTATGCGTCTGAACTT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGa                                                                                    >  4:828772/1‑149 (MQ=255)
                                                                       actttATGCGTCTAACTTT‑CTCTCATGGTGACTATGAAGAACCCGGTCAGCGCATTCAAAATCGGCGAAAC‑GTGATGAGAAAAATT‑GGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAg                                                                                  <  4:544160/147‑1 (MQ=255)
                                                                          tttggCGTCTCACTTT‑CTCTCAGGGTCAATATGCAAAACACGGCCCGCCAATACAAAATCTGCACAAC‑GTAAAGAGAAAAACT‑GGAAAAACCGCCGCTAGCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCg                                                                               <  1:883168/146‑1 (MQ=255)
                                                                          ttATGCTTCTGATCTT‑CTCTAGGGGTGACTATGCAAAACACGGTCAGCGCTTTAAAATTCTGCCCAAC‑GGGATGAGAAAAATG‑GAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTgg                                                                                     <  1:323040/144‑1 (MQ=255)
                                                                          ttATGCGTCTAAACCT‑CGCTAAGGGTGATGATAAAGAACGCGGCCAGCTCATAAATATTCTGCGAAAG‑TTGATGAGAAAAAGT‑GGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGaa                                                                                   <  3:518335/146‑1 (MQ=255)
                                                                          ttAGGCGTCTAACCTT‑CTCTCAGGGTGACTATGCAGAACCCGGTCAGCGCATTAAAAACCGGC‑CAACGGTGAGGAAAAAAACT‑GGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGc                                                                                <  4:789321/149‑1 (MQ=255)
                                                                            aTGCGTCTGAACTT‑CTCTCATGGCGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCgg                                                                              >  4:551955/1‑149 (MQ=255)
                                                                              gCGTCTGAACTT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAac                                                                          >  4:567146/1‑151 (MQ=255)
                                                                               cttctaacttt‑CTCTTATGGTGACTAGGAAAAACACGGCCAGCGCATTCAGAATCGGCGCAAC‑GGGATGA‑AAAAACTGGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAac                                                                          <  4:1282781/141‑1 (MQ=255)
                                                                                gTCTGAACTT‑CTCTCATGGTGACTATGCAGACCACGGTCAGCTCATTCAGAATCTGCGAAAC‑GTGATGAGAAAATCT‑GAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTa                                                                            <  4:1236546/147‑1 (MQ=255)
                                                                                  ctaaccTT‑CTCTCAGGGTGACTCCGCAAAACACGGACAGCTCTTCCATAATCTGCGCAAG‑GTGATGAGAAAAACT‑GATAAACCCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTa                                                                            <  3:349546/140‑1 (MQ=255)
                                                                                    gAACTT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGtt                                                                    >  4:1506309/1‑151 (MQ=255)
                                                                                    gAACTT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGtt                                                                    >  2:269710/1‑151 (MQ=255)
                                                                                    gAACTT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGtt                                                                    >  4:1144443/1‑151 (MQ=255)
                                                                                     acttt‑CTCTCTGGGTGATTATGAAGAACACGGTCACCTCATACAAAATCTGCAAAAC‑GAAAAGAGAAAAACG‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGt                                                                     <  4:1822704/146‑1 (MQ=255)
                                                                                     accTT‑CTCTCATGGTGACTATGCAGAACACGGCCAGCGCATTCAAAATCTGCGAAAG‑GTGAGGAGAAAAACT‑GAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGtt                                                                    <  4:528933/148‑1 (MQ=255)
                                                                                        tt‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTtctc                                                                >  4:1171594/1‑151 (MQ=255)
                                                                                            tctctgggTGACTAAGAAAAACACGGCCAGCCAATTAAAAATCTGCAAAAC‑GTAAGGAGAAAATGT‑GAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTtct                                                                 <  1:1142437/141‑1 (MQ=255)
                                                                                            tctcATGGGGACTATGCAAAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTtct                                                                 <  2:29523/147‑1 (MQ=255)
                                                                                               cATGGTGATTAGGAAAAACAGGGTAAGCGAATTCAAAATCTGCGCAAC‑GTGAGGAGAAAAACG‑GAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGaaa                                                           <  2:501429/150‑1 (MQ=255)
                                                                                               cATGGTGACTATGCATAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGaa                                                            >  4:678531/1‑149 (MQ=255)
                                                                                                aTGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGc                                                         >  4:1654385/1‑151 (MQ=255)
                                                                                                 tGGTGACTATGCAGAACCCGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGAGGAGAAAACCT‑GGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGc                                                         <  1:152573/150‑1 (MQ=255)
                                                                                                  ggTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACC‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTg                                                       >  4:304305/1‑151 (MQ=255)
                                                                                                  ggTGACTATGCAGAACACGGCCAGCGCATTAAAAATCTGCCAAAC‑GTGATGAAAAAAACT‑GAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGaaa                                                           <  4:1847536/147‑1 (MQ=255)
                                                                                                  ggTGACTAGGCAGAACACGGTCGGCGCATTAAGAATCTGCGCAAC‑GTGAGGAGAAAAAGT‑GAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAg                                                          <  4:837427/148‑1 (MQ=255)
                                                                                                  ggTGACTAGGCAAAACCGGGCCAGCGCTTTAAGAATCGGCGAAAC‑GTGAGGAGAAAAACG‑GGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGaaa                                                           <  4:1588291/147‑1 (MQ=255)
                                                                                                  ggTGACTAGCAAAAACACGGTCAGCACATTCAAAATCTGCGCAAC‑GTGATGAGAAAAACT‑GAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTg                                                       <  2:349621/151‑1 (MQ=255)
                                                                                                  ggTGAATATAAAAAACGCGGTCAGCGCATACAGACTCTGCCCACC‑GTGATGAAAAAAACT‑GATAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTg                                                       <  4:731613/151‑1 (MQ=255)
                                                                                                         atcaaaaccacGGTCAGCGCATCCAAAATCTCCGCAAC‑GTGATGAAAAAAACT‑GAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAg                                                  <  2:55426/141‑1 (MQ=255)
                                                                                                         agcaaaaaacaggTCAGCGCATTCAGAATCTGCGCAAT‑GGTAAGAATAAATCT‑AGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAg                                                  <  3:516555/138‑1 (MQ=255)
                                                                                                         aagaaaaccacGGTCAGCCAATTAAGAACCGGCGCAAG‑GTGAAGAAAAAAGCT‑GGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAg                                                  <  1:1725914/141‑1 (MQ=255)
                                                                                                         aTGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGAGGAGAAAAACT‑GATAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGa                                                 <  4:960756/150‑1 (MQ=255)
                                                                                                         aTGCAAAACAGGGCCAGCGCATCCAGATTCTGCACAAC‑GTGATGAAAAAAACT‑GAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAg                                                  <  1:453668/149‑1 (MQ=255)
                                                                                                             aGAACCTGGCCAGCCAATTCAGATTGCGCCCAAT‑GTAATGACAAAAAGT‑GAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCt                                              <  4:237526/149‑1 (MQ=255)
                                                                                                             aGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCttt                                            >  4:666119/1‑151 (MQ=255)
                                                                                                               aaCCCGGCCAGCGCATTAAGAATCTGCCAAAC‑GTGATGAGAAAAACT‑GGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTAc                                          <  3:194729/151‑1 (MQ=255)
                                                                                                               aaCAGGGTGAGCGCATTAAGAATCCGCGCAAC‑GTGGAGAAAAAAGCG‑GAAAAACCCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCttt                                            <  1:1432197/149‑1 (MQ=255)
                                                                                                               aaCACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTAc                                          >  2:477210/1‑151 (MQ=255)
                                                                                                                 cacGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACtt                                        >  4:1625962/1‑151 (MQ=255)
                                                                                                                    ggccggcgcATTAAAAACCGGCGCAAG‑GGAAGGAAAAAAACG‑GAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTca                                     <  1:1004743/146‑1 (MQ=255)
                                                                                                                          cgcATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTg                               >  4:738120/1‑151 (MQ=255)
                                                                                                                           gcATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTg                               >  4:1613548/1‑150 (MQ=255)
                                                                                                                            cATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGAt                             >  1:69498/1‑151 (MQ=255)
                                                                                                                            cATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGAt                             >  1:1342796/1‑151 (MQ=255)
                                                                                                                             aTTCAAAATCTGCGCAAC‑GTGATGAGAAAAACT‑GGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACt                                <  4:1601733/147‑1 (MQ=255)
                                                                                                                                cAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAAt                         >  1:1119123/1‑151 (MQ=255)
                                                                                                                                 agatccggcgaAAC‑GTGATAAGAAAACCG‑GAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATc                        <  2:299638/144‑1 (MQ=255)
                                                                                                                                 aaaaccggcacAAC‑GAGAAGAGAAAAACG‑GGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATaaa                          <  1:304198/139‑1 (MQ=255)
                                                                                                                                 aaaaTCTGCGCAAC‑GTGATGAGAAAAACT‑GGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATaaa                          <  4:402144/147‑1 (MQ=255)
                                                                                                                                      cTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGtt                    >  4:1728853/1‑150 (MQ=255)
                                                                                                                                       tGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTAt                  >  4:854467/1‑151 (MQ=255)
                                                                                                                                        ggctaacttttataaGCAAAATT‑GATAATACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGGCCTTTACTTTCACCACTGATAAATCTGTTAt                  <  4:1287955/137‑1 (MQ=255)
                                                                                                                                         cgcAAC‑GTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCg                >  4:245282/1‑151 (MQ=255)
                                                                                                                                            aaC‑TTGAGGAGAAAAACG‑GTTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCa             <  4:1885336/151‑1 (MQ=255)
                                                                                                                                             aaggTGATGAGAAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGc              <  4:358440/147‑1 (MQ=255)
                                                                                                                                              c‑GTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAAc           >  1:1856763/1‑151 (MQ=255)
                                                                                                                                                 tgatgaGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAAc           >  1:480302/1‑149 (MQ=255)
                                                                                                                                                 tgatgaGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAAc           >  1:842528/1‑149 (MQ=255)
                                                                                                                                                        cAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATg  >  4:874039/1‑151 (MQ=255)
                                                                                                                                                        cAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATg  >  3:429477/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                       
AATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAC‑GTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATG  >  NC_000913/1755712‑1756012

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