Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,967,771 C→A 6.9% E222* (GAG→TAG)  ptsP ← phosphoenolpyruvate‑protein phosphotransferase PtsP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,967,7710CA 6.9% 212.3 / 3.6 64E222* (GAG→TAG) ptsPphosphoenolpyruvate‑protein phosphotransferase PtsP
Reads supporting (aligned to +/- strand):  ref base C (32/27);  new base A (2/3);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGC  >  NC_000913/2967623‑2967918
                                                                                                                                                    |                                                                                                                                                    
gCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCg                                                                                                                                                     <  1:486290/150‑1 (MQ=255)
 cccACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGt                                                                                                                                                    <  2:1345236/150‑1 (MQ=255)
 cccACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGt                                                                                                                                                    >  1:436402/1‑150 (MQ=255)
 cccACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGt                                                                                                                                                    <  1:728130/150‑1 (MQ=255)
         gCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGct                                                                                                                                           >  2:299647/1‑150 (MQ=255)
           cACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGcttc                                                                                                                                         >  1:723381/1‑150 (MQ=255)
           cACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGcttc                                                                                                                                         <  2:723381/150‑1 (MQ=255)
           cACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGC‑AAACGAGTTTCGCCGctac                                                                                                                                         <  1:1162508/150‑3 (MQ=255)
               gAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTcc                                                                                                                                     >  1:1315498/1‑150 (MQ=255)
               gAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGC‑AAACGAGTTTCGCCGctacagca                                                                                                                                     <  1:718385/150‑7 (MQ=255)
                 gCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAg                                                                                                                                   >  2:172429/1‑150 (MQ=255)
                  ccTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAgc                                                                                                                                  <  2:473764/150‑1 (MQ=255)
                    tttATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAgcgc                                                                                                                                >  2:1344472/1‑150 (MQ=255)
                      tATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGccc                                                                                                                              <  2:1337002/150‑1 (MQ=255)
                        tCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGTAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCg                                                                                                                            >  1:368897/1‑150 (MQ=255)
                         cAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGC‑AAACGAGTTTCGCCGct                                                                                                                                           >  2:291070/1‑134 (MQ=255)
                         cAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGC‑AAACGAGTTTCGCCGct                                                                                                                                           <  1:291070/134‑1 (MQ=255)
                         cAACCTCCGCAAACAATTCGCGACGCAGACTGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCgg                                                                                                                           >  2:530775/1‑150 (MQ=255)
                            ccTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCa                                                                                                                        <  1:1125901/150‑1 (MQ=255)
                                ggCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCg                                                                                                                    >  2:74642/1‑150 (MQ=255)
                                ggCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCg                                                                                                                    >  1:544479/1‑150 (MQ=255)
                                   aaaCAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTc                                                                                                                 <  2:39591/150‑1 (MQ=255)
                                    aaCAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTcg                                                                                                                >  2:1326868/1‑150 (MQ=255)
                                      cAATTCGCGACGCAGCCGGGTGTCCTAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGACGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTACAGCGCCCCGGTCAGTCGTTcgcg                                                                                                              >  2:670629/1‑150 (MQ=255)
                                      cAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTcgcg                                                                                                              >  2:379697/1‑150 (MQ=255)
                                        aTTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGtt                                                                                                            >  1:637429/1‑150 (MQ=255)
                                           cgcgACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCa                                                                                                         >  1:239445/1‑150 (MQ=255)
                                           cgcgACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGAGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCa                                                                                                         >  2:790621/1‑150 (MQ=255)
                                              gACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCagag                                                                                                      <  1:172429/150‑1 (MQ=255)
                                               aCGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGc                                                                                                     <  1:1142044/150‑1 (MQ=255)
                                               aCGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGc                                                                                                     <  1:1407221/150‑1 (MQ=255)
                                                             ccGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTg                                                                                       >  1:263940/1‑150 (MQ=255)
                                                              cGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCTGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGa                                                                                      >  2:969312/1‑150 (MQ=255)
                                                                  aGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTAGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGcc                                                                                  >  2:106788/1‑150 (MQ=255)
                                                                    cAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTg                                                                                <  2:436402/150‑1 (MQ=255)
                                                                     aGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGa                                                                               <  1:404173/150‑1 (MQ=255)
                                                                      ggTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGAt                                                                              >  2:884732/1‑150 (MQ=255)
                                                                                      gAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTaa                                                              >  1:937717/1‑150 (MQ=255)
                                                                                        aaaTAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAg                                                            <  1:530775/150‑1 (MQ=255)
                                                                                        aaaTAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAg                                                            <  2:544479/150‑1 (MQ=255)
                                                                                         aaTAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAgg                                                           >  1:408335/1‑150 (MQ=255)
                                                                                          aTAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGt                                                          >  1:708116/1‑150 (MQ=255)
                                                                                              ccgccgTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACg                                                      <  1:299647/150‑1 (MQ=255)
                                                                                               cgccgcTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTAGATGCCTGATACACCTGTTCCATTAAAGGTAACGt                                                     >  1:1147843/1‑150 (MQ=255)
                                                                                                 ccgTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTgg                                                   <  1:1344472/150‑1 (MQ=255)
                                                                                                     ttCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATc                                               <  1:682063/150‑1 (MQ=255)
                                                                                                      tCTTTTTGTGCACCGGCGGCAATGGGTTTGGTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCATCGCCCCGGTCAGTCGTTCGCGTTCCTGAGCCTGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGGATcc                                              <  2:1028389/150‑1 (MQ=255)
                                                                                                        tttttGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAGCGTGGCATCCTg                                            >  1:1145081/1‑150 (MQ=255)
                                                                                                          tttGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGcc                                          <  1:379697/150‑1 (MQ=255)
                                                                                                             gtgCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGc                                       <  2:1315498/150‑1 (MQ=255)
                                                                                                              tgCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGcc                                      >  1:928650/1‑150 (MQ=255)
                                                                                                                  ccGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCg                                  >  1:72769/1‑150 (MQ=255)
                                                                                                                  ccGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCg                                  >  1:1320010/1‑150 (MQ=255)
                                                                                                                           aaGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGcc                         <  1:1326868/150‑1 (MQ=255)
                                                                                                                            aGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCca                        >  1:1182583/1‑150 (MQ=255)
                                                                                                                              cGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCcaca                      <  1:236483/150‑1 (MQ=255)
                                                                                                                                      gTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTc                                   >  1:723738/1‑129 (MQ=255)
                                                                                                                                      gTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTc                                   <  2:723738/129‑1 (MQ=255)
                                                                                                                                      gTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGtgc              <  1:1345965/150‑1 (MQ=255)
                                                                                                                                      gTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGtgc              >  1:410073/1‑150 (MQ=255)
                                                                                                                                         gcggcgAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGtgctgc           >  1:644062/1‑150 (MQ=255)
                                                                                                                                             cgAAACTCTTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGt       <  2:708116/150‑1 (MQ=255)
                                                                                                                                                aaCTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAAt    <  2:239445/150‑1 (MQ=255)
                                                                                                                                                  cTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATgc  >  1:1341986/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTT‑GCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGC  >  NC_000913/2967623‑2967918

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