Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,058,585 C→T 30.4% intergenic (‑175/+81) serA ← / ← rpiA D‑3‑phosphoglycerate dehydrogenase/ribose 5‑phosphate isomerase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,058,5850CT30.4% 163.6 / 88.5 122intergenic (‑175/+81)serA/rpiAD‑3‑phosphoglycerate dehydrogenase/ribose 5‑phosphate isomerase, constitutive
Reads supporting (aligned to +/- strand):  ref base C (46/37);  new base T (24/13);  total (72/50)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.56e-01

Only 100 of 125 total aligned reads displayed.

GGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTT  >  NC_000913/3058486‑3058683
                                                                                                   |                                                                                                  
ggCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGAc                                                                                                    <  2:1285771/100‑1 (MQ=255)
 gCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACa                                                                                                   <  2:1913296/100‑1 (MQ=255)
     tatGAACGTTTGCGTCGCGATGGTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTg                                                                                               <  2:1716364/100‑1 (MQ=255)
      atGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGATAATTGc                                                                                              >  2:891265/1‑100 (MQ=255)
       tGAACGTTTGCGTCGCGATGTTGAGGAAATATAACGCAAAGGTAATTCTTGGATGAATATACGGGTTTAGAGGATTTATTTAGAGCAATCGACAAttcc                                                                                              >  2:1351393/1‑97 (MQ=255)
        gAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCt                                                                                            <  1:140644/100‑1 (MQ=255)
         aaCGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGATAATTGCCTg                                                                                           <  2:2423249/100‑1 (MQ=255)
           cGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGt                                                                                         >  1:1597916/1‑100 (MQ=255)
            gTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTa                                                                                        <  2:216613/100‑1 (MQ=255)
               tGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGATAATTGCCTGGTaaaa                                                                                     >  2:1029477/1‑100 (MQ=255)
               tGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTaaaa                                                                                     >  2:2363532/1‑100 (MQ=255)
                gCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAg                                                                                    <  1:1100447/100‑1 (MQ=255)
                    cgcgATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGATAATTGCCTGGTAAAAGCGTg                                                                                <  2:1921725/100‑1 (MQ=255)
                    cgcgATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTg                                                                                <  2:104254/100‑1 (MQ=255)
                     gcgATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGa                                                                               >  1:2187203/1‑100 (MQ=255)
                      cgATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGac                                                                              <  2:2141125/100‑1 (MQ=255)
                      cgATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGac                                                                              <  1:1208355/100‑1 (MQ=255)
                       gATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGaca                                                                             <  2:1558422/100‑1 (MQ=255)
                       gATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGaca                                                                             >  2:2955543/1‑100 (MQ=255)
                         tGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGAcaca                                                                           >  1:2162710/1‑100 (MQ=255)
                          gTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATGTGCTGGATTAGAGGATTTATTTAGAGCAATCGAAAATTGCCTGGTAAAAGCGTGAcaca                                                                           >  2:2555474/1‑99 (MQ=255)
                           ttGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGATAATTGCCTGGTAAAAGCGTGACACATg                                                                         <  1:1414294/100‑1 (MQ=255)
                            tGAGGAAATATCACGCAACGGTAATTATTGGATGAATATGCTGGTTTATAGGATATATTTAGAGGAATCGACAATTGCGTGGTAAAAGCGTGAAACATg                                                                         >  2:700293/1‑99 (MQ=255)
                            tGAGGAAATATCACGCAAAGGTAATTCTTGGTTGACTCTGTTGTTATAAAGGAAAAATTTAGAAAAATCGACAATTGCCTTGTAAAAGCGTGACACATGt                                                                        >  1:62501/1‑100 (MQ=255)
                            tGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGATAATTGCCTGGTAAAAGCGTGACACATGt                                                                        >  2:1928973/1‑100 (MQ=255)
                            tGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGATAATTGCCTGGTAAAAGCGTGACACATGt                                                                        >  1:758536/1‑100 (MQ=255)
                            tGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGt                                                                        >  2:767048/1‑100 (MQ=255)
                             gAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGATAATTGCCTGGTAAAAGCGTGACACATg                                                                         >  1:2847308/1‑98 (MQ=255)
                               ggAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGTAATCGACAATTGCCTGGTAAAAGCGTGACAAATGTCAc                                                                     >  1:982334/1‑100 (MQ=255)
                               ggAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCAc                                                                     >  1:2357368/1‑100 (MQ=255)
                                       cACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCAGGTAAAAGCGTGACACATGTCACCAAATTTa                                                             >  1:454801/1‑100 (MQ=255)
                                          gCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATTGATAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTaa                                                            <  1:2295673/98‑1 (MQ=255)
                                          gCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGATAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATg                                                          <  2:395568/100‑1 (MQ=255)
                                          gCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATg                                                          <  2:1281535/100‑1 (MQ=255)
                                          gCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATg                                                          <  2:1718238/100‑1 (MQ=255)
                                          gCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGAAAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTa                                                             >  2:745049/1‑97 (MQ=255)
                                           cAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGATAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGa                                                         >  2:2289306/1‑100 (MQ=255)
                                               ggTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAagag                                                     >  1:1747701/1‑100 (MQ=255)
                                                   aTTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAAtt                                                 >  2:349928/1‑100 (MQ=255)
                                                     tCTTGGATGAATATGCTGGTTTAGAGGATTTAATTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAAtttt                                               >  1:2891211/1‑100 (MQ=255)
                                                      cTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAAttttt                                              <  1:894335/100‑1 (MQ=255)
                                                      cTTGGATGAATATGCTGGATTAGAGGATTAATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAAtttt                                               >  2:1471498/1‑99 (MQ=255)
                                                       ttGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGATAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAAtttttn                                             <  1:2369443/100‑2 (MQ=255)
                                                       ttGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAAtttttt                                             <  2:2955693/100‑1 (MQ=255)
                                                       ccgatctaaTATGCTGGTTTAGAGGATTTATTTAGAGCAATCGATAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAAtttttt                                             <  2:375718/93‑1 (MQ=255)
                                                         ggATGAATATGCTGGTTTAGAGGAGTTATTTAGAGCAATCGATAATTGCATGGTAAAAGCGTGACACATGTAACCAAATTTAATGAAGAGAATTTTTTTa                                           >  2:1422852/1‑100 (MQ=255)
                                                          gATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTaa                                          >  2:798607/1‑100 (MQ=255)
                                                           aTGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAAc                                         >  1:2220883/1‑100 (MQ=255)
                                                           aTGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATATTTTTAAc                                         >  1:1547999/1‑100 (MQ=255)
                                                            tGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAAc                                         <  1:2521511/99‑1 (MQ=255)
                                                            tGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACg                                        >  2:1642204/1‑100 (MQ=255)
                                                             gAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACgg                                       <  1:2104955/100‑1 (MQ=255)
                                                              aaTATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACggg                                      >  1:1750942/1‑100 (MQ=255)
                                                                tatGCTGGTTTAGAGGATTTATTTAGAGCAATCGATAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACggggg                                    >  1:1425744/1‑100 (MQ=255)
                                                                tatGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACg                                        >  1:515303/1‑96 (MQ=255)
                                                                 atGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGa                                   >  1:1466605/1‑100 (MQ=255)
                                                                  tGCTGGTTTAGAGGATTTATTTAGAGCAATCGATAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAg                                  >  2:1897588/1‑100 (MQ=255)
                                                                   gCTGGTTTAGAGGATTTATTTAGAGCAATCGATAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACTGGGGa                                   >  2:806678/1‑98 (MQ=255)
                                                                   gCTGGTTTAGAGGATTTATTTAGAGCAATCGATAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAg                                  >  1:470247/1‑99 (MQ=255)
                                                                   gCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAgg                                 >  2:1568236/1‑100 (MQ=255)
                                                                    cTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAgg                                 >  1:2782384/1‑99 (MQ=255)
                                                                     tagTTTAGAGGATTTACTTAGAGCAATCGACAATTGCCTCGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAAGt                                <  2:2121137/97‑1 (MQ=255)
                                                                     tGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGaacct                               <  1:2174973/100‑5 (MQ=255)
                                                                         ttAGAGGATTTATTTAGAGCAATCGATAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTcccc                           <  1:479031/100‑1 (MQ=255)
                                                                          tAGAGGATTTATTTAGAGCAATCGATAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTccc                            >  1:2587346/1‑98 (MQ=255)
                                                                          tAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTccccc                          <  2:2043824/100‑1 (MQ=255)
                                                                           agagGATTTATTTAGAGCAATCGATAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCg                         <  1:1159208/100‑1 (MQ=255)
                                                                           agagGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCg                         <  1:1133242/100‑1 (MQ=255)
                                                                           agagGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAAAAATTTTTTAACGGGGGAGGTTCCCCCg                         >  1:1204902/1‑100 (MQ=255)
                                                                            gagGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGt                        >  1:2601733/1‑100 (MQ=255)
                                                                             agGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTc                       <  2:674030/100‑1 (MQ=255)
                                                                              ggATTTATTTAGAGCAATCGATAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAACCTCCCCCGtta                      <  2:1273452/100‑3 (MQ=255)
                                                                              ggATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCa                      >  1:2971734/1‑100 (MQ=255)
                                                                              ggATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCa                      >  2:1468390/1‑100 (MQ=255)
                                                                               gATTTATTTAGAGCAATCGATAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAg                     <  2:1637418/100‑1 (MQ=255)
                                                                               gATTTATTTAGAGCAATCGATAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAg                     <  2:2061921/100‑1 (MQ=255)
                                                                               gATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAg                     >  1:104212/1‑100 (MQ=255)
                                                                               gATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAg                     >  1:1555354/1‑100 (MQ=255)
                                                                                 tttatttaGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGAt                   >  2:406975/1‑100 (MQ=255)
                                                                                    atttaGAGCAATCGATAATTGCCTGGTAAAAGAGTGACACATGTCACCAAATTTAATGAAGAAAATTTTTTTAACGGGGGAGGTTCACCAGTCAGATCAt                >  1:2889262/1‑100 (MQ=255)
                                                                                    atttaGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCAt                >  2:287397/1‑100 (MQ=255)
                                                                                     tttaGAGCAATCGACAATTGCCTGGTAAAATCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTGAACGGGGGAACCTCCCCCGttaaaaaaat               <  1:1255467/100‑10 (MQ=37)
                                                                                     tttaGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCAtn               <  2:1819685/100‑2 (MQ=255)
                                                                                       taGAGCAATCGATAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTc             >  2:1067984/1‑100 (MQ=255)
                                                                                       taGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTc             <  1:2759667/100‑1 (MQ=255)
                                                                                       taGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTc             <  2:2891179/100‑1 (MQ=255)
                                                                                        agagCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTca            <  1:1592796/100‑1 (MQ=255)
                                                                                        agagCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTca            >  1:2567732/1‑100 (MQ=255)
                                                                                        agagCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTca            <  2:1419983/100‑1 (MQ=255)
                                                                                         gagCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTcac           <  2:1709232/100‑1 (MQ=255)
                                                                                           gCAATCGATAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACaa         >  2:2156279/1‑100 (MQ=255)
                                                                                           gCAATCGATAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCACCGTCAGATCATTTCACaa         >  1:1240492/1‑100 (MQ=255)
                                                                                           gCAATCGATAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGAGGGAGGTTCCCCCGTCAGATCATTTCACaa         >  2:2156331/1‑100 (MQ=255)
                                                                                           gCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACaa         >  1:2450334/1‑100 (MQ=255)
                                                                                            cAATCGATAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGAGGAGGTTCCCCAGTCAGATCATTTCACAAt        >  1:2400833/1‑100 (MQ=255)
                                                                                            cAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAAt        <  2:104211/100‑1 (MQ=255)
                                                                                                 gACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGttt   <  2:2357341/100‑1 (MQ=255)
                                                                                                 gACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGttt   >  1:2468624/1‑100 (MQ=255)
                                                                                                   cAATTGCCTGGTAATAGCGTGACACATGTTACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGtttt  <  1:1862234/99‑1 (MQ=255)
                                                                                                   cAATTGCCTGGTAAAAGCGTGAAACATATCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGtttt  >  2:710158/1‑99 (MQ=255)
                                                                                                   |                                                                                                  
GGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTT  >  NC_000913/3058486‑3058683

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