Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 968,384 T→C S266P (TCT→CCT)  ycaH → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913968,3840TC99.2% 406.5 / NA 128S266P (TCT→CCT) ycaHhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (67/60);  minor base G (1/0);  total (68/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 130 total aligned reads displayed.

GCGGGGATTGGGCATCCGCCGCGCTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGTCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGAAATGCCGGGCCTTT  >  NC_000913/968294‑968473
                                                                                          |                                                                                         
gCGGGGATTGGGCATCCGCCGCGCTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCt                                                                                         <  2:1810630/93‑1 (MQ=255)
   gggATTGGGCATCCGCCGCGCTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTg                                                                                      >  2:230125/1‑93 (MQ=255)
   gggATTGGGCATCCGCCGCGCTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTg                                                                                      >  2:3391457/1‑93 (MQ=255)
        tGGGCATCCGCCGCGCTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCa                                                                                 >  2:101033/1‑93 (MQ=255)
         gggCATCCGCCGCTCTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCAt                                                                                <  1:34800/93‑1 (MQ=255)
         gggCATCCGCCGCGCTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCAt                                                                                <  1:3707373/93‑1 (MQ=255)
           gCATCCGCCGCGCTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGc                                                                              <  2:3582436/93‑1 (MQ=255)
            cATCCGCCGCGCTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCg                                                                             <  1:2144926/93‑1 (MQ=255)
              tCCGCCGCGCTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGa                                                                           <  2:4580941/93‑1 (MQ=255)
              tCCGCCGCGCTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGa                                                                           <  2:3023789/93‑1 (MQ=255)
              tCCGCCGCGCTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGa                                                                           <  1:2773131/93‑1 (MQ=255)
                cgccgcGCTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCAGATCATCAGCCTTCGAACCATGCGGATg                                                                         >  1:521683/1‑93 (MQ=255)
                  ccgcTCTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTc                                                                       <  2:2060544/93‑1 (MQ=255)
                   cgcgCTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCa                                                                      <  1:2953919/93‑1 (MQ=255)
                   cgcgCTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCa                                                                      >  1:2836184/1‑93 (MQ=255)
                   cgcgCTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCa                                                                      <  1:139693/93‑1 (MQ=255)
                    gcgcTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAg                                                                     <  2:1620595/93‑1 (MQ=255)
                    gcgcTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAg                                                                     >  1:3498482/1‑93 (MQ=255)
                     cgcTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGt                                                                    >  2:2323956/1‑93 (MQ=255)
                      gcTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTCCCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTg                                                                   <  1:4026419/93‑1 (MQ=255)
                      gcTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTg                                                                   >  2:4374881/1‑93 (MQ=255)
                      gcTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTg                                                                   >  2:1187583/1‑93 (MQ=255)
                       cTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGc                                                                  <  2:2089138/93‑1 (MQ=255)
                       cTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGc                                                                  >  2:1371206/1‑93 (MQ=255)
                       cTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGc                                                                  >  1:4452842/1‑93 (MQ=255)
                         tttttGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGt                                                                >  1:515716/1‑93 (MQ=255)
                            ttGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTgg                                                             <  2:1588524/93‑1 (MQ=255)
                             tGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGt                                                            >  1:2580521/1‑93 (MQ=255)
                              gCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTa                                                           <  1:2763599/93‑1 (MQ=255)
                              gCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTa                                                           <  2:4330123/93‑1 (MQ=255)
                              gCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTa                                                           >  2:4429911/1‑93 (MQ=255)
                               ccACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTaa                                                          <  2:2610879/93‑1 (MQ=255)
                                cACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAg                                                         <  2:2074203/93‑1 (MQ=255)
                                 aCGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAgc                                                        <  2:379883/93‑1 (MQ=255)
                                 aCGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAgc                                                        <  1:1099849/93‑1 (MQ=255)
                                      gAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCgg                                                   >  1:1655838/1‑93 (MQ=255)
                                      gAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCgg                                                   >  2:577881/1‑93 (MQ=255)
                                      gAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCgg                                                   >  2:62427/1‑93 (MQ=255)
                                       aaGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCggg                                                  <  2:4015917/93‑1 (MQ=255)
                                       aaGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCggg                                                  >  1:1517992/1‑93 (MQ=255)
                                        aGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGc                                                 >  2:2834350/1‑93 (MQ=255)
                                        aGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGc                                                 >  2:1021275/1‑93 (MQ=255)
                                         gATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCa                                                >  1:296537/1‑93 (MQ=255)
                                         gATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCa                                                >  1:3729844/1‑93 (MQ=255)
                                         gATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCa                                                >  2:1338315/1‑93 (MQ=255)
                                         gATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCa                                                >  1:2391431/1‑93 (MQ=255)
                                           tgtgtgGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCaaa                                              <  1:1986285/93‑1 (MQ=255)
                                           tgtgtgGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCaaa                                              >  1:3970157/1‑93 (MQ=255)
                                            gtgtgGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAAc                                             <  2:2193714/93‑1 (MQ=255)
                                            gtgtgGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAAc                                             <  2:2318324/93‑1 (MQ=255)
                                            gtgtgGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAAc                                             <  1:1986650/93‑1 (MQ=255)
                                              gtgGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGc                                           >  1:1094914/1‑93 (MQ=255)
                                              gtgGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGc                                           <  1:3797934/93‑1 (MQ=255)
                                              gtgGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGc                                           >  2:104622/1‑93 (MQ=255)
                                               tgGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCt                                          <  1:1699254/93‑1 (MQ=255)
                                               tgGCGTACAACCGCAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCt                                          <  2:3927155/93‑1 (MQ=255)
                                                 gCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTgg                                        <  1:2625123/93‑1 (MQ=255)
                                                   gTACAACCGGAAAAATGTGTACCGCTGGCCGATTATCAGGCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAACCGCTGGTa                                      >  1:802246/1‑93 (MQ=255)
                                                   gTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTa                                      >  2:39539/1‑93 (MQ=255)
                                                    tACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTATGCGCCGGGCAAACGCTGGTaa                                     >  2:2640917/1‑93 (MQ=255)
                                                    tACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTaa                                     >  2:4241186/1‑93 (MQ=255)
                                                    tACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTaa                                     >  1:3246171/1‑93 (MQ=255)
                                                     aCAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAAt                                    <  2:1025073/93‑1 (MQ=255)
                                                     aCAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAAt                                    <  1:510297/93‑1 (MQ=255)
                                                       aaCCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGa                                  >  1:4604436/1‑93 (MQ=255)
                                                            gAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGaa                             >  1:3407930/1‑93 (MQ=255)
                                                            gAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGaa                             >  1:1746679/1‑93 (MQ=255)
                                                              aaaaTGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGaaaa                           <  2:2544201/93‑1 (MQ=255)
                                                               aaaTGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGaaaaa                          >  2:557984/1‑93 (MQ=255)
                                                               aaaTGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGaaaaa                          <  1:4497324/93‑1 (MQ=255)
                                                                 aTGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGa                        <  1:4553132/93‑1 (MQ=255)
                                                                 aTGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGa                        <  1:4374587/93‑1 (MQ=255)
                                                                  tgtgTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGAt                       >  2:1661122/1‑93 (MQ=255)
                                                                   gtgtACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATg                      <  1:1094042/93‑1 (MQ=255)
                                                                   gtgtACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATg                      <  1:2374643/93‑1 (MQ=255)
                                                                    tgtACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGTTAATGACTGAAAAAGATGc                     <  2:1876236/93‑1 (MQ=255)
                                                                    tgtACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGc                     <  2:4080761/93‑1 (MQ=255)
                                                                      tACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCgg                   <  1:3427947/93‑1 (MQ=255)
                                                                      tACCGCTGGCCGATCAGCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGTATACGCTGGTAATGACTGAAAAATATGCgg                   <  2:2487391/93‑1 (MQ=255)
                                                                          gCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGaa               <  2:1214991/93‑1 (MQ=255)
                                                                           cTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGaaa              >  1:2943916/1‑93 (MQ=255)
                                                                            tGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGAAAt             >  2:2021880/1‑93 (MQ=255)
                                                                            tGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGAAAt             >  2:4629777/1‑93 (MQ=255)
                                                                              gCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGAAATGc           <  2:1395006/93‑1 (MQ=255)
                                                                               ccGATCATCAGCCTTTGAACCATGCGGATGTCAGTTCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGAAATGcc          <  2:2419388/93‑1 (MQ=255)
                                                                               cagATCAGCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGTTAAGCGCTGGGCAAACGCGGGTAATGGCTGAAAAATATTCCGTGAAACGcc          <  2:3901550/91‑1 (MQ=255)
                                                                                cGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGAAATGCCg         >  1:313655/1‑93 (MQ=255)
                                                                                 gATCATCAGCCTTTGAACCATGCGTATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGAAATGCCgg        <  2:3323200/93‑1 (MQ=255)
                                                                                 gATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGAAATGCCgg        >  1:844136/1‑93 (MQ=255)
                                                                                  atcatcAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGAAATGCCggg       >  2:1611045/1‑93 (MQ=255)
                                                                                  atcatcAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGAAATGCCggg       <  1:4177419/93‑1 (MQ=255)
                                                                                   tcatcaGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGAAATGCCGGGc      >  1:2201818/1‑93 (MQ=255)
                                                                                   tcatcaGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGAAATGCCGGGc      >  2:2649199/1‑93 (MQ=255)
                                                                                    catcaGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGAAATGCCGGGcc     <  2:4627667/93‑1 (MQ=255)
                                                                                     atcaGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGAAATGCCGGGCCt    >  1:1587860/1‑93 (MQ=255)
                                                                                      tcaGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGAAATGCCGGGCCtt   >  1:4267575/1‑93 (MQ=255)
                                                                                      tcaGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGAAATGCCGGGCCtt   >  1:1936467/1‑93 (MQ=255)
                                                                                       caGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGAAATGCCGGGCCttt  <  2:1028218/93‑1 (MQ=255)
                                                                                       caGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGAAATGCCGGGCCttt  <  1:1029505/93‑1 (MQ=255)
                                                                                       caGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGAAATGCCGGGCCttt  <  1:2302564/93‑1 (MQ=255)
                                                                                          |                                                                                         
GCGGGGATTGGGCATCCGCCGCGCTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGTCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGAAATGCCGGGCCTTT  >  NC_000913/968294‑968473

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