Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,966,582 C→T 27.2% 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,5820CT27.2% 156.7 / 65.8 96G618D (GGC→GAC) ptsPPEP‑protein phosphotransferase enzyme I; GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base C (32/37);  new base T (15/11);  total (48/48)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 101 total aligned reads displayed.

AGAGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGT  >  NC_000913/2966483‑2966681
                                                                                                   |                                                                                                   
agagCTCGTAACATTGCAGGATTAAGACTGTCATACATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGc                                                                                                     >  2:1366135/1‑100 (MQ=255)
 gagCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGGTCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGcc                                                                                                    <  2:1903455/100‑1 (MQ=255)
 gagCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGcc                                                                                                    >  1:2237505/1‑100 (MQ=255)
  agCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCa                                                                                                   <  2:1925708/100‑1 (MQ=255)
    cTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAAc                                                                                                 <  2:712181/100‑1 (MQ=255)
    cTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAAc                                                                                                 <  2:1676959/100‑1 (MQ=255)
     tCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAAcn                                                                                                <  2:1791709/100‑2 (MQ=255)
       gTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACaga                                                                                              >  2:703667/1‑100 (MQ=255)
             ttGCAGGATGAATACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGaa                                                                                        <  1:686535/100‑1 (MQ=255)
             ttGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGaa                                                                                        <  1:574830/100‑1 (MQ=255)
             tggCAGGACGAAGGCTGTCATCGATGTTCGCCACCCGTGTATTGTTGCGATCAACGGCCAGAATGGAGTGAGTCAGATCGTTGGTGTCAACAGAGATg                                                                                          <  2:531957/96‑1 (MQ=255)
              tGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGaaa                                                                                       >  2:2189685/1‑100 (MQ=255)
               gCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAAt                                                                                      <  1:1328004/100‑1 (MQ=255)
                cAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATc                                                                                     >  1:499319/1‑100 (MQ=255)
                 aGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCg                                                                                    <  2:764497/100‑1 (MQ=255)
                     tGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGAccc                                                                                <  1:2192472/100‑1 (MQ=255)
                      gAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGAccc                                                                                >  2:1769814/1‑99 (MQ=255)
                      gAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCg                                                                               >  2:215398/1‑100 (MQ=255)
                       aaGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGc                                                                              >  2:1578926/1‑100 (MQ=255)
                        acaCTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCt                                                                             <  1:242437/98‑1 (MQ=255)
                        aGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCt                                                                             <  1:176400/100‑1 (MQ=255)
                             gTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAGCGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGcn                                                                        <  2:2212652/100‑2 (MQ=255)
                             gTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGcc                                                                        <  2:2248117/100‑1 (MQ=255)
                             gTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGcc                                                                        <  2:2315308/100‑1 (MQ=255)
                               cATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAg                                                                      >  1:1893107/1‑100 (MQ=255)
                               cATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAg                                                                      >  1:1310842/1‑100 (MQ=255)
                               cATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAg                                                                      <  2:1589296/100‑1 (MQ=255)
                                 tAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGAt                                                                    >  2:386577/1‑100 (MQ=255)
                                 tAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGACAGCATGTATTGAGTCAGATCGTTGGTGCCACCAGAGATGACATCGACCCGCTTTGCCAGAt                                                                    >  2:1260265/1‑100 (MQ=255)
                                    aTGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCg                                                                 <  2:2161101/100‑1 (MQ=255)
                                     tGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCgg                                                                >  2:1561589/1‑100 (MQ=255)
                                      gtccGCCACCCGTGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgc                                                               <  1:326472/97‑1 (MQ=255)
                                      gTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCTAGATGCGgc                                                               <  1:2170838/100‑1 (MQ=255)
                                      gTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCgga                                                               <  2:321262/100‑2 (MQ=255)
                                       ttCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgca                                                              <  1:1830039/100‑1 (MQ=255)
                                       ttCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgca                                                              <  2:1752177/100‑1 (MQ=255)
                                        tCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgcag                                                             >  1:1799808/1‑100 (MQ=255)
                                        tCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgcag                                                             >  1:839565/1‑100 (MQ=255)
                                         cGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgcagc                                                            >  2:1723541/1‑100 (MQ=255)
                                          gCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgcagca                                                           >  2:1883541/1‑100 (MQ=255)
                                          gCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgcagca                                                           >  1:1506018/1‑100 (MQ=255)
                                          gCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgcagca                                                           <  2:2237486/100‑1 (MQ=255)
                                               ccGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAAt                                                      <  1:710018/100‑1 (MQ=255)
                                               ccGGGTATTGTGGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAAt                                                      <  1:1457415/100‑1 (MQ=255)
                                                cGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATa                                                     <  1:475412/100‑1 (MQ=255)
                                                 gggTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATAc                                                    >  2:588422/1‑100 (MQ=255)
                                                  ggTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATAcc                                                   <  1:95186/100‑1 (MQ=255)
                                                      ttgttgCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTg                                               <  1:1110188/100‑1 (MQ=255)
                                                       tgttgCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGa                                              >  2:1520801/1‑100 (MQ=255)
                                                        gttgCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGAc                                             >  1:848564/1‑100 (MQ=255)
                                                        gttgCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGAc                                             >  1:1626475/1‑100 (MQ=255)
                                                        gttgCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGAc                                             >  1:651112/1‑100 (MQ=255)
                                                        gttgCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGAc                                             >  1:454911/1‑100 (MQ=255)
                                                        gttgCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGAc                                             >  1:1976831/1‑100 (MQ=255)
                                                         ttgCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACg                                            <  1:1055833/100‑1 (MQ=255)
                                                            cGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGc                                          >  1:713185/1‑99 (MQ=255)
                                                            cGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCa                                         >  1:680289/1‑100 (MQ=255)
                                                            cGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCa                                         <  1:1635861/100‑1 (MQ=255)
                                                            cGATCAACGGCCAGAATGTATTGAGGCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCa                                         <  1:1077279/100‑1 (MQ=255)
                                                             gATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCAc                                        <  1:739412/100‑1 (MQ=255)
                                                               tCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACtt                                      <  2:893467/100‑1 (MQ=255)
                                                                cAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTc                                     <  1:1822254/100‑1 (MQ=255)
                                                                cAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTc                                     <  1:994979/100‑1 (MQ=255)
                                                                  aCGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCa                                   >  1:2352578/1‑100 (MQ=255)
                                                                    ggCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGc                                 >  1:543095/1‑100 (MQ=255)
                                                                     gCCAGAATGTCTTGAGACCGATCGTCGGTGTAAACAGAGATGAACTCGACCAGCTTTGCCAGATGCGGCAGCATAAAGACGGTTGAGGGCACTGccgtc                                 <  2:2201200/99‑4 (MQ=255)
                                                                     gCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCa                                >  2:931644/1‑100 (MQ=255)
                                                                      ccAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCat                               >  2:2109604/1‑100 (MQ=255)
                                                                      ccAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCat                               >  1:771358/1‑100 (MQ=255)
                                                                      ccAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCat                               >  1:2052863/1‑100 (MQ=255)
                                                                      ccAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCat                               >  1:2052511/1‑100 (MQ=255)
                                                                       cAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCatg                              <  2:238756/100‑1 (MQ=255)
                                                                       cAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCat                               <  2:1309239/99‑1 (MQ=255)
                                                                        aGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCatga                             <  1:1127203/100‑1 (MQ=255)
                                                                        aGAATGTATTGAGTCAGATCGTTGGTGCCAACATAGATGAAATCGACCCGTTTTGCCAGATGCGGCAGCATAAATACCATTTACGGCACTTCCAGCat                               <  2:583996/98‑1 (MQ=255)
                                                                                ttGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATa                     <  2:431495/100‑1 (MQ=255)
                                                                                 tGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATAc                    >  1:852088/1‑100 (MQ=255)
                                                                                 tGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATAc                    >  2:1721889/1‑100 (MQ=255)
                                                                                  gAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACg                   <  1:703670/100‑1 (MQ=255)
                                                                                  gAGTCAGATCGTTGGTGAAAAAAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATaag                   >  2:1758225/1‑98 (MQ=255)
                                                                                   aGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGt                  >  1:853057/1‑100 (MQ=255)
                                                                                   aGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGt                  >  1:1178218/1‑100 (MQ=255)
                                                                                    gTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTg                 >  1:691652/1‑100 (MQ=255)
                                                                                     tCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTgg                <  1:801180/100‑1 (MQ=255)
                                                                                     tCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTgg                >  1:2271986/1‑100 (MQ=255)
                                                                                      cAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGt               >  1:1098541/1‑100 (MQ=255)
                                                                                       aGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGtt              >  2:1570051/1‑100 (MQ=255)
                                                                                        gATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGttt             >  1:336591/1‑100 (MQ=255)
                                                                                        gATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGttt             >  2:77884/1‑100 (MQ=255)
                                                                                        gATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGttt             >  2:154440/1‑100 (MQ=255)
                                                                                            gTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGa         >  1:91763/1‑100 (MQ=255)
                                                                                            gTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGa         >  1:1454970/1‑100 (MQ=255)
                                                                                            gTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGa         >  2:1745589/1‑100 (MQ=255)
                                                                                               ggTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAAtt      >  2:279413/1‑100 (MQ=255)
                                                                                               ggTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAAtt      <  2:835490/100‑1 (MQ=255)
                                                                                                gTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAAttt     <  2:1178213/100‑1 (MQ=255)
                                                                                                 tGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTc    <  1:723504/100‑1 (MQ=255)
                                                                                                  gCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCg   <  1:1745601/100‑1 (MQ=255)
                                                                                                   ccAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGt  >  1:2029961/1‑100 (MQ=255)
                                                                                                   ccAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGt  >  1:1707754/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
AGAGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGT  >  NC_000913/2966483‑2966681

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