Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,755,705 (A)5→4 coding (8/1413 nt) pykF → pyruvate kinase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,755,7010A.98.4% 573.6 / 1.6 125coding (4/1413 nt)pykFpyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base . (60/63);  total (61/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

Only 100 of 139 total aligned reads displayed.

CACCACCACTTTCGTAATACCGGATTCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCT  >  NC_000913/1755604‑1755795
                                                                                                 |                                                                                              
 accaccacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGaaaa                                                                                             <  1:3045232/100‑1 (MQ=255)
 accaccacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGaaaa                                                                                             >  1:202913/1‑100 (MQ=255)
  ccaccacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGaaaag                                                                                            >  1:1262718/1‑99 (MQ=255)
  ccaccacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGaaaag                                                                                            <  1:14508/100‑2 (MQ=255)
   caccacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGaaaaga                                                                                           <  2:535326/100‑3 (MQ=255)
   caccacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGaaaaga                                                                                           >  1:2821413/1‑98 (MQ=255)
    accacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGaaaagac                                                                                          >  2:426351/1‑97 (MQ=255)
     ccacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGaaaagacc                                                                                         <  2:76340/100‑5 (MQ=255)
     ccacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGaaaagacc                                                                                         >  2:3238893/1‑96 (MQ=255)
     ccacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGaaaaga                                                                                           <  2:1677916/98‑3 (MQ=255)
        cTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCaa                                                                                      >  2:712930/1‑99 (MQ=255)
          ttCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAAt                                                                                   >  1:888666/1‑100 (MQ=255)
              tAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTATAGACTAAGACTGTCATGAAAAGACCAAAATTGtt                                                                               >  2:2887150/1‑100 (MQ=255)
              tAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGtt                                                                               >  1:2210161/1‑100 (MQ=255)
              tAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGtt                                                                               >  1:2187993/1‑100 (MQ=255)
              gctcttccgATCCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGGCTAAGACTGTCATGAAAAGACCAAAATTGtt                                                                               <  2:871637/92‑1 (MQ=255)
               aaTACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGttt                                                                              >  1:318803/1‑100 (MQ=255)
               aaTACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGttt                                                                              <  2:1501957/100‑1 (MQ=255)
                aTACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTg                                                                             <  1:2591553/100‑1 (MQ=255)
                 tACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGc                                                                            >  1:1271025/1‑100 (MQ=255)
                 tACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGc                                                                            <  1:1818150/100‑1 (MQ=255)
                 tACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGc                                                                            >  2:2268789/1‑100 (MQ=255)
                  aCCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCa                                                                           <  1:426352/100‑1 (MQ=255)
                  aCCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCa                                                                           >  2:651221/1‑100 (MQ=255)
                   ccGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCAc                                                                          >  1:2898999/1‑100 (MQ=255)
                    cGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCAcc                                                                         <  2:875381/100‑1 (MQ=255)
                        ttCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCg                                                                     >  1:621193/1‑100 (MQ=255)
                        ttCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCg                                                                     <  1:3192018/100‑1 (MQ=255)
                          cGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCAGCATCGGa                                                                   >  2:332110/1‑100 (MQ=255)
                           gCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGAc                                                                  <  2:1706264/100‑1 (MQ=255)
                           gCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGAc                                                                  <  2:312231/100‑1 (MQ=255)
                            cTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGAcc                                                                 >  1:734071/1‑100 (MQ=255)
                            cTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGAcc                                                                 >  2:518215/1‑100 (MQ=255)
                            cTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGAcc                                                                 <  1:588077/100‑1 (MQ=255)
                            cTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGAcc                                                                 <  1:2268795/100‑1 (MQ=255)
                             tttCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCg                                                                >  1:3179594/1‑100 (MQ=255)
                             tttCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCg                                                                >  1:1592296/1‑100 (MQ=255)
                               tCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGaa                                                              >  2:250235/1‑100 (MQ=255)
                                ccGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGaaa                                                             >  1:3258698/1‑100 (MQ=255)
                                  gacAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCAGCATCGGACCGAAAAc                                                           <  1:332111/98‑1 (MQ=255)
                                  gacAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAAc                                                           <  1:3134161/98‑1 (MQ=255)
                                    cAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCg                                                         <  2:2025486/100‑1 (MQ=255)
                                      gTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGaa                                                       >  1:2182432/1‑100 (MQ=255)
                                       tGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGaa                                                       <  1:410618/100‑1 (MQ=255)
                                       tGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAAt                                                      <  1:380971/100‑1 (MQ=255)
                                         cgcCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCt                                                    >  2:2341301/1‑100 (MQ=255)
                                          gcCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTg                                                   >  2:3400561/1‑100 (MQ=255)
                                           cccAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGa                                                  <  2:104621/100‑1 (MQ=255)
                                              aGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAaga                                               >  2:2270751/1‑100 (MQ=255)
                                               gAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAagag                                              >  2:2773738/1‑100 (MQ=255)
                                                 aaGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGAt                                            <  2:2256468/100‑1 (MQ=255)
                                                 aaGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGAt                                            <  1:2777368/100‑1 (MQ=255)
                                                 aaGCAAGTGTCGCCCATCCGTGGCAACGTAAAGACTAGTACGGTTATGAAATGTCCAAAATTGTTTGCACCATCGTACCGAAAAGCGAATCTGAAGAGAt                                            >  2:1132209/1‑100 (MQ=255)
                                                  aGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATg                                           >  1:833739/1‑100 (MQ=255)
                                                    cAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGtt                                         <  1:1021540/100‑1 (MQ=255)
                                                     aaGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTa                                        <  1:3400574/100‑1 (MQ=255)
                                                     aaGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTa                                        <  1:2436144/100‑1 (MQ=255)
                                                       gTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGc                                      <  2:1262715/100‑1 (MQ=255)
                                                        tttCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCt                                     <  2:1204334/100‑1 (MQ=255)
                                                             cccATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAAt                                >  1:957721/1‑100 (MQ=255)
                                                             cccATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAAt                                >  2:1882525/1‑100 (MQ=255)
                                                               cATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGc                              <  1:2220623/100‑1 (MQ=255)
                                                                aTCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCt                             >  2:2879556/1‑100 (MQ=255)
                                                                aTCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCt                             >  2:1015789/1‑100 (MQ=255)
                                                                 tCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTg                            <  2:202912/100‑1 (MQ=255)
                                                                   cTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGa                          <  2:3258685/100‑1 (MQ=255)
                                                                   cTTCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGa                          <  2:2898988/100‑1 (MQ=255)
                                                                    ttCTCAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGAc                         >  1:1209033/1‑100 (MQ=255)
                                                                        cAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTg                     >  2:1751530/1‑100 (MQ=255)
                                                                        cAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTg                     >  2:1130143/1‑100 (MQ=255)
                                                                        cAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTg                     >  1:1763646/1‑100 (MQ=255)
                                                                        cAACTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTg                     <  2:621192/100‑1 (MQ=255)
                                                                         aaCTTAAAGACTAATACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTgg                    <  2:2312371/100‑1 (MQ=255)
                                                                         aaCTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTgg                    <  2:1209030/100‑1 (MQ=255)
                                                                         aaCTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTgg                    <  2:318802/100‑1 (MQ=255)
                                                                           cTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCa                  <  2:2187987/100‑1 (MQ=255)
                                                                           cTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCa                  <  2:734071/100‑1 (MQ=255)
                                                                           cTTAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCa                  <  1:1838952/100‑1 (MQ=255)
                                                                            tcatgacagtagactGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCAt                 >  2:3318596/11‑100 (MQ=255)
                                                                             tAAAGACTAAGACTGTCATGAAAAGACCAAAATTGTTTACACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCAt                 >  2:2098096/1‑99 (MQ=255)
                                                                              aaaGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGa               <  1:2887161/100‑1 (MQ=255)
                                                                              aaaGACTAAGACTGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGa               >  2:426065/1‑100 (MQ=255)
                                                                                 gactaagactGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACg            >  1:2031216/1‑100 (MQ=255)
                                                                                 gactaagactGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACg            >  1:2131155/1‑100 (MQ=255)
                                                                                 acagtagactGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACg            <  1:3318609/95‑1 (MQ=255)
                                                                                  actaagactGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACg            >  2:1446025/1‑100 (MQ=255)
                                                                                   ctaagactGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGtt          >  1:685496/1‑100 (MQ=255)
                                                                                   ctaagactGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGtt          <  1:2270757/100‑1 (MQ=255)
                                                                                     aagactGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTAt        <  2:313655/100‑1 (MQ=255)
                                                                                     aagactGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTAt        >  2:2632136/1‑100 (MQ=255)
                                                                                     aagactGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTAt        <  1:2537177/100‑1 (MQ=255)
                                                                                     aagactGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTAt        >  1:3124406/1‑100 (MQ=255)
                                                                                     aagactGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTAt        >  1:3253154/1‑100 (MQ=255)
                                                                                        actGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCg     <  1:101108/100‑1 (MQ=255)
                                                                                         ctGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGt    >  1:556366/1‑100 (MQ=255)
                                                                                          tGTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTc   <  1:1015790/100‑1 (MQ=255)
                                                                                           gTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCt  >  2:527741/1‑100 (MQ=255)
                                                                                           gTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCt  >  1:2563149/1‑100 (MQ=255)
                                                                                           gTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCt  <  2:2761283/100‑1 (MQ=255)
                                                                                           gTCATGAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCt  <  2:2631354/100‑1 (MQ=255)
                                                                                                 |                                                                                              
CACCACCACTTTCGTAATACCGGATTCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCT  >  NC_000913/1755604‑1755795

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