Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,966,582 C→T 100% 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,5820CT99.3% 458.8 / ‑4.6 142G618D (GGC→GAC) ptsPPEP‑protein phosphotransferase enzyme I; GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (71/69);  total (72/70)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.17e-01

Only 100 of 142 total aligned reads displayed.

AGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATT  >  NC_000913/2966485‑2966677
                                                                                                 |                                                                                               
agCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACTCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCa                                                                                               <  1:3524582/100‑1 (MQ=255)
   tggttaCATGGCAAGATGAGGACTGTCATAAATGTTTGCCAGCAGGGGGTTTTTGCGATCAACGGTAAGAATGTATTGATGCACATAGTTGATGTCaat                                                                                             <  2:419036/94‑2 (MQ=255)
     gTAACATTGCAGGATGAAGACTGTCATAAATGTTTGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAAc                                                                                             <  1:900577/97‑1 (MQ=255)
     gTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACaga                                                                                          <  1:1462469/100‑1 (MQ=255)
       aaCATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAgaga                                                                                        <  2:2761725/100‑1 (MQ=255)
        aCATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGAt                                                                                       >  2:464882/1‑100 (MQ=255)
         cATTGAAGGATTAAGACTGTCATAAATGTTGGTAATCCGGATATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATTGTTGGTGTTAACAGAGATg                                                                                      <  2:1735499/100‑1 (MQ=255)
          aTTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGTGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATg                                                                                      <  1:2593496/99‑1 (MQ=255)
          aTTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGACCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGa                                                                                     >  1:2083842/1‑100 (MQ=255)
           ttGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTAAGATCGTTGGTGTCAACAGAGATGaa                                                                                    >  1:3784258/1‑100 (MQ=255)
            tGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGaaa                                                                                   <  1:757313/100‑1 (MQ=255)
             gCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAAt                                                                                  <  2:2209065/100‑1 (MQ=255)
              cAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATc                                                                                 <  1:3315776/100‑1 (MQ=255)
              cAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATc                                                                                 >  1:1588880/1‑100 (MQ=255)
               aGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCg                                                                                >  2:2497165/1‑100 (MQ=255)
                  aTGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGAcc                                                                             >  1:2019246/1‑100 (MQ=255)
                  aTGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGAcc                                                                             <  2:1145765/100‑1 (MQ=255)
                    gAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCg                                                                           <  2:1666357/100‑1 (MQ=255)
                    gAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCATTGGTGTCAACAGAGATGAAATCGACCCg                                                                           >  1:3776324/1‑100 (MQ=255)
                     aaGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGCTCAACGTCAAGAATGTATTGAATCAAATCATTGTAGTCAATAAAGATGAAATCGACACg                                                                           >  1:50621/1‑99 (MQ=255)
                      aGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCt                                                                         >  1:1200045/1‑100 (MQ=255)
                       gACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCtt                                                                        >  2:82165/1‑100 (MQ=255)
                         cTGTCATAAATGTTCGCCACCCGTTTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTg                                                                      <  1:1242033/100‑1 (MQ=255)
                         cTGTCATAAATGTTCGCCACCCGGGTATTGTTGCTATTAATGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTg                                                                      <  1:323622/100‑1 (MQ=255)
                         cTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTg                                                                      >  2:3208380/1‑100 (MQ=255)
                         cTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAAAAGAGATGAAATCGACCCGCTTTg                                                                      >  1:461190/1‑100 (MQ=255)
                         cTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGACAACAGAGATGAAATCGACCCGCTTTg                                                                      >  1:683063/1‑100 (MQ=255)
                           gTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGcc                                                                    >  1:3817635/1‑100 (MQ=255)
                           gTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCAACCCGCTTTGcc                                                                    >  1:127394/1‑100 (MQ=255)
                           gTCATAAATGTTCGCCACCAGGGGATTATTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGcc                                                                    >  2:2426497/1‑100 (MQ=255)
                             cATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAg                                                                  >  1:3763351/1‑100 (MQ=255)
                             cATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAg                                                                  >  1:2204515/1‑100 (MQ=255)
                              aTAAATGTTCGCCACCCGGGTATTGTTGCGATCCACGGCCAGAATGTATTGAATCAGATCGTTGGTGTCAACAAAGATGAAATCGACCCGCTTTGCCAGa                                                                 >  1:2948299/1‑100 (MQ=255)
                              aTAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGa                                                                 >  2:3060093/1‑100 (MQ=255)
                                  aTGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGc                                                              <  1:3843920/99‑1 (MQ=255)
                                   tGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGc                                                              >  2:895634/1‑98 (MQ=255)
                                    gTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCgg                                                            >  2:1932437/1‑99 (MQ=255)
                                    gTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgc                                                           <  2:467594/100‑1 (MQ=255)
                                     ttCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgca                                                          >  2:125966/1‑100 (MQ=255)
                                       ctccACCCGGGTATTTTTGCGATCAACTTCCAGAATGTATTGATTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgcagc                                                        <  1:2275681/98‑1 (MQ=255)
                                        gCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgcagca                                                       <  1:2933055/100‑1 (MQ=255)
                                        gCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgcagca                                                       <  2:3499447/100‑1 (MQ=255)
                                        gCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGCGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATACGgca                                                          >  2:2552669/1‑97 (MQ=255)
                                         ccACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGCTGAAATCGACCCGCTTTGCCAGATGCGgcagca                                                       >  2:3036370/1‑99 (MQ=255)
                                          catcccggTATTGTTGCTATCAACAGCCAGAATGTATTGAGTAAGATCGTTGGTGTCAAAAGAGATGAAATCGACCCGCTTTGTCAGATGCGGCAGCAt                                                      >  2:3199149/7‑99 (MQ=255)
                                          cACCCGGGTATTGTTGCGATCAACGGCCAGAATGTCTTGAGTCAGATCATTAGTGTCAACAGAGAAAAAAACGACCCGCTTTGCCAGATGCGGCAGCATa                                                     >  1:2009747/1‑100 (MQ=255)
                                             ccGGGTATTGTTTCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAAt                                                  <  1:664799/100‑1 (MQ=255)
                                              cGGGTATTTTTGCGATCAACGGCCAGATTGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGTTAGATGCGGCAGCATAAAt                                                  <  2:2182951/99‑1 (MQ=255)
                                              cGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATa                                                 <  2:1384249/100‑1 (MQ=255)
                                              cGGGTATTGTTGCGATAAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATa                                                 <  2:734630/100‑1 (MQ=255)
                                               gggTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATAc                                                <  2:1937060/100‑1 (MQ=255)
                                                ggTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATAcc                                               <  1:2749161/100‑1 (MQ=255)
                                                ggTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATAcc                                               <  2:1106484/100‑1 (MQ=255)
                                                  tATTGTTGCGATCAACGGCCAGAATGTATTGAGGCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCAt                                             <  2:3854922/100‑1 (MQ=255)
                                                       ttgCGATCAACGGCCAGAATGTATTGATTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACg                                        <  1:2497188/100‑1 (MQ=255)
                                                           gATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCa                                     >  1:1928505/1‑99 (MQ=255)
                                                             tCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACtt                                  >  1:3472672/1‑100 (MQ=255)
                                                               aaCGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTcc                                >  1:933778/1‑100 (MQ=255)
                                                                aCGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCa                               >  1:3044518/1‑100 (MQ=255)
                                                                  ggCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGc                             >  1:143921/1‑100 (MQ=255)
                                                                  ggCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGc                             <  1:1164016/100‑1 (MQ=255)
                                                                   gCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCa                            <  2:3045674/100‑1 (MQ=255)
                                                                   gCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCa                            >  1:1729189/1‑100 (MQ=255)
                                                                     cAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCa                            <  1:1100078/98‑1 (MQ=255)
                                                                         aTGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGc                      <  1:2726402/100‑1 (MQ=255)
                                                                          tGTATTGAGTCAGATCGTTGGTGTCAACAGTGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTTACGGCACTTCCAGCATGATGcg                     <  2:18337/100‑2 (MQ=255)
                                                                          tGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGcn                     <  2:3307542/100‑2 (MQ=255)
                                                                          tGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGcc                     <  1:311815/100‑1 (MQ=255)
                                                                           gTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCg                    <  1:3208411/100‑1 (MQ=255)
                                                                            tATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGCTGCGGCAGCATAAATACCATTCACGGCACTTCCAGCatcatga                      >  2:3846897/1‑96 (MQ=255)
                                                                            tATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTACCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGa                   >  1:3073648/1‑100 (MQ=255)
                                                                            tATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCCACCCGCCTTGCCACATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGa                   <  2:2312186/100‑1 (MQ=255)
                                                                             ttccgatctgATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGAt                  <  1:2942993/91‑1 (MQ=255)
                                                                              ttGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGATCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATa                 <  1:3161476/100‑1 (MQ=255)
                                                                              ttGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGTACTTCCAGCATGATGCCGATa                 >  1:3444794/1‑100 (MQ=255)
                                                                              ttGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGTACTTCCAGCATGATGCCGATa                 >  1:3444812/1‑100 (MQ=255)
                                                                              ttGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATa                 <  2:683055/100‑1 (MQ=255)
                                                                               tGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATAc                >  1:3469892/1‑100 (MQ=255)
                                                                               tGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCAGATa                 >  2:3223190/1‑99 (MQ=255)
                                                                                gAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACg               >  1:3581314/1‑100 (MQ=255)
                                                                                gAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACg               >  2:3789707/1‑100 (MQ=255)
                                                                                gAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACg               >  1:3923209/1‑100 (MQ=255)
                                                                                  gTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTg             >  2:1573968/1‑100 (MQ=255)
                                                                                    cAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGt           >  2:3515769/1‑100 (MQ=255)
                                                                                     aGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGtt          >  1:3572739/1‑100 (MQ=255)
                                                                                     aGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGtt          >  2:1082635/1‑100 (MQ=255)
                                                                                     aGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGt           <  2:3469853/99‑1 (MQ=255)
                                                                                     aGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCACATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGtt          <  2:111043/100‑1 (MQ=255)
                                                                                      gATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATagatcggaa         >  2:2942967/1‑92 (MQ=255)
                                                                                      gATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGttt         >  2:3575204/1‑100 (MQ=255)
                                                                                      gATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGACACTTCCAGCATGATGCCGATACGTGGttt         >  1:555086/1‑100 (MQ=255)
                                                                                       aTCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTACAGCATGATGCCGATACGTGGtttn        <  2:2019228/100‑2 (MQ=255)
                                                                                        tCGTTGGTGTCAACAGAGATGAAATCGACTCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTgg       <  1:1848389/100‑1 (MQ=255)
                                                                                        tCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTTCCAGATGCGGCAGCATAAATACCATTGACGGCACTTTCAGCATGATGCCGATACGTGGTTTgg       <  1:2386872/100‑1 (MQ=255)
                                                                                        tCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTgg       <  2:338941/100‑1 (MQ=255)
                                                                                        tCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTgg       <  1:2389214/100‑1 (MQ=255)
                                                                                          gTTGGGGTCAACAGAGATGAAATGGACCCGCTTTGCCAGATGCGGCAGCATAAATAGCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTGTGGGa     <  2:2778308/100‑1 (MQ=255)
                                                                                            tGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAAt   >  1:1887310/1‑100 (MQ=255)
                                                                                            tGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAAt   >  1:1510668/1‑100 (MQ=255)
                                                                                             ggTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAAtt  >  2:2306052/1‑100 (MQ=255)
                                                                                                 |                                                                                               
AGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATT  >  NC_000913/2966485‑2966677

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