Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,966,582 C→T G618D (GGC→GAC)  ptsP ← PEP‑protein phosphotransferase enzyme I; GAF domain containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,966,5820CT98.2% 351.2 / ‑2.6 108G618D (GGC→GAC) ptsPPEP‑protein phosphotransferase enzyme I; GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (60/46);  total (61/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.20e-01

Only 100 of 109 total aligned reads displayed.

GCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTT  >  NC_000913/2966486‑2966678
                                                                                                |                                                                                                
gCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCaa                                                                                               >  2:1841465/1‑100 (MQ=255)
 cTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAAc                                                                                              <  2:1338840/100‑1 (MQ=255)
 cTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAAc                                                                                              <  2:2738696/100‑1 (MQ=255)
  tCGTAACATTGTAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACa                                                                                             <  1:1769771/100‑1 (MQ=255)
  tCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACa                                                                                             >  1:3022259/1‑100 (MQ=255)
  tCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGACCAGAATGTATTGAGTCAGATCGTTGGTGTCAACa                                                                                             >  2:580183/1‑100 (MQ=255)
     tAACATTGCAGGATGAAGACTGTCATAAACGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACagag                                                                                          <  2:1968820/100‑1 (MQ=255)
      aaCATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAgaga                                                                                         <  2:2715927/100‑1 (MQ=255)
       aCATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGAt                                                                                        >  1:2459121/1‑100 (MQ=255)
       aCATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGAt                                                                                        >  2:1796106/1‑100 (MQ=255)
        cATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATg                                                                                       >  2:638897/1‑100 (MQ=255)
        cATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATg                                                                                       <  1:1660766/100‑1 (MQ=255)
        cATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCACAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATg                                                                                       <  2:1739822/100‑1 (MQ=255)
         aTTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGa                                                                                      >  1:1361622/1‑100 (MQ=255)
           tGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGaaa                                                                                    <  2:3099271/100‑1 (MQ=255)
             cAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAAt                                                                                   >  2:1165840/1‑99 (MQ=255)
             cAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATc                                                                                  <  2:127396/100‑1 (MQ=255)
             cAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATc                                                                                  <  1:2461976/100‑1 (MQ=255)
              aGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGTGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCg                                                                                 >  1:2283377/1‑100 (MQ=255)
              aGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCg                                                                                 <  1:2847752/100‑1 (MQ=255)
                gATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGAc                                                                               >  2:673159/1‑100 (MQ=255)
                 aTGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGAcc                                                                              >  1:3107400/1‑100 (MQ=255)
                  tGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGAccc                                                                             >  1:2454892/1‑100 (MQ=255)
                  tGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGAccc                                                                             >  1:2959501/1‑100 (MQ=255)
                   gAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCg                                                                            >  2:364856/1‑100 (MQ=255)
                   gAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCg                                                                            >  1:236466/1‑100 (MQ=255)
                    aaGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCg                                                                            >  2:891168/1‑99 (MQ=255)
                    aaGACTGTCATAAATGTTCACCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGc                                                                           >  1:2239517/1‑100 (MQ=255)
                     aGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGCTGGTGTCAACAGAGATGAAATCGACCCGCt                                                                          <  1:462689/100‑1 (MQ=255)
                         tGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCtttcc                                                                      <  1:2212794/100‑3 (MQ=255)
                          gTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGcc                                                                     >  1:2233434/1‑100 (MQ=255)
                            cATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAg                                                                   >  1:2436/1‑100 (MQ=255)
                            cATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAg                                                                   >  2:582332/1‑100 (MQ=255)
                            cATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAg                                                                   >  2:383978/1‑100 (MQ=255)
                              tAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGAt                                                                 <  1:2535536/100‑1 (MQ=255)
                              tAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGAt                                                                 >  1:39226/1‑100 (MQ=255)
                               aaaTGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATg                                                                <  1:647527/100‑1 (MQ=255)
                               aaaTGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATg                                                                >  2:1699207/1‑100 (MQ=255)
                               aaaTGTTCGCCACCCGGGTATTGTTGCGATAAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGAACCGCTTTGCCAGATg                                                                >  2:1477285/1‑100 (MQ=255)
                               aaaTGTTCGCCACCCGGGTATTGTTGCGAACAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTGAACAGAGATGAAATCGACCCGCTTTGCCAGATg                                                                <  2:1967484/100‑1 (MQ=255)
                                aaTGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGc                                                               <  2:1874872/100‑1 (MQ=255)
                                aaTGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGc                                                               <  2:290060/100‑1 (MQ=255)
                                aaTGTTCGCCACCCGGGTATTGATGCGATCAGCGGCCTGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGc                                                               <  2:2203527/100‑1 (MQ=255)
                                 aTGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGcn                                                              <  2:1691934/100‑2 (MQ=255)
                                 aTGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCg                                                              <  1:1525952/100‑1 (MQ=255)
                                   gTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgc                                                            >  2:1473213/1‑100 (MQ=255)
                                   gTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTACCAGATGCGgc                                                            >  1:1822108/1‑100 (MQ=255)
                                     tCGCCACCCGGGTATTGTTGTGATCAACGGCCAGAATGTATTCAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgcag                                                          <  2:2283351/100‑1 (MQ=255)
                                     tCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgcag                                                          <  1:1414528/100‑1 (MQ=255)
                                      cGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgcagc                                                         <  1:86833/100‑1 (MQ=255)
                                       gCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgcagca                                                        <  2:245320/100‑1 (MQ=255)
                                       gCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgcagca                                                        >  2:2900447/1‑100 (MQ=255)
                                       gCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgcagca                                                        >  1:2663498/1‑100 (MQ=255)
                                       gCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgcagca                                                        >  2:1690199/1‑100 (MQ=255)
                                       gCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgcagca                                                        >  2:1329288/1‑100 (MQ=255)
                                        ccACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCAt                                                       >  1:1535326/1‑100 (MQ=255)
                                         cACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATa                                                      >  1:1585452/1‑100 (MQ=255)
                                           cccGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATaaa                                                    >  2:1873624/1‑100 (MQ=255)
                                             cGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATa                                                  <  1:231229/100‑1 (MQ=255)
                                                gTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCa                                               <  1:1841487/100‑1 (MQ=255)
                                                 tATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCAt                                              <  2:1586923/100‑1 (MQ=255)
                                                  aTTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCAtt                                             <  2:2959466/100‑1 (MQ=255)
                                                   ttgttgCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTg                                            >  2:2505612/1‑100 (MQ=255)
                                                     gttgCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGAc                                          >  2:2943012/1‑100 (MQ=255)
                                                       tgCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACgg                                        >  2:2063209/1‑100 (MQ=255)
                                                       tgCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACgg                                        >  1:766877/1‑100 (MQ=255)
                                                         cGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCa                                      <  1:1189249/100‑1 (MQ=255)
                                                          gATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCAtgacacgggc                                     >  2:39319/1‑91 (MQ=255)
                                                          gATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCAc                                     >  2:2824707/1‑100 (MQ=255)
                                                                cGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAg                               >  1:1690762/1‑100 (MQ=255)
                                                                cGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGACAGCATAAATACCATTGACGGCACTTCCAg                               >  1:2850875/1‑100 (MQ=255)
                                                                 ggCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGc                              <  2:2893242/100‑1 (MQ=255)
                                                                  gCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGc                              >  2:2979599/1‑99 (MQ=255)
                                                                  gCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCa                             >  1:2163849/1‑100 (MQ=255)
                                                                    cAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCatg                           >  1:1382982/1‑100 (MQ=255)
                                                                     aGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCCTTGACGGCACTTCCAgatag                           >  2:1845524/1‑95 (MQ=255)
                                                                     aGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCatga                          >  2:2521166/1‑100 (MQ=255)
                                                                     aGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCatga                          >  2:1616625/1‑100 (MQ=255)
                                                                     aGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCatga                          >  2:2348206/1‑100 (MQ=255)
                                                                       aaTGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCatgatg                        >  2:2504798/1‑100 (MQ=255)
                                                                       aaTGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCatgatg                        >  2:1648463/1‑100 (MQ=255)
                                                                       aaTGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCatgatg                        >  2:2504614/1‑100 (MQ=255)
                                                                         tGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGcn                      <  1:2996873/100‑2 (MQ=255)
                                                                          gTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCg                     <  2:1026208/100‑1 (MQ=255)
                                                                           tATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGa                    >  2:1730126/1‑100 (MQ=255)
                                                                               gAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACg                >  2:782158/1‑100 (MQ=255)
                                                                                  tCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTgg             >  2:854572/1‑100 (MQ=255)
                                                                                    aGATCGTTGGTGTCAACAGAGATGACATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGtt           <  2:2438403/100‑1 (MQ=255)
                                                                                     gATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGttt          >  2:1153344/1‑100 (MQ=255)
                                                                                      aTCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTg         <  2:3022223/100‑1 (MQ=255)
                                                                                      aTCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTg         <  2:1160270/100‑1 (MQ=255)
                                                                                      aTCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTg         <  2:373399/100‑1 (MQ=255)
                                                                                       tCGTTGGTGTCAACAGAGATGAATTCTACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTg         <  1:580192/99‑1 (MQ=255)
                                                                                       tCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTgg        <  1:344426/100‑1 (MQ=255)
                                                                                       tCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTgg        <  2:236460/100‑1 (MQ=255)
                                                                                        cGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTggg       <  2:2984348/100‑1 (MQ=255)
                                                                                         gTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGa      >  1:2546984/1‑100 (MQ=255)
                                                                                         gTTGGTGTCAACAGAGATGAAATCGACCCGCTTTACCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGa      <  2:1822086/100‑1 (MQ=255)
                                                                                            ggTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAAtt   >  1:850629/1‑100 (MQ=255)
                                                                                             gtgtCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAAttt  <  2:148325/100‑1 (MQ=255)
                                                                                                |                                                                                                
GCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTT  >  NC_000913/2966486‑2966678

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: