Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,484,088 (C)6→7 coding (887/1023 nt) yheT → putative hydrolase YheT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,484,0821.C98.7% 951.0 / 2.8 235coding (881/1023 nt)yheTputative hydrolase YheT
Reads supporting (aligned to +/- strand):  ref base . (0/3);  new base C (112/120);  total (112/123)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 267 total aligned reads displayed.

GCCATGCCGATGCTGAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCT‑‑CCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATCCGCAAATGTGGCTGGAGTCA  >  NC_000913/3483982‑3484176
                                                                                                      |                                                                                                
 ccATGCCGATGCTGAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGtct‑‑c                                                                                                <  1:4548112/101‑1 (MQ=255)
  cATGCCGATGCTGAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCT‑‑cc                                                                                              >  2:4548545/1‑101 (MQ=255)
   aTGCCGATGCTGAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCT‑‑ccc                                                                                             >  1:4548546/1‑101 (MQ=255)
      ccGATGCTGAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCT‑‑cccccc                                                                                          >  2:4548550/1‑101 (MQ=255)
      ccGATGCTGAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCT‑‑cccccc                                                                                          <  2:4548196/101‑1 (MQ=255)
      ccGATGCTGAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCT‑‑cccccc                                                                                          <  2:4548208/101‑1 (MQ=255)
       cGATGCTGAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCT‑‑ccccccc                                                                                         <  1:4547871/101‑2 (MQ=255)
        gATGCTGAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCT‑‑cccccccg                                                                                        >  1:4548554/1‑99 (MQ=255)
        gATGCTGAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCT‑‑cccccccg                                                                                        <  1:4548222/101‑3 (MQ=255)
         aTGCTGAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCT‑‑cccccccgc                                                                                       >  2:4548555/1‑98 (MQ=255)
             tGAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGt                                                                                    >  1:4548559/1‑101 (MQ=255)
             tGAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGt                                                                                    <  2:4547948/101‑1 (MQ=255)
                 ccGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAg                                                                                <  2:4548170/101‑1 (MQ=255)
                 ccGGATCGCCAAACCGACGCTGATTATTCACGCAAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAg                                                                                <  1:4548165/101‑1 (MQ=255)
                     aTCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCT‑‑cccccc                                                                                          >  2:4548562/1‑86 (MQ=255)
                      tCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCa                                                                           <  2:4548307/101‑1 (MQ=255)
                      tCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCa                                                                           >  1:4548563/1‑101 (MQ=255)
                      tCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCa                                                                           >  2:4548565/1‑101 (MQ=255)
                      tCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCa                                                                           <  1:4548311/101‑1 (MQ=255)
                        gCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAac                                                                         >  2:4548568/1‑101 (MQ=255)
                         ccAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAact                                                                        <  1:4548182/101‑1 (MQ=255)
                          cAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCCGCAGGTGGAGTATCAact                                                                        <  1:4548139/101‑1 (MQ=255)
                           aaaCCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAactga                                                                      >  1:4548573/1‑101 (MQ=255)
                           aaaCCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAactga                                                                      >  1:4548571/1‑101 (MQ=255)
                            aaCCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAactgac                                                                     <  1:4548262/101‑1 (MQ=255)
                                gACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGaa                                                                 <  1:4548281/101‑1 (MQ=255)
                                gACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGaa                                                                 >  1:4548577/1‑101 (MQ=255)
                                   gCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACAt                                                              >  2:4548584/1‑101 (MQ=255)
                                   gCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACAt                                                              <  1:4548156/101‑1 (MQ=255)
                                   gCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACAt                                                              >  2:4548582/1‑101 (MQ=255)
                                    cTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATg                                                             <  1:4548447/101‑1 (MQ=255)
                                     tGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAact                                                                        <  2:4548137/89‑1 (MQ=255)
                                     tGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATgg                                                            >  1:4548588/1‑101 (MQ=255)
                                     tGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATgg                                                            >  2:4548587/1‑101 (MQ=255)
                                      gATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGc                                                           <  1:4548147/101‑1 (MQ=255)
                                      gATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGc                                                           <  1:4548357/101‑1 (MQ=255)
                                         tattCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGt                                                        >  1:4548591/1‑101 (MQ=255)
                                         tattCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGt                                                        <  2:4548213/101‑1 (MQ=255)
                                         tattCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGt                                                        <  1:4548075/101‑1 (MQ=255)
                                            tCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCAt                                                     <  1:4548054/101‑1 (MQ=255)
                                              aCGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGt                                                   >  1:4548595/1‑101 (MQ=255)
                                                gCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTg                                                 <  2:4548111/101‑1 (MQ=255)
                                                 ccAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTgg                                                >  2:4548599/1‑101 (MQ=255)
                                                  cAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGc                                               >  1:4548603/1‑101 (MQ=255)
                                                  cAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGc                                               <  1:4548150/101‑1 (MQ=255)
                                                   aaaGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAg                                                                                <  2:4548060/67‑1 (MQ=255)
                                                   aaaGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCt                                              >  1:4548605/1‑101 (MQ=255)
                                                   aaaGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCt                                              >  2:4548606/1‑101 (MQ=255)
                                                    aaGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCtt                                             <  1:4548083/101‑1 (MQ=255)
                                                    aaGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCtt                                             <  1:4548325/101‑1 (MQ=255)
                                                    aaGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCtt                                             >  2:4548609/1‑101 (MQ=255)
                                                       aCGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTAt                                          <  1:4548189/101‑1 (MQ=255)
                                                         gATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTg                                        <  2:4548577/101‑1 (MQ=255)
                                                         gATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTg                                        <  1:4548288/101‑1 (MQ=255)
                                                          aTCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTgg                                       >  1:4548615/1‑101 (MQ=255)
                                                          aTCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTgg                                       <  2:4548131/101‑1 (MQ=255)
                                                           tCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGc                                      >  2:4548616/1‑101 (MQ=255)
                                                            ccGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCaa                                                                          >  2:4548620/1‑64 (MQ=255)
                                                            ccGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCg                                     <  1:4548246/101‑1 (MQ=255)
                                                             cGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGGGGAg                                                                                >  1:4548621/1‑57 (MQ=255)
                                                              gTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGt                                   <  2:4548221/101‑1 (MQ=255)
                                                              gTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGt                                   <  2:4548227/101‑1 (MQ=255)
                                                              gTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGt                                   <  2:4548243/101‑1 (MQ=255)
                                                               tttATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAactgact                                                                    <  1:4548484/67‑1 (MQ=255)
                                                               tttATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTa                                  >  1:4548624/1‑101 (MQ=255)
                                                                 tATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACa                                <  1:4547976/101‑1 (MQ=255)
                                                                 tATGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACa                                >  1:4548628/1‑101 (MQ=255)
                                                                  aTGGATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACAt                               >  2:4548630/1‑101 (MQ=255)
                                                                    ggATCATCAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTAc                                 >  1:4548633/1‑97 (MQ=255)
                                                                      atcatcAGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTAATTGGCGGTACATTACt                           <  2:4548337/101‑1 (MQ=255)
                                                                       tcatcaGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACtt                          <  1:4548251/101‑1 (MQ=255)
                                                                        catcaGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTc                         >  2:4548635/1‑101 (MQ=255)
                                                                        catcaGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTc                         <  1:4548371/101‑1 (MQ=255)
                                                                         atcaGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCa                        >  1:4548637/1‑101 (MQ=255)
                                                                          tcaGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCAt                       >  2:4548639/1‑101 (MQ=255)
                                                                           caGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATc                      <  1:4548301/101‑1 (MQ=255)
                                                                           caGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATc                      <  1:4548403/101‑1 (MQ=255)
                                                                            aGGTGATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATcc                     <  2:4548273/101‑1 (MQ=255)
                                                                             ggTGATCCCGAAACCGGAAAGTCGC‑‑CCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATCCg                    <  1:4548204/101‑1 (MQ=255)
                                                                                gATCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATCCGCaa                 >  1:4548643/1‑101 (MQ=255)
                                                                                 aTCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATCCGCaaa                >  1:4548647/1‑101 (MQ=255)
                                                                                  tCCCGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATCCGCAAAt               >  2:4548648/1‑101 (MQ=255)
                                                                                    ccGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATCCGCAAAtgt             <  1:4548350/101‑1 (MQ=255)
                                                                                    ccGAAACCGGAAAGTCT‑‑CACCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATCCGCAAAtgt             <  1:4548133/101‑1 (MQ=255)
                                                                                     cGAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATCCGCAAAtgtg            >  1:4548651/1‑101 (MQ=255)
                                                                                      gAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATCCGCAAAtg              <  1:4548177/98‑1 (MQ=255)
                                                                                      gAAACCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATCCGCAAATGTgg           >  1:4548655/1‑101 (MQ=255)
                                                                                       aaaCCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATCCGCAAATGTGGc          >  2:4548658/1‑101 (MQ=255)
                                                                                       aaaCCGGAAAGTCTCCCCCCCGCAGGTGGAGCATCAACTGACTGAACATGGCGGTCATGTTg                                                 >  1:4548659/1‑62 (MQ=255)
                                                                                        aaCCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATCCGCAAATGTGGCt         <  2:4548347/101‑1 (MQ=255)
                                                                                         aCCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATCCGCAAATGTGGCTg        <  2:4548446/101‑1 (MQ=255)
                                                                                         aCCGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATCCGCAAATGTGGCTg        <  1:4548201/101‑1 (MQ=255)
                                                                                          ccGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATCCGCAAATGTGGCTgg       >  2:4548662/1‑101 (MQ=255)
                                                                                          ccGGAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATCCGCAAATGTGGCTgg       >  1:4548664/1‑101 (MQ=255)
                                                                                            ggAAAGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGGACATTACTTCAt                       >  2:4548668/1‑83 (MQ=255)
                                                                                              aaaGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATCCGCAAATGTGGCTgg       >  2:4548672/1‑97 (MQ=255)
                                                                                              aaaGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATCCGCAAATGTGGCTGGAGTc   >  1:4548671/1‑101 (MQ=255)
                                                                                              aaaGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATCCGCAAATGTGGCTGGAGTc   <  2:4548367/101‑1 (MQ=255)
                                                                                               aaGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATCCGCAAATGTGGCTGGAGTCa  <  2:4548287/101‑1 (MQ=255)
                                                                                               aaGTCTCCCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATCCGCAAATGTGGCTGGAGTCa  <  2:4548418/101‑1 (MQ=255)
                                                                                                      |                                                                                                
GCCATGCCGATGCTGAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCT‑‑CCCCCCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATCCGCAAATGTGGCTGGAGTCA  >  NC_000913/3483982‑3484176

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