Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,966,413 A→T 78.9% 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,4130AT78.9% 207.7 / 48.8 99D674E (GAT→GAAptsPPEP‑protein phosphotransferase enzyme I; GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base A (9/12);  new base T (46/32);  total (55/44)
Fisher's exact test for biased strand distribution p-value = 2.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.18e-01

CAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGAAGACT  >  NC_000913/2966316‑2966511
                                                                                                 |                                                                                                  
cAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTc                                                                                                  >  1:643804/1‑100 (MQ=255)
 aaTGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTcg                                                                                                 >  2:903570/1‑100 (MQ=255)
 aaTGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTcg                                                                                                 <  1:955017/100‑1 (MQ=255)
    gcgcCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTcgccc                                                                                              >  2:1547874/1‑100 (MQ=255)
     cgccgcAGCAGGTATTTTGCCCGCGCTACTGAACGTCCGTTCATAGACAGATGGCGATTCCCAAGCCCAATGAGGATTGCCACGCACATGGGTTcgcccg                                                                                             <  2:1206665/100‑1 (MQ=255)
     cgccgcAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTcgcccg                                                                                             >  2:1244941/1‑100 (MQ=255)
     cgccgcAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTcgcccg                                                                                             >  2:1484417/1‑100 (MQ=255)
        cgcAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCa                                                                                          >  1:1468944/1‑100 (MQ=255)
           agcagGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATAGGTTCGCCCGCCAt                                                                                         >  1:1574978/1‑98 (MQ=255)
            gcagGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCAtttn                                                                                      <  1:1460710/100‑2 (MQ=255)
            gcagGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTc                                                                                      >  2:61026/1‑100 (MQ=255)
            gcagGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTc                                                                                      <  2:2308067/100‑1 (MQ=255)
              agGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCAc                                                                                    >  2:1930726/1‑100 (MQ=255)
               ggTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCATTTCAcc                                                                                   <  2:1310677/100‑1 (MQ=255)
               ggTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCATTTCAcc                                                                                   <  1:840462/100‑1 (MQ=255)
                gTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCg                                                                                  <  2:1915166/100‑1 (MQ=255)
                gTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCg                                                                                  >  2:2235099/1‑100 (MQ=255)
                 tATTTTGCCCGCGCTACAGAACGTCCGTTCATATACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGAGTTCGCCCGACATTTCACCGc                                                                                 >  1:852039/1‑100 (MQ=255)
                 tATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCATTTCACCGc                                                                                 >  2:720909/1‑100 (MQ=255)
                   ttttGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGcac                                                                               <  2:798818/100‑1 (MQ=255)
                    tttGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGcaca                                                                              >  1:1839713/1‑100 (MQ=255)
                      tGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACaaa                                                                            >  1:1800304/1‑100 (MQ=255)
                         ccGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACgg                                                                         >  1:1112798/1‑100 (MQ=255)
                          cgcgCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGga                                                                        <  2:1914892/100‑1 (MQ=255)
                              cTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATc                                                                    <  2:1606498/100‑1 (MQ=255)
                              cTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATc                                                                    >  2:2293899/1‑100 (MQ=255)
                              cTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCATTTCACCGCACAAACGGAGATc                                                                    <  2:294762/100‑1 (MQ=255)
                              cTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCATTTCACCGCACAAACGGAGATc                                                                    >  1:630454/1‑100 (MQ=255)
                               tACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCg                                                                   >  2:1922640/1‑100 (MQ=255)
                                aCAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACACACGCAGATCg                                                                   >  1:92128/1‑99 (MQ=255)
                                aCAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGa                                                                  >  1:1345634/1‑100 (MQ=255)
                                aCAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGa                                                                  >  1:191280/1‑100 (MQ=255)
                                    aaCGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACTGCACAAACGGAGATCGATTcc                                                              >  2:1456650/1‑100 (MQ=255)
                                     aCGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCa                                                             <  2:1261008/100‑1 (MQ=255)
                                      cGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCAt                                                            <  2:1864539/100‑1 (MQ=255)
                                        tCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGt                                                          <  1:2561279/100‑1 (MQ=255)
                                          cGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTAt                                                        >  2:1672971/1‑100 (MQ=255)
                                              cATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTcc                                                    <  2:2227317/100‑1 (MQ=255)
                                               aTAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCg                                                   <  2:2260076/100‑1 (MQ=255)
                                               aTAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCg                                                   >  2:680412/1‑100 (MQ=255)
                                                tAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGc                                                  >  1:1288718/1‑100 (MQ=255)
                                                tAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGc                                                  >  2:2408936/1‑100 (MQ=255)
                                                tAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGc                                                  <  2:519082/100‑1 (MQ=255)
                                                tAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGc                                                  >  2:955271/1‑100 (MQ=255)
                                                tAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGc                                                  <  2:1287218/100‑1 (MQ=255)
                                                tAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCATTTCACCGCACAAACGGAGACCGATTCCATGCATTTCCg                                                   >  2:2161575/1‑99 (MQ=255)
                                                tAGACAGATGGCGATACCCAAGCCCAATAAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGc                                                  >  1:2475229/1‑100 (MQ=255)
                                                 aGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCt                                                 <  2:920295/100‑1 (MQ=255)
                                                  gACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCtt                                                >  2:2123986/1‑100 (MQ=255)
                                                    cAGATGACGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCATTTCACCGCACAAACGGAGATCGATACCATGTATTTCCGCTTcc                                              <  1:1471763/100‑1 (MQ=255)
                                                       aTGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCgg                                           >  2:2158979/1‑100 (MQ=255)
                                                       aTGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCgg                                           <  1:1595947/100‑1 (MQ=255)
                                                        tGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACACACGGAGATCGATTCCATGTATTTCCGCTTCCCggg                                          <  1:867196/100‑1 (MQ=255)
                                                         ggCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGc                                         <  1:747656/100‑1 (MQ=255)
                                                         ggCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGc                                         <  1:2303208/100‑1 (MQ=255)
                                                           cGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGa                                       >  2:463234/1‑100 (MQ=255)
                                                            gATACCCAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGat                                      <  2:1112802/100‑1 (MQ=255)
                                                              tACCCCAGCCCAATGAGGATTGCAACGCACTTGGGTTCGCCCGCCATTTCACCGCACAAACGTAGATCGATTCCATGTATTTCCGCTTCCCGGGCGa                                       <  1:847028/97‑1 (MQ=255)
                                                              tACCCAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGatca                                    >  2:1550628/1‑100 (MQ=255)
                                                                cccAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGatcatc                                  >  2:1555593/1‑100 (MQ=255)
                                                                cccAAGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGatcatc                                  >  2:2287747/1‑100 (MQ=255)
                                                                cccAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGatcatc                                  <  1:2139471/100‑1 (MQ=255)
                                                                cccAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGatcagc                                  >  1:1568149/1‑98 (MQ=255)
                                                                   aaGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGcc                               <  2:125603/100‑1 (MQ=255)
                                                                    aGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCa                              <  2:603674/100‑1 (MQ=255)
                                                                    aGCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCa                              <  1:2224656/100‑1 (MQ=255)
                                                                     gCCCAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCag                             >  2:1993795/1‑100 (MQ=255)
                                                                        cAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGc                          >  1:1052525/1‑100 (MQ=255)
                                                                        cAATGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGAAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGc                          >  1:2282401/1‑100 (MQ=255)
                                                                        cAATGAGGATTGCCACGCACATGGGATCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGc                          <  1:1263339/100‑1 (MQ=255)
                                                                         aaTGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCt                         >  2:1558738/1‑100 (MQ=255)
                                                                           tGAGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCg                       <  1:1930716/100‑1 (MQ=255)
                                                                             aGGATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTa                     >  2:74466/1‑100 (MQ=255)
                                                                              ggATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTaa                    <  1:2149103/100‑1 (MQ=255)
                                                                               gATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTatc                   <  2:1746348/100‑3 (MQ=255)
                                                                               gATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAAc                   <  2:41784/100‑1 (MQ=255)
                                                                               gATTGCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAAc                   <  2:2173685/100‑1 (MQ=255)
                                                                               gATTGCCACGCACATGGGTTCCGCCGCCATTGCACCGTACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTaan                   <  1:1748713/100‑2 (MQ=255)
                                                                                   gCCACGCACATGGGTTCGCCCGCCATTTCACCGCACACACGGAGATCGATTCCATGTGTTTCCGCTTCCCGGGGGCTCACCGCAAGAGCCCGGAACa                  >  1:540308/1‑97 (MQ=255)
                                                                                   gCCACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTg               >  2:2513477/1‑100 (MQ=255)
                                                                                   gCCACGCACATGGGATCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTg               >  2:1357374/1‑100 (MQ=255)
                                                                                    ccACGCACCGGGGTTCGCCCGCAATTTCTCCGCACCCACGGAGATCGATTCCATGTATTTCCGATTCCCGGGCGATCATCGCCAGAGCGCGTAACAtt                <  1:824733/98‑1 (MQ=255)
                                                                                     cACGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTACCATTGc              >  1:1605904/1‑99 (MQ=255)
                                                                                     cACGCACATGGGATCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCa             <  2:2421337/100‑1 (MQ=255)
                                                                                     cACGCACATGGGATCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGATTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCa             >  1:2078140/1‑100 (MQ=255)
                                                                                      aCGCACATGGGTTCGCCCGCCATTTCCCCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAg            <  1:369501/100‑1 (MQ=255)
                                                                                      aCGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAg            <  2:1812922/100‑1 (MQ=255)
                                                                                      aCGCACATGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAg            <  2:1459873/100‑1 (MQ=255)
                                                                                       cGCACATGGGATCGCCCGCCATTTCACCGCGCACACGGAGATCGATTCCAGGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATGGCAgg           <  1:2113529/100‑1 (MQ=255)
                                                                                         cacaTGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGAt         >  2:438030/1‑100 (MQ=255)
                                                                                         cacaTGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGAt         >  1:289330/1‑100 (MQ=255)
                                                                                           caTGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAg            >  1:597849/1‑95 (MQ=255)
                                                                                           caTGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGACCATCGCCAGAGCTCGTAACATTGCAGGATGa       >  2:1839195/1‑100 (MQ=255)
                                                                                           caTGGGATCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGa       <  1:1749822/100‑1 (MQ=255)
                                                                                             tGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGaa      >  2:1772570/1‑99 (MQ=255)
                                                                                             tGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGAAg     >  1:1475120/1‑100 (MQ=255)
                                                                                             tGGGTTCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCCCGTAACATTGCCGGATGAAg     >  2:1772501/1‑100 (MQ=255)
                                                                                              gggATCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGAAGa    >  2:104895/1‑100 (MQ=255)
                                                                                                gATCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGAAGACt  >  2:2505533/1‑100 (MQ=255)
                                                                                                 |                                                                                                  
CAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCATTTCACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGAAGACT  >  NC_000913/2966316‑2966511

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