Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_tpiA 1,756,356 A→C K220T (AAA→ACA)  pykF → pyruvate kinase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA1,756,3560AC98.7% 263.7 / NA 75K220T (AAA→ACA) pykFpyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (39/35);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAT  >  NC_000913_3_pae_tpiA/1756209‑1756496
                                                                                                                                                   |                                                                                                                                            
cTGAAAAAGACTAACAGGACCTGCTCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCacaa                                                                                                                                            <  2:94803/150‑1 (MQ=255)
  gAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATc                                                                                                                                          <  4:479229/150‑1 (MQ=255)
     aaaGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGaa                                                                                                                                       <  4:17612/150‑1 (MQ=255)
       aGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAAACCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGaaaa                                                                                                                                     <  2:682286/150‑1 (MQ=255)
        gACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAAc                                                                                                                                    <  2:228692/150‑1 (MQ=255)
          cAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCa                                                                                                                                  <  3:433178/150‑1 (MQ=255)
          cAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCa                                                                                                                                  <  2:592912/150‑1 (MQ=255)
           aaaCAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAg                                                                                                                                 >  1:837866/1‑150 (MQ=255)
                ggACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATAAAAAAAAAAAAAAACCAGGAAgg                                                                                                                            <  4:355296/150‑1 (MQ=255)
                 gACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAGCATCCACATCATCTTCACAATCGAAAACCAg                                                                                                                                 >  1:501254/1‑144 (MQ=255)
                 gACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTTCACAATCGAAAACCAg                                                                                                                                 <  2:501254/144‑1 (MQ=255)
                 gACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGc                                                                                                                           >  4:328766/1‑150 (MQ=255)
                    cTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTc                                                                                                                        <  1:418978/150‑1 (MQ=255)
                            tGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACtt                                                                                                                <  1:616531/150‑1 (MQ=255)
                                tGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTcgac                                                                                                            <  1:531651/150‑1 (MQ=255)
                                       aaGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATcc                                                                                                     <  4:594386/150‑1 (MQ=255)
                                         ggCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTc                                                                                                   <  3:83767/150‑1 (MQ=255)
                                          gCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAgg                                                                                                                            <  3:32906/124‑1 (MQ=255)
                                          gCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAgg                                                                                                                            >  4:32906/1‑124 (MQ=255)
                                            gTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGaa                                                                                                >  3:418717/1‑150 (MQ=255)
                                             tAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTc                                                                                                                        >  4:196169/1‑125 (MQ=255)
                                             tAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTc                                                                                                                        <  3:196169/125‑1 (MQ=255)
                                                         ctgctTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGca                                                                                   <  4:78206/150‑1 (MQ=255)
                                                            ctTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGcatca                                                                                >  4:243734/1‑150 (MQ=255)
                                                              tCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATg                                                                              >  2:479484/1‑150 (MQ=255)
                                                                cTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCCCGAAGCCTCTGACGGCATCATGGt                                                                            >  4:368376/1‑150 (MQ=255)
                                                                   tATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGcc                                                                                                                          <  3:246036/101‑1 (MQ=255)
                                                                   tATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGcc                                                                                                                          >  4:246036/1‑101 (MQ=255)
                                                                      tCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGtcgccac                                                                      >  1:661664/1‑144 (MQ=255)
                                                                         tAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTgg                                                                   >  2:837461/1‑150 (MQ=255)
                                                                          aaGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGc                                                                  >  3:355136/1‑150 (MQ=255)
                                                                            gCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTgg                                                                   >  2:239336/1‑147 (MQ=255)
                                                                            gCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTgg                                                                   <  1:239336/147‑1 (MQ=255)
                                                                            gCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTg                                                                    <  3:84111/146‑1 (MQ=255)
                                                                            gCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTg                                                                    >  4:84111/1‑146 (MQ=255)
                                                                            gCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGa                                                                >  1:80633/1‑150 (MQ=255)
                                                                            gCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGa                                                                >  1:343283/1‑150 (MQ=255)
                                                                                tCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTg                                                            <  4:628/150‑1 (MQ=255)
                                                                                 cTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTgg                                                           <  4:493370/150‑1 (MQ=255)
                                                                                  tGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTggg                                                          >  2:769405/1‑150 (MQ=255)
                                                                                  tGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTggg                                                          >  2:34299/1‑150 (MQ=255)
                                                                                   gACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGgt                                                         >  2:601767/1‑150 (MQ=255)
                                                                                    aCGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGt                                                                     >  4:365058/1‑137 (MQ=255)
                                                                                    aCGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGt                                                                     <  3:365058/137‑1 (MQ=255)
                                                                                    aCGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGgtg                                                        <  1:302137/150‑1 (MQ=255)
                                                                                      gTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGaa                                                                                                                              >  4:24568/1‑78 (MQ=255)
                                                                                      gTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGaa                                                                                                                              <  3:24568/78‑1 (MQ=255)
                                                                                        tATCGACATCCGTGAGCACCGGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAACTCCTCGAAGCCTCTGACGGCATCATGGTTGCTCGTGCCGACCTGGGTGCAGa                                                    >  4:551879/1‑150 (MQ=255)
                                                                                        tATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGa                                                    >  3:440040/1‑150 (MQ=255)
                                                                                         aTCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGaa                                                   >  2:702450/1‑150 (MQ=255)
                                                                                          tCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGaaa                                                  >  3:512330/1‑150 (MQ=255)
                                                                                           cGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAAt                                                 >  2:206911/1‑150 (MQ=255)
                                                                                            gAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATc                                                >  1:559920/1‑150 (MQ=255)
                                                                                             aaaTCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATcc                                               >  1:428734/1‑150 (MQ=255)
                                                                                                 ccGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGt                                           <  2:428734/150‑1 (MQ=255)
                                                                                                        cACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAgaagaa                                    <  4:355136/150‑1 (MQ=255)
                                                                                                         aCCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTcga                                                                                                             >  2:493509/1‑76 (MQ=255)
                                                                                                         aCCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTcga                                                                                                             <  1:493509/76‑1 (MQ=255)
                                                                                                         aCCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTa                                                      >  1:920452/1‑131 (MQ=255)
                                                                                                         aCCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTa                                                      <  2:920452/131‑1 (MQ=255)
                                                                                                         aCCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGaagaag                                   >  1:114268/1‑150 (MQ=255)
                                                                                                           cTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAGTCGAAAACCTGGAAGTCCTCAACAACTTCGACGAAATCCGCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGtt                                 >  3:391547/1‑150 (MQ=255)
                                                                                                           cTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGt                                           >  2:120915/1‑140 (MQ=255)
                                                                                                           cTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGt                                           <  1:120915/140‑1 (MQ=255)
                                                                                                           cTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGtt                                 <  1:212340/150‑1 (MQ=255)
                                                                                                            tGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTa                                >  1:187256/1‑150 (MQ=255)
                                                                                                            tGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTa                                <  3:243734/150‑1 (MQ=255)
                                                                                                                aGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACgg                                                                                     <  4:362396/93‑1 (MQ=255)
                                                                                                                aGCGCACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACgg                                                                                     >  3:362396/1‑93 (MQ=255)
                                                                                                                 gcgcACGGCGGCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTc                           <  3:600690/150‑1 (MQ=255)
                                                                                                                      cggcggCGAAAACATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTGGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCa                      >  1:570972/1‑150 (MQ=255)
                                                                                                                              aaaaCATCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAgatga              <  1:702450/150‑1 (MQ=255)
                                                                                                                                   aTCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGACGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATc         >  1:836714/1‑150 (MQ=255)
                                                                                                                                   aTCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATc         <  1:837461/150‑1 (MQ=255)
                                                                                                                                    tCCACATCATCTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCATCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCg        <  3:328766/150‑1 (MQ=255)
                                                                                                                                          tcatcTCCACAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAt  >  2:965594/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
CTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAT  >  NC_000913_3_pae_tpiA/1756209‑1756496

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: