Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,966,413 A→T D674E (GAT→GAA 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,4130AT99.1% 391.7 / ‑4.3 113D674E (GAT→GAAptsPPEP‑protein phosphotransferase enzyme I; GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (58/54);  total (59/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

Only 100 of 113 total aligned reads displayed.

CAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGAAG  >  NC_000913/2966316‑2966508
                                                                                                 |                                                                                               
cAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTc                                                                                               <  2:1865345/100‑1 (MQ=255)
  aTGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTcgc                                                                                             >  2:1505732/1‑100 (MQ=255)
   tGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTcgcc                                                                                            >  1:2240188/1‑100 (MQ=255)
    gtgccgcAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTcgccc                                                                                           <  1:1917173/98‑1 (MQ=255)
     cgccgcAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTcgcccg                                                                                          <  2:2807935/100‑1 (MQ=255)
     cgccgcAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTcgcccg                                                                                          >  1:35957/1‑100 (MQ=255)
      gccgcAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTcgcccgc                                                                                         >  1:2022205/1‑100 (MQ=255)
      gccgcAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTcgcccgc                                                                                         >  1:1838630/1‑100 (MQ=255)
         gcagcaGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCAt                                                                                      <  2:1064054/100‑1 (MQ=255)
          cagcagGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCAtt                                                                                     >  1:1500124/1‑100 (MQ=255)
            gcagGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTc                                                                                   <  2:1315985/100‑1 (MQ=255)
            gcagGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTc                                                                                   <  1:231781/100‑1 (MQ=255)
                gTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCg                                                                               <  2:1600863/100‑1 (MQ=255)
                 tATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGc                                                                              <  2:2750671/100‑1 (MQ=255)
                  aTTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGca                                                                             >  1:2543128/1‑100 (MQ=255)
                  aTTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGca                                                                             >  1:2436244/1‑100 (MQ=255)
                  aTTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGca                                                                             >  1:2523719/1‑100 (MQ=255)
                       gCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAAc                                                                        >  1:2176305/1‑100 (MQ=255)
                        cccGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACg                                                                       <  2:1850975/100‑1 (MQ=255)
                        cccGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACg                                                                       >  2:677997/1‑100 (MQ=255)
                         ccGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACgg                                                                      <  2:1806248/100‑1 (MQ=255)
                          cgcgCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGga                                                                     <  2:1920124/100‑1 (MQ=255)
                          cgcgCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGga                                                                     >  1:2150472/1‑100 (MQ=255)
                            cgcTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGgaga                                                                   >  2:127503/1‑100 (MQ=255)
                             gcTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGAt                                                                  <  2:675902/100‑1 (MQ=255)
                              cTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATc                                                                 <  1:2094773/100‑1 (MQ=255)
                              cTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATc                                                                 <  1:1664416/100‑1 (MQ=255)
                              cTACAGAACGACCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATc                                                                 <  2:1898381/100‑1 (MQ=255)
                                aCAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGa                                                               >  1:2053374/1‑100 (MQ=255)
                                 cAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGAt                                                              <  1:1093893/100‑1 (MQ=255)
                                   gAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTc                                                            >  1:1500783/1‑100 (MQ=255)
                                    aaCGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTcc                                                           <  2:2512414/100‑1 (MQ=255)
                                     aCGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTccn                                                          <  1:1339025/100‑2 (MQ=255)
                                     aCGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCa                                                          <  2:673601/100‑1 (MQ=255)
                                         ccGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTa                                                      >  2:441931/1‑100 (MQ=255)
                                         ccGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTa                                                      <  2:1454307/100‑1 (MQ=255)
                                         ccGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTa                                                      <  1:127503/100‑1 (MQ=255)
                                          ngTTCATAGACAGATGGCGATACCCAAGCCCAATNAGGNTTGCNACNCACANGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTAt                                                     >  1:1325574/2‑100 (MQ=255)
                                           gTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTAtt                                                    <  2:1792342/100‑1 (MQ=255)
                                           gTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTAtt                                                    >  1:792444/1‑100 (MQ=255)
                                             tCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTc                                                  >  1:551424/1‑100 (MQ=255)
                                               aTAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCg                                                <  1:2829681/100‑1 (MQ=255)
                                                tAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGc                                               >  1:2756301/1‑100 (MQ=255)
                                                tAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGc                                               <  2:327049/100‑1 (MQ=255)
                                                tAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGc                                               <  1:257572/100‑1 (MQ=255)
                                                 aGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCt                                              >  2:170109/1‑100 (MQ=255)
                                                 aGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCt                                              >  2:848925/1‑100 (MQ=255)
                                                 aGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCt                                              <  1:244067/100‑1 (MQ=255)
                                                   aCAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTc                                            >  2:1049881/1‑100 (MQ=255)
                                                   aCAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTc                                            >  2:1051487/1‑100 (MQ=255)
                                                    cAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTcc                                           >  2:542635/1‑100 (MQ=255)
                                                    cAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTc                                            >  2:539707/1‑99 (MQ=255)
                                                     aGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGCGAACGATTCCATGTATTTCCGCTTccc                                          >  2:872708/1‑100 (MQ=255)
                                                      gATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCg                                         <  2:2142703/100‑1 (MQ=255)
                                                      gATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCg                                         >  1:911436/1‑100 (MQ=255)
                                                       atgtccaTACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCgg                                        <  2:2469815/94‑1 (MQ=255)
                                                       aTGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGTGTTCGCCGGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCgg                                        >  1:2616398/1‑100 (MQ=255)
                                                       aTGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCgg                                        >  2:318966/1‑100 (MQ=255)
                                                       aTGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCgg                                        <  1:1277614/100‑1 (MQ=255)
                                                         ggCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGc                                      <  1:405717/100‑1 (MQ=255)
                                                         ggCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCACGTATTTCCGCTTCCCGGGc                                      >  2:1516567/1‑100 (MQ=255)
                                                             aTACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGatc                                  <  2:2873983/100‑1 (MQ=255)
                                                               aCCCCAGCCCAATGAGGATTGCCACGCACATGGGTTAGCCAGCCATTTCACCGCACAAACGGAGATCGAATCAATGTATTTCCACTTCCCGGGCGa                                    >  2:2825843/1‑96 (MQ=255)
                                                               aCCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGatcat                                >  2:703167/1‑100 (MQ=255)
                                                                cccAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGatcatc                               >  1:1271267/1‑100 (MQ=255)
                                                                cccAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGatcatc                               >  2:684497/1‑100 (MQ=255)
                                                                cccAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGatcatc                               >  2:2583296/1‑100 (MQ=255)
                                                                   aaGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGcc                            >  2:1129291/1‑100 (MQ=255)
                                                                    aGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCa                           <  1:819932/100‑1 (MQ=255)
                                                                     gCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCag                          <  2:32117/100‑1 (MQ=255)
                                                                        cAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGc                       >  1:1299941/1‑100 (MQ=255)
                                                                        cAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGc                       <  2:2176292/100‑1 (MQ=255)
                                                                          aTGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATGATCGcccgggaag                     <  1:1190811/100‑8 (MQ=255)
                                                                           tGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCg                    <  2:309956/100‑1 (MQ=255)
                                                                              ggATTGCCACGCACATGGGATCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTaa                 >  1:208280/1‑100 (MQ=255)
                                                                                aTTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACa               >  1:2661054/1‑100 (MQ=255)
                                                                                   gCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTg            <  2:1838617/100‑1 (MQ=255)
                                                                                   gCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTg            >  1:215429/1‑100 (MQ=255)
                                                                                   gCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTg            >  1:646668/1‑100 (MQ=255)
                                                                                   gCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTg            >  2:112015/1‑100 (MQ=255)
                                                                                    ccACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGc           <  2:35957/100‑1 (MQ=255)
                                                                                     cACGCACATGGGTTCGCCCGCCATTTCATCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCa          <  1:1505740/100‑1 (MQ=255)
                                                                                     cACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCa          <  1:1817022/100‑1 (MQ=255)
                                                                                     cACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCa          <  1:287980/100‑1 (MQ=255)
                                                                                     cACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCa          >  2:492180/1‑100 (MQ=255)
                                                                                      aCGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAg         <  1:492182/100‑1 (MQ=255)
                                                                                       cGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAgg        >  1:259937/1‑100 (MQ=255)
                                                                                       cGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAgg        <  1:2107597/100‑1 (MQ=255)
                                                                                       cGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAg         <  2:2202910/99‑1 (MQ=255)
                                                                                        gCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGa       <  1:2226724/100‑1 (MQ=255)
                                                                                        gCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGa       <  2:1271262/100‑1 (MQ=255)
                                                                                          acaTGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATg     <  2:291611/100‑1 (MQ=255)
                                                                                           caTGGGTTCGCCCGCCATTTCACCGCACAAACGGCGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCCCAAGAGCCCGTAACATTGCAGGATGa    >  1:2836491/1‑100 (MQ=255)
                                                                                           caTGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGa    >  1:319729/1‑100 (MQ=255)
                                                                                           caTGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGa    <  2:1500116/100‑1 (MQ=255)
                                                                                           caTGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGa    >  1:906845/1‑100 (MQ=255)
                                                                                           caTGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGa    <  2:1431385/100‑1 (MQ=255)
                                                                                             tGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGAAg  <  2:428790/100‑1 (MQ=255)
                                                                                             tGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGAAg  >  2:2416122/1‑100 (MQ=255)
                                                                                             tGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGAAg  >  1:1786893/1‑100 (MQ=255)
                                                                                                 |                                                                                               
CAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGAAG  >  NC_000913/2966316‑2966508

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