Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,756,377 T→C L227P (CTC→CCC)  pykF → pyruvate kinase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,756,3770TC100.0% 229.9 / NA 62L227P (CTC→CCC) pykFpyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (36/26);  total (36/26)

GATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTT  >  NC_000913/1756231‑1756519
                                                                                                                                                  |                                                                                                                                              
gATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCcaa                                                                                                                                             <  1:299363/150‑1 (MQ=255)
     ttGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCcaaca                                                                                                                                           <  2:793370/147‑1 (MQ=255)
       ggTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTcg                                                                                                                                     <  1:441723/151‑1 (MQ=255)
          tGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTcgacg                                                                                                                                  >  2:1030941/1‑151 (MQ=255)
             gAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGaaa                                                                                                                               <  1:631804/151‑1 (MQ=255)
               aCAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATc                                                                                                                             <  2:1104384/151‑1 (MQ=255)
                  aGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTc                                                                                                                          <  1:353964/151‑1 (MQ=255)
                       tAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAg                                                                                                                      >  1:11288/1‑150 (MQ=255)
                       tAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGc                                                                                                                     >  1:693172/1‑151 (MQ=255)
                       tAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGc                                                                                                                     >  2:596115/1‑151 (MQ=255)
                         gACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCct                                                                                                                   >  1:1126829/1‑151 (MQ=255)
                         gACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCct                                                                                                                   >  1:511646/1‑151 (MQ=255)
                           cTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCctct                                                                                                                 >  1:1196149/1‑151 (MQ=255)
                            tttGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTg                                                                                                                >  2:498626/1‑151 (MQ=255)
                                ttgCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACgg                                                                                                            >  2:754590/1‑151 (MQ=255)
                                   ctgctTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGcat                                                                                                         >  1:37991/1‑151 (MQ=255)
                                    tgctTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGcatc                                                                                                        >  1:347636/1‑151 (MQ=255)
                                        tCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATgg                                                                                                    >  2:1090359/1‑151 (MQ=255)
                                             tATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTgcg                                                                                               >  2:90700/1‑151 (MQ=255)
                                                tCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGt                                                                                            >  1:246466/1‑151 (MQ=255)
                                                 cGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTg                                                                                           >  1:314149/1‑151 (MQ=255)
                                                   tAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGc                                                                                         >  2:177380/1‑151 (MQ=255)
                                                      gCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGAc                                                                                      >  1:1039632/1‑151 (MQ=255)
                                                       cGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCACGGTTGCGCGTGGCGAcc                                                                                     >  2:470372/1‑151 (MQ=255)
                                                        gTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCt                                                                                    >  2:187501/1‑151 (MQ=255)
                                                        gTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCt                                                                                    >  1:440102/1‑151 (MQ=255)
                                                        gTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCt                                                                                    <  2:1048386/151‑1 (MQ=255)
                                                                 ttATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGgtgt                                                                              >  2:199276/1‑148 (MQ=255)
                                                                  tATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGa                                                                           <  1:39663/150‑1 (MQ=255)
                                                                   aTCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGaaa                                                                         <  1:199274/151‑1 (MQ=255)
                                                                     cGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAAt                                                                        <  1:71472/150‑1 (MQ=255)
                                                                      gAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATcc                                                                      >  2:390002/1‑151 (MQ=255)
                                                                       aaaTCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATccc                                                                     <  1:652300/151‑1 (MQ=255)
                                                                           ccGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTa                                                                 <  1:598458/151‑1 (MQ=255)
                                                                             gTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTaga                                                               <  2:541153/151‑1 (MQ=255)
                                                                                aGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTagaaga                                                            <  1:789203/151‑1 (MQ=255)
                                                                                       gAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATc                                                     >  1:990779/1‑151 (MQ=255)
                                                                                        aaaGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCt                                                    >  1:950171/1‑151 (MQ=255)
                                                                                         aaGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCtt                                                   <  2:1163539/151‑1 (MQ=255)
                                                                                          aGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTc                                                  <  1:106719/151‑1 (MQ=255)
                                                                                             gcACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGcc                                               >  1:87691/1‑151 (MQ=255)
                                                                                               aCGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGccc                                              >  1:580074/1‑150 (MQ=255)
                                                                                                 ggcggcGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGccc                                              >  1:26862/1‑148 (MQ=255)
                                                                                                    ggcgAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCaga                                           >  2:495604/1‑148 (MQ=255)
                                                                                                      cgAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGTAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAgatg                                      <  1:341317/151‑1 (MQ=255)
                                                                                                       gAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAgatga                                     <  2:335289/151‑1 (MQ=255)
                                                                                                         aaaCATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAgatgatg                                   >  2:975497/1‑151 (MQ=255)
                                                                                                         aaaCATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAgatgatg                                   >  2:587949/1‑151 (MQ=255)
                                                                                                         aaaCATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAgatgatg                                   <  2:927864/151‑1 (MQ=255)
                                                                                                                  catcatCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGaaaaa                          <  2:440105/151‑1 (MQ=255)
                                                                                                                      atcTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGt                       >  2:764504/1‑150 (MQ=255)
                                                                                                                      atcTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGt                       >  1:456909/1‑150 (MQ=255)
                                                                                                                          ccAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTAt                     <  2:652161/148‑1 (MQ=255)
                                                                                                                           cAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGt                 <  1:1056729/151‑1 (MQ=255)
                                                                                                                              aaTCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGc               >  1:538002/1‑150 (MQ=255)
                                                                                                                              aaTCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCa              <  1:1165584/151‑1 (MQ=255)
                                                                                                                                  gAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACg            <  1:756114/149‑1 (MQ=255)
                                                                                                                                       ccAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTc     <  2:350360/151‑1 (MQ=255)
                                                                                                                                       ccAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTc     <  2:988217/151‑1 (MQ=255)
                                                                                                                                         aGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGt   >  2:1105959/1‑151 (MQ=255)
                                                                                                                                         aGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGt   >  1:1189531/1‑151 (MQ=255)
                                                                                                                                          ggAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGtt  >  2:398822/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
GATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTT  >  NC_000913/1756231‑1756519

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