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,5820CT97.5% 503.6 / ‑0.3 157G618D (GGC→GAC) ptsPPEP‑protein phosphotransferase enzyme I; GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (83/70);  total (85/72)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 159 total aligned reads displayed.

AGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTT  >  NC_000913/2966485‑2966678
                                                                                                 |                                                                                                
agCCCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCa                                                                                                <  1:288463/100‑1 (MQ=255)
   tCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGACTGTATTGAGTCAGATCGTTGGTGTCAACa                                                                                             >  2:951457/1‑100 (MQ=255)
   tCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCAATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACa                                                                                             <  2:3818014/100‑1 (MQ=255)
     gTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACaga                                                                                           >  1:330989/1‑100 (MQ=255)
     gTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACaga                                                                                           >  1:1888901/1‑100 (MQ=255)
      tAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACagag                                                                                          >  2:1505397/1‑100 (MQ=255)
      tAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACaga                                                                                           >  1:1802472/1‑99 (MQ=255)
      tAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACagag                                                                                          >  2:478648/1‑100 (MQ=255)
       aaCATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAgaga                                                                                         <  2:1947057/100‑1 (MQ=255)
       aaCATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAgaga                                                                                         >  2:1264539/1‑100 (MQ=255)
       aaCATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAgaga                                                                                         <  2:1758440/100‑1 (MQ=255)
        aCATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGAt                                                                                        >  1:746681/1‑100 (MQ=255)
        aCATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGAt                                                                                        >  2:1656002/1‑100 (MQ=255)
         cATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATg                                                                                       >  2:1441546/1‑100 (MQ=255)
         cATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATg                                                                                       >  1:3365734/1‑100 (MQ=255)
         cATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATg                                                                                       >  2:761797/1‑100 (MQ=255)
           ttGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGaa                                                                                     >  2:2800125/1‑100 (MQ=255)
              cAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATc                                                                                  <  1:1312283/100‑1 (MQ=255)
               aGGATGAAGACTGGCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCg                                                                                 <  2:1502832/100‑1 (MQ=255)
                ggATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGa                                                                                >  1:1914884/1‑100 (MQ=255)
                 gATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGAc                                                                               <  2:519525/100‑1 (MQ=255)
                      aGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCt                                                                          >  2:374208/1‑100 (MQ=255)
                      aGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGACCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCt                                                                          <  1:1995689/100‑1 (MQ=255)
                         cgatcttAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTAAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTg                                                                       <  1:660612/94‑1 (MQ=255)
                         cTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTg                                                                       >  2:1191319/1‑100 (MQ=255)
                         cTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTg                                                                       <  1:2075565/100‑1 (MQ=255)
                           gTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGcn                                                                     <  1:1013349/100‑2 (MQ=255)
                           gTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGcc                                                                     >  1:4126719/1‑100 (MQ=255)
                            tCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCa                                                                    >  1:2358749/1‑100 (MQ=255)
                             cATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAg                                                                   <  2:3329685/100‑1 (MQ=255)
                             cATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAg                                                                   >  2:3889111/1‑100 (MQ=255)
                             cATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAg                                                                   >  1:1145529/1‑100 (MQ=255)
                               tAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGAt                                                                 <  1:1809095/100‑1 (MQ=255)
                               tAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTAAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTgagatcg                                                                 >  2:660609/1‑94 (MQ=255)
                                aaaTGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATg                                                                >  2:811854/1‑100 (MQ=255)
                                 aaTGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGc                                                               >  1:196969/1‑100 (MQ=255)
                                  aTGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCg                                                              <  2:196966/100‑1 (MQ=255)
                                    gTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgc                                                            >  1:604109/1‑100 (MQ=255)
                                    gTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgc                                                            >  1:600365/1‑100 (MQ=255)
                                     ttCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAGATCGACCCGCTTTGCCAGATGCGgca                                                           >  1:1753003/1‑100 (MQ=255)
                                      tagCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgcag                                                          <  2:2529770/98‑1 (MQ=255)
                                       cGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgcagc                                                         >  2:3899227/1‑100 (MQ=255)
                                       cGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgcagc                                                         >  2:1332250/1‑100 (MQ=255)
                                       cGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgcagc                                                         >  2:2974917/1‑100 (MQ=255)
                                        gCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgcagca                                                        >  2:3034016/1‑100 (MQ=255)
                                          cACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATa                                                      >  1:2499434/1‑100 (MQ=255)
                                          cACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGATGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATa                                                      <  1:3729710/100‑1 (MQ=255)
                                             ccGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAAt                                                   <  2:1731163/100‑1 (MQ=255)
                                             ccGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAAt                                                   <  1:731874/100‑1 (MQ=255)
                                             ccGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAAt                                                   <  2:8577/100‑1 (MQ=255)
                                             ccGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAAt                                                   <  1:3167017/100‑1 (MQ=255)
                                               gggTATTGTTGCGATCATCGGCCAGAATGGATTGACGCAGATCGTTGGTGTCAGGAGAGATGAGATCGTACCGCCTGGGCAGCTGCGGCAGCATAAAt                                                   <  2:2812958/98‑1 (MQ=255)
                                               gggTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATAc                                                 <  2:1818765/100‑1 (MQ=255)
                                                 gTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCa                                               <  1:883036/100‑1 (MQ=255)
                                                   aTTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCAtt                                             >  2:3152925/1‑100 (MQ=255)
                                                   aTTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCAtt                                             >  1:997696/1‑100 (MQ=255)
                                                   aTTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCAtt                                             <  1:3056492/100‑1 (MQ=255)
                                                      gttgCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGa                                           >  2:4051561/1‑99 (MQ=255)
                                                      gttgCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGAc                                          >  2:3990458/1‑100 (MQ=255)
                                                       ttgCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACg                                         <  2:2344717/100‑1 (MQ=255)
                                                        tgCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACgg                                        >  1:303698/1‑100 (MQ=255)
                                                        tgCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACgg                                        >  2:1599455/1‑100 (MQ=255)
                                                        tgCGACCACCGGCCAGAATGTATTGAGGCAGATCGTAGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACg                                         <  1:2998607/99‑1 (MQ=255)
                                                         gCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGc                                       >  2:3315204/1‑100 (MQ=255)
                                                         gCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGc                                       <  2:375221/100‑1 (MQ=255)
                                                          cGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCa                                      >  1:385607/1‑100 (MQ=255)
                                                             tCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACtt                                   <  2:110294/100‑1 (MQ=255)
                                                               aaCGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTcc                                 <  1:3490864/100‑1 (MQ=255)
                                                               aaCGGCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTcc                                 >  2:1063348/1‑100 (MQ=255)
                                                                 cGGCCATAATGTATTGAGTCAGATCGTTGGTGTCAACAGATATGAAATCGACCCGCTTTGCCATATGCGGCAGCATCAATACCATTGACGGCACTTCCAg                               <  1:3541600/100‑1 (MQ=255)
                                                                   gCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCa                             >  1:490229/1‑100 (MQ=255)
                                                                   gCCAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCa                             <  1:1656014/100‑1 (MQ=255)
                                                                    ccAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCat                            <  1:2585918/100‑1 (MQ=255)
                                                                    ccAGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCat                            >  1:3098588/1‑100 (MQ=255)
                                                                      aGAATGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCatga                          >  2:3201911/1‑100 (MQ=255)
                                                                        aaTGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCatgat                         <  2:3082544/99‑1 (MQ=255)
                                                                          tGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCCCTTCCAGCATGATGcc                      >  1:3662704/1‑100 (MQ=255)
                                                                          tGTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGcc                      <  1:1449579/100‑1 (MQ=255)
                                                                           gTATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGAGGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCg                     >  1:54831/1‑100 (MQ=255)
                                                                            tATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGa                    >  1:314387/1‑100 (MQ=255)
                                                                            tATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGa                    >  1:2002688/1‑100 (MQ=255)
                                                                            tATTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGa                    >  1:8059/1‑100 (MQ=255)
                                                                             aTTGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGAt                   >  1:3456677/1‑100 (MQ=255)
                                                                              ttGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATa                  <  2:385603/100‑1 (MQ=255)
                                                                              ttGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATa                  <  1:94786/100‑1 (MQ=255)
                                                                              gtGAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGTCACTTCCAGCATGATGCCGATa                  <  2:2425371/99‑1 (MQ=255)
                                                                               tGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATAc                 >  1:3941880/1‑100 (MQ=255)
                                                                                gAGTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACg                >  2:1075123/1‑100 (MQ=255)
                                                                                  ggcAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTg              <  1:3599071/98‑1 (MQ=255)
                                                                                  gTCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGTATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTg              >  2:2861958/1‑100 (MQ=255)
                                                                                   tCAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTgg             <  1:501108/100‑1 (MQ=255)
                                                                                    cAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGt            >  2:1649064/1‑100 (MQ=255)
                                                                                    cAGATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGt            >  1:2486822/1‑100 (MQ=255)
                                                                                      gATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGttt          >  2:416796/1‑100 (MQ=255)
                                                                                      gATCGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGttt          >  1:2472774/1‑100 (MQ=255)
                                                                                         cGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTggg       <  1:1063353/100‑1 (MQ=255)
                                                                                         cGTTGGTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTggg       <  1:4171606/100‑1 (MQ=255)
                                                                                             ggTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAAtt   <  2:1919662/100‑1 (MQ=255)
                                                                                             ggTGTCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAAt    >  1:2678602/1‑99 (MQ=255)
                                                                                              gtgtCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAAttt  >  1:2117847/1‑100 (MQ=255)
                                                                                                 |                                                                                                
AGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTT  >  NC_000913/2966485‑2966678

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: