Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,966,582 C→T 73.1% 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,5820CT73.1% 161.1 / 75.6 110G618D (GGC→GAC) ptsPPEP‑protein phosphotransferase enzyme I; GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base C (11/18);  new base T (35/45);  total (47/63)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.98e-01

Only 100 of 112 total aligned reads displayed.

CTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCG  >  NC_000913/2966487‑2966680
                                                                                               |                                                                                                  
cTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAAc                                                                                                <  1:2850679/100‑1 (MQ=255)
   gTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACaga                                                                                             <  1:2498160/100‑1 (MQ=255)
    tAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACagag                                                                                            >  2:103017/1‑100 (MQ=255)
    tAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTCAGTAAGATCGTTGGTGCCAACagag                                                                                            <  1:2228443/100‑1 (MQ=255)
     aaCATTGCATGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAgaga                                                                                           <  2:2451297/100‑1 (MQ=255)
     aaCATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAgaga                                                                                           <  2:728803/100‑1 (MQ=255)
        gatctaGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGa                                                                                        <  2:251570/95‑1 (MQ=255)
        aTTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGa                                                                                        >  1:53965/1‑100 (MQ=255)
           gCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAAt                                                                                     <  1:1324688/100‑1 (MQ=255)
             aGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGaagatc                                                                                   >  1:251568/1‑96 (MQ=255)
              ggATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGa                                                                                  >  2:184438/1‑100 (MQ=255)
               gATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGAc                                                                                 <  2:1771975/100‑1 (MQ=255)
               gATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGAc                                                                                 <  1:2134635/100‑1 (MQ=255)
               gATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGAc                                                                                 >  2:1978183/1‑100 (MQ=255)
                  gAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCg                                                                              >  1:172619/1‑100 (MQ=255)
                     gACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCtt                                                                           >  2:2355937/1‑100 (MQ=255)
                       cTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTg                                                                         <  1:213691/100‑1 (MQ=255)
                        tGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGc                                                                        <  1:2825129/100‑1 (MQ=255)
                        tGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGc                                                                        <  1:1124454/100‑1 (MQ=255)
                         gTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGcc                                                                       <  1:917383/100‑1 (MQ=255)
                          tCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCa                                                                      <  1:267255/100‑1 (MQ=255)
                          tCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCa                                                                      <  2:471418/100‑1 (MQ=255)
                           cATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAg                                                                     >  2:1714464/1‑100 (MQ=255)
                            aTAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGa                                                                    >  1:730131/1‑100 (MQ=255)
                            aTAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGa                                                                    >  2:2813241/1‑100 (MQ=255)
                              aaaTGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCTGCTTTGCCAGacg                                                                  <  2:1973339/100‑3 (MQ=255)
                              aaaTGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATg                                                                  >  1:917995/1‑100 (MQ=255)
                               aaTGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGc                                                                 <  2:931602/100‑1 (MQ=255)
                                aTGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCg                                                                <  2:595107/100‑1 (MQ=255)
                                aTGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCg                                                                <  2:909596/100‑1 (MQ=255)
                                aTGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCg                                                                <  1:1710270/100‑1 (MQ=255)
                                   ttCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgca                                                             <  2:1037968/100‑1 (MQ=255)
                                    tCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgcag                                                            <  2:1350341/100‑1 (MQ=255)
                                    tCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgcag                                                            <  2:1950790/100‑1 (MQ=255)
                                      gCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgcagca                                                          <  2:33852/100‑1 (MQ=255)
                                      gCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgcagca                                                          <  1:1148344/100‑1 (MQ=255)
                                      gCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgcagca                                                          >  1:1462211/1‑100 (MQ=255)
                                      gCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgcagca                                                          >  2:1857866/1‑100 (MQ=255)
                                      gCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgcagca                                                          >  2:210118/1‑100 (MQ=255)
                                       ccACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATTAAATCGAGCCGCTTTGCCAGATGCGGCAGCAt                                                         <  1:1674603/100‑1 (MQ=255)
                                       ccACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCAt                                                         >  2:2521019/1‑100 (MQ=255)
                                        cACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATa                                                        <  2:2663555/100‑1 (MQ=255)
                                        cACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATa                                                        <  1:86296/100‑1 (MQ=255)
                                        cACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATa                                                        >  1:147408/1‑100 (MQ=255)
                                          cccGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATaaa                                                      >  2:1201072/1‑100 (MQ=255)
                                            cGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATa                                                    <  2:1865295/100‑1 (MQ=255)
                                            cGGGGATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCTGCCGCAtctggc                                                    <  1:591296/100‑6 (MQ=255)
                                              ggTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATAcc                                                  <  1:1885656/100‑1 (MQ=255)
                                              ggTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATAcc                                                  <  1:1201070/100‑1 (MQ=255)
                                              ggTATTGTTGCGACCAACGGCCAGAATGGATTGACTCAGATCGTTAGTGTCCACAGAGCCGACATCGACCCGCTTTGCCAGATGCGGCCGCATCAATAcc                                                  <  1:1884990/100‑1 (MQ=255)
                                                  ttgttgCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTg                                              >  2:958249/1‑100 (MQ=255)
                                                    gttgCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGAc                                            <  2:2623286/100‑1 (MQ=255)
                                                    gttgCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGAc                                            >  1:482332/1‑100 (MQ=255)
                                                     ttgCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACg                                           <  2:120327/100‑1 (MQ=255)
                                                      tgCGATCAACGTCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACgg                                          <  2:190849/100‑1 (MQ=255)
                                                       gCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGc                                         <  2:2660323/100‑1 (MQ=255)
                                                       gCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGc                                         <  1:1104601/100‑1 (MQ=255)
                                                       gCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATAGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGc                                         >  1:1866581/1‑100 (MQ=255)
                                                        cGATCAGCGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGc                                         <  1:826677/99‑1 (MQ=255)
                                                         gATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCAc                                       >  1:1536363/1‑100 (MQ=255)
                                                         gATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGAACCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCAc                                       >  1:1398563/1‑100 (MQ=255)
                                                          aTCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACt                                      >  1:788231/1‑100 (MQ=255)
                                                          aTCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGACAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACt                                      >  1:1550677/1‑100 (MQ=255)
                                                            cAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTc                                    >  2:209693/1‑100 (MQ=255)
                                                            cAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTc                                    <  1:213082/100‑1 (MQ=255)
                                                             aaCGGCCAGGATGTATTGCGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAACACCATTGACGGCACTTcc                                   <  2:1781234/100‑1 (MQ=255)
                                                               cGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAg                                 <  2:318302/100‑1 (MQ=255)
                                                               cGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAg                                 >  2:1441575/1‑100 (MQ=255)
                                                                ggCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGc                                >  2:1900146/1‑100 (MQ=255)
                                                                 gCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCa                               >  2:184509/1‑100 (MQ=255)
                                                                 gCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCa                               >  2:226281/1‑100 (MQ=255)
                                                                  ccAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCat                              >  1:1939005/1‑100 (MQ=255)
                                                                        tGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGcc                        <  1:1445256/100‑1 (MQ=255)
                                                                         gTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCg                       <  1:209691/100‑1 (MQ=255)
                                                                          tATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCg                       <  2:54821/99‑1 (MQ=255)
                                                                          tATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGa                      >  2:1169643/1‑100 (MQ=255)
                                                                          tATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACACGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGa                      >  1:1585145/1‑100 (MQ=255)
                                                                           aTTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGAt                     >  2:746547/1‑100 (MQ=255)
                                                                           aTTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGAt                     <  2:2858953/100‑1 (MQ=255)
                                                                           aTTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGAt                     <  2:611972/100‑1 (MQ=255)
                                                                             tGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATAc                   >  1:1085502/1‑100 (MQ=255)
                                                                              gAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACg                  >  2:307546/1‑100 (MQ=255)
                                                                               aGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGt                 >  1:1316901/1‑100 (MQ=255)
                                                                               aGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGt                 >  1:2624063/1‑100 (MQ=255)
                                                                                gTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTg                <  2:650662/100‑1 (MQ=255)
                                                                                  cAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGt              >  1:2428185/1‑100 (MQ=255)
                                                                                   aGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGt              >  2:1194439/1‑99 (MQ=255)
                                                                                    gATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACAATTGACGGCACTTCCAGCATGATGCCGATACGTGGttt            >  1:2465312/1‑100 (MQ=255)
                                                                                      tCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTgg          <  1:1515212/100‑1 (MQ=255)
                                                                                      tCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGGTTg           >  1:2117777/1‑99 (MQ=255)
                                                                                      tCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTgg          <  2:1832779/100‑1 (MQ=255)
                                                                                       cGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTggg         <  1:1698574/100‑1 (MQ=255)
                                                                                       cGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTggg         <  2:1380756/100‑1 (MQ=255)
                                                                                       cGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTggg         <  2:178378/100‑1 (MQ=255)
                                                                                       cGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGACGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTggg         <  1:1019189/100‑1 (MQ=255)
                                                                                         ttGGTGTCAACAGAGCTGACATCGACCCGCTTTGCCAGATGCGGCAGCAGAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGaa       <  2:584284/100‑1 (MQ=255)
                                                                                         ttGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGaa       <  2:245689/100‑1 (MQ=255)
                                                                                           ggTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAAtt     <  2:786876/100‑1 (MQ=255)
                                                                                            gtgtCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAAttt    >  2:1746044/1‑100 (MQ=255)
                                                                                            gtgtCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAAttt    <  1:1501805/100‑1 (MQ=255)
                                                                                               |                                                                                                  
CTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCG  >  NC_000913/2966487‑2966680

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: