Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,757,006 G→A A437T (GCT→ACT)  pykF → pyruvate kinase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,757,0060GA100.0% 239.9 / NA 74A437T (GCT→ACT) pykFpyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (34/40);  total (34/40)

ATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAAT  >  NC_000913/1756875‑1757130
                                                                                                                                    |                                                                                                                            
aTGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCa                                                                                                                       >  1:1716063/1‑139 (MQ=255)
aTGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCa                                                                                                                       >  1:1683138/1‑139 (MQ=255)
  gCCACCATCCTGGCACTGGCCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCaga                                                                                                                     <  2:217933/139‑1 (MQ=255)
  gCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCaga                                                                                                                     <  2:767339/139‑1 (MQ=255)
   ccaccaTCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGGACTGAGCAATGGCGATGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGTTTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCagag                                                                                                                    <  1:1226332/139‑1 (MQ=255)
       caTCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGt                                                                                                                <  1:1859859/139‑1 (MQ=255)
       caTCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGt                                                                                                                <  2:365412/139‑1 (MQ=255)
          ccTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTg                                                                                                             <  1:1337754/139‑1 (MQ=255)
           cTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTgg                                                                                                            <  2:1683138/139‑1 (MQ=255)
            tGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGc                                                                                                           >  2:951665/1‑139 (MQ=255)
                   gACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACaaa                                                                                                    >  1:133423/1‑139 (MQ=255)
                    accaccAACGAAAGAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACaaa                                                                                                    >  1:1874814/1‑139 (MQ=255)
                    accaccAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCagag                                                                                                                    <  1:857212/122‑1 (MQ=255)
                    accaccAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAg                                                                                                   >  2:74077/1‑139 (MQ=255)
                    accaccAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAg                                                                                                   >  2:211302/1‑139 (MQ=255)
                    accaccAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAg                                                                                                   >  1:1758658/1‑139 (MQ=255)
                    accaccAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAg                                                                                                   >  1:1296736/1‑139 (MQ=255)
                    accaccAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTAATGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCagag                                                                                                                    >  2:857212/1‑122 (MQ=255)
                          aaCGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACg                                                                                             <  2:1867720/139‑1 (MQ=255)
                             gAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGATGTTg                                                                                          <  1:978697/139‑1 (MQ=255)
                                       cTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGt                                                                                <  2:1875242/139‑1 (MQ=255)
                                       cTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGt                                                                                <  1:951665/139‑1 (MQ=255)
                                               ttGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTg                                                                        <  2:1613590/139‑1 (MQ=255)
                                                       gAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTgg                                                                          <  2:648052/129‑1 (MQ=255)
                                                       gAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTgg                                                                          >  1:648052/1‑129 (MQ=255)
                                                         gCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCa                                                                                                                       >  1:1851041/1‑82 (MQ=255)
                                                         gCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCa                                                                                                                       <  2:1851041/82‑1 (MQ=255)
                                                         gCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTAcc                                                              >  2:748985/1‑139 (MQ=255)
                                                         gCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTAcc                                                              >  2:180007/1‑139 (MQ=255)
                                                            aaGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAg                                                           <  2:10490/139‑1 (MQ=255)
                                                            aaGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAg                                                           <  1:180007/139‑1 (MQ=255)
                                                                  ttGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCagag                                                                                                                    >  1:1284246/1‑76 (MQ=255)
                                                                  ttGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCagag                                                                                                                    <  2:1284246/76‑1 (MQ=255)
                                                                  ttGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCac                                                     >  1:792664/1‑139 (MQ=255)
                                                                    gtgCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGa                                                                                               >  1:1587061/1‑95 (MQ=255)
                                                                    gtgCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGa                                                                                               >  1:1370964/1‑95 (MQ=255)
                                                                    gtgCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGa                                                                                               <  2:1587061/95‑1 (MQ=255)
                                                                    gtgCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGa                                                                                               <  2:1370964/95‑1 (MQ=255)
                                                                                 ttAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCAt                                       <  2:1874814/139‑1 (MQ=255)
                                                                                 ttAAAGAGATCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATc                                      <  2:133423/139‑1 (MQ=255)
                                                                                      gagaTCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTc                                 <  2:792664/139‑1 (MQ=255)
                                                                                         aTCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACg                              <  1:74077/139‑1 (MQ=255)
                                                                                         aTCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACg                              <  2:1716063/139‑1 (MQ=255)
                                                                                           cACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACTAAGGTGACGTTGTAGTTATGGTTTCTgg                                                                          >  2:1011481/1‑93 (MQ=255)
                                                                                           cACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAATGTGACGtt                                                                                           >  1:1098038/1‑76 (MQ=255)
                                                                                           cACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGtt                                                                                           <  2:1098038/76‑1 (MQ=255)
                                                                                           cACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTgg                                                                          <  1:1011481/93‑1 (MQ=255)
                                                                                           cACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCa                                        <  1:174932/127‑1 (MQ=255)
                                                                                           cACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCa                                        >  2:174932/1‑127 (MQ=255)
                                                                                                   cTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCa                                                                      >  1:53855/1‑89 (MQ=255)
                                                                                                   cTGATGATTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCa                                                                      <  2:53855/89‑1 (MQ=255)
                                                                                                     gatgatTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACAGCATCTGTTCACGTCCTGTAATAtt                  >  2:1044636/1‑139 (MQ=255)
                                                                                                      atgatTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTg                 >  2:1723878/1‑139 (MQ=255)
                                                                                                        gatTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTc                                 >  1:817735/1‑121 (MQ=255)
                                                                                                        gatTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTc                                 <  2:817735/121‑1 (MQ=255)
                                                                                                        gatTTCTACCGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCt               >  1:1540247/1‑139 (MQ=255)
                                                                                                                ccGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGaa       <  1:211302/139‑1 (MQ=255)
                                                                                                                ccGTCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGaa       <  2:1758658/139‑1 (MQ=255)
                                                                                                                   tCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGcacac                                                                                                       <  2:824686/40‑1 (MQ=255)
                                                                                                                   tCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGcacac                                                                                                       >  1:824686/1‑40 (MQ=255)
                                                                                                                   tCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACg                                                                                             <  2:656290/50‑1 (MQ=255)
                                                                                                                   tCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACg                                                                                             >  1:656290/1‑50 (MQ=255)
                                                                                                                   tCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCa                                >  1:1460079/1‑111 (MQ=255)
                                                                                                                   tCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCa                                <  2:1460079/111‑1 (MQ=255)
                                                                                                                   tCT‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGaatta    >  1:1417434/1‑139 (MQ=255)
                                                                                                                     t‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGc                                                       >  2:63964/1‑86 (MQ=255)
                                                                                                                     t‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGc                                                       <  1:63964/86‑1 (MQ=255)
                                                                                                                     t‑GGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGaattaat  <  2:1540247/139‑1 (MQ=255)
                                                                                                                       gggTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGc                                                       <  2:1208921/85‑1 (MQ=255)
                                                                                                                       gggCAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTCCCGAGCGGc                                                       >  1:1208921/1‑85 (MQ=255)
                                                                                                                             aGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAg                                                           <  2:223228/75‑1 (MQ=255)
                                                                                                                             aGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAg                                                           >  1:223228/1‑75 (MQ=255)
                                                                                                                             aGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTaa                                               <  2:490141/87‑1 (MQ=255)
                                                                                                                             aGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTaa                                               >  1:490141/1‑87 (MQ=255)
                                                                                                                                    |                                                                                                                            
ATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCT‑GGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAAT  >  NC_000913/1756875‑1757130

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