Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,894,054 G→T A62D (GCC→GAC)  yieG ← putative membrane/transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,894,0540GT99.1% 1504.4 / ‑1.2 424A62D (GCC→GAC) yieGputative membrane/transport protein
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base T (199/221);  total (201/223)
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 438 total aligned reads displayed.

ATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTAAAACCGGCGATC  >  NC_000913/3893964‑3894143
                                                                                          |                                                                                         
aTAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGtcg                                                                                         >  2:2513517/1‑93 (MQ=255)
 tAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGtcgc                                                                                        >  2:4402421/1‑93 (MQ=255)
  aGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCGGCGATCAGACAGGTAGTGACGAAGACGTCGCt                                                                                       <  1:4216189/93‑1 (MQ=255)
   gCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTg                                                                                      <  2:1995921/93‑1 (MQ=255)
   gCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTg                                                                                      <  1:969742/93‑1 (MQ=255)
      ggTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTa                                                                                   >  2:5188815/1‑93 (MQ=255)
        tGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATc                                                                                 >  1:5258731/1‑93 (MQ=255)
        tGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATc                                                                                 >  1:7027469/1‑93 (MQ=255)
         gCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATcc                                                                                <  2:468857/93‑1 (MQ=255)
          ccAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCa                                                                               <  1:3213046/93‑1 (MQ=255)
            aGTGCAACTGGCAGGTGGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATg                                                                             <  2:3616786/93‑1 (MQ=255)
             gTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGc                                                                            >  2:3453612/1‑93 (MQ=255)
             gTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGc                                                                            >  2:5117022/1‑93 (MQ=255)
              tGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCATACAGGTAGTGACGAAGACGTCGCTGGTATCCATGcc                                                                           >  2:5814327/1‑93 (MQ=255)
              tGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGcc                                                                           <  1:83409/93‑1 (MQ=255)
               gCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCa                                                                          >  2:899756/1‑93 (MQ=255)
               gCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCa                                                                          <  2:7157891/93‑1 (MQ=255)
                 aaCTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGc                                                                        >  1:5548319/1‑93 (MQ=255)
                   cTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCaa                                                                      <  1:437647/93‑1 (MQ=255)
                   cTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCaa                                                                      >  1:7149948/1‑93 (MQ=255)
                     ggCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACg                                                                    >  1:786704/1‑93 (MQ=255)
                     gacAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGAAACg                                                                    <  2:4095081/91‑1 (MQ=255)
                       cAGGTTGGCAAGCAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGcc                                                                  <  1:6002931/93‑1 (MQ=255)
                       cAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGcc                                                                  <  2:3168311/93‑1 (MQ=255)
                        aGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCa                                                                 >  2:6120661/1‑93 (MQ=255)
                           ttggCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGa                                                              >  2:5255634/4‑93 (MQ=255)
                           ttggCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGa                                                              >  2:5380048/4‑93 (MQ=255)
                            tggCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGaa                                                             >  2:815423/3‑93 (MQ=255)
                             ggCAAACAGTCCCATCATAATACTGCCGTATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGCATCCATGCCAGCAACGCCAAGAAt                                                            >  2:1206612/2‑93 (MQ=255)
                             ggCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAAt                                                            >  1:1459764/2‑93 (MQ=255)
                             ggCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAAt                                                            <  1:3619691/92‑1 (MQ=255)
                                aaaCAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTg                                                         <  1:1712643/93‑1 (MQ=255)
                                aaaCAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTg                                                         >  2:1706272/1‑93 (MQ=255)
                                  aCAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCg                                                       >  1:7710122/1‑93 (MQ=255)
                                  aCAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCg                                                       >  1:3783493/1‑93 (MQ=255)
                                    aGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCggg                                                     <  1:661933/93‑1 (MQ=255)
                                    aGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCggg                                                     >  1:152555/1‑93 (MQ=255)
                                    aGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCggg                                                     <  2:7620694/93‑1 (MQ=255)
                                      tCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGtt                                                   >  1:7405724/1‑93 (MQ=255)
                                       cccATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTa                                                  <  1:2241535/93‑1 (MQ=255)
                                       cccATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTa                                                  >  1:6264531/1‑93 (MQ=255)
                                        ccATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTaa                                                 >  2:7883521/1‑93 (MQ=255)
                                           tcatAATACTGCCGAATTCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACaa                                              <  1:3399756/93‑1 (MQ=255)
                                           tcatAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACaa                                              >  2:4747922/1‑93 (MQ=255)
                                           tcatAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACaa                                              <  2:2342624/93‑1 (MQ=255)
                                           tcatAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACaa                                              <  2:6766593/93‑1 (MQ=255)
                                             ataataCTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACaaaa                                            >  1:5058758/1‑93 (MQ=255)
                                              taataCTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGTCTGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGCGTTAACaaaaa                                           >  1:2093662/1‑93 (MQ=255)
                                              taataCTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACaaaaa                                           <  1:1809096/93‑1 (MQ=255)
                                              taataCTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACaaaaa                                           >  2:6962641/1‑93 (MQ=255)
                                               aataCTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAAc                                          <  1:5626873/93‑1 (MQ=255)
                                                  aCTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGAt                                       <  1:6951437/93‑1 (MQ=255)
                                                   cTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATg                                      >  2:3319430/1‑93 (MQ=255)
                                                   cTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATg                                      >  1:3020159/1‑93 (MQ=255)
                                                    tGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATga                                     >  1:2675187/1‑92 (MQ=255)
                                                    tGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGt                                     >  2:2797087/1‑93 (MQ=255)
                                                    tGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGt                                     <  1:105899/93‑1 (MQ=255)
                                                       cGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTaaa                                  <  2:206229/93‑1 (MQ=255)
                                                       cGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTaaa                                  <  1:758437/93‑1 (MQ=255)
                                                        gAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAAc                                 <  2:1483051/93‑1 (MQ=255)
                                                        gAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAAc                                 <  2:6059668/93‑1 (MQ=255)
                                                        gAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAAc                                 >  1:13585/1‑93 (MQ=255)
                                                        gAACGCAGCGATCAGACAGGGAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAAc                                 <  1:4704951/93‑1 (MQ=255)
                                                         aaTGCAGCGGTCAGACATGTAGTGACGAAGACCGCGCTGGTATCCATGCCAGCAACGCCAAGGCTTTTCGGGTTAACCAACACGATGTAAAcc                                >  2:7401878/1‑93 (MQ=255)
                                                         aaTGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAAcc                                <  1:4912117/93‑1 (MQ=255)
                                                         aaTGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAAcc                                >  2:5621946/1‑93 (MQ=255)
                                                          aTGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCa                               >  2:6262892/1‑93 (MQ=255)
                                                           tGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCAt                              >  2:4861217/1‑93 (MQ=255)
                                                           tGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCAt                              >  1:6428123/1‑93 (MQ=255)
                                                           tGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCAt                              <  2:3778745/93‑1 (MQ=255)
                                                           tGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCAt                              <  1:594158/93‑1 (MQ=255)
                                                           tGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCAt                              >  2:7675532/1‑93 (MQ=255)
                                                             cAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCg                            <  2:7522158/93‑1 (MQ=255)
                                                               gCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTc                          >  2:1860768/1‑93 (MQ=255)
                                                                cGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCa                         >  2:1960075/1‑93 (MQ=255)
                                                                 gATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAg                        <  1:6880282/93‑1 (MQ=255)
                                                                 gATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAg                        >  1:2467642/1‑93 (MQ=255)
                                                                 gATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAg                        >  1:5764158/1‑93 (MQ=255)
                                                                  aTCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAgg                       >  2:1319360/1‑93 (MQ=255)
                                                                   tCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGa                      >  2:6441809/1‑93 (MQ=255)
                                                                   tCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGa                      >  2:3059026/1‑93 (MQ=255)
                                                                    cAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGaa                     >  2:297377/1‑93 (MQ=255)
                                                                    cAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGaa                     <  1:5474191/93‑1 (MQ=255)
                                                                     aGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAg                    <  2:6364675/93‑1 (MQ=255)
                                                                       aCAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAggt                  <  1:4068039/93‑1 (MQ=255)
                                                                       aCAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGACGTAAACCATCGTCAGGAAggg                  >  2:1403256/1‑92 (MQ=255)
                                                                         aGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAggtgg                >  1:7806200/1‑93 (MQ=255)
                                                                         aGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAggtgg                >  1:6078086/1‑93 (MQ=255)
                                                                         aGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAggtgg                <  1:1217799/93‑1 (MQ=255)
                                                                         aGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAggtgg                <  1:6742784/93‑1 (MQ=255)
                                                                         aGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAggtgg                >  1:7238589/1‑93 (MQ=255)
                                                                          ggTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAggtggt               >  2:2272802/1‑93 (MQ=255)
                                                                              gTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTaaaa           <  2:6584061/93‑1 (MQ=255)
                                                                              gTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTaaaa           >  2:7323094/1‑93 (MQ=255)
                                                                                gACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTAAAAcc         <  1:4814285/93‑1 (MQ=255)
                                                                                 aCGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTAAAACCg        >  1:7685304/1‑93 (MQ=255)
                                                                                    aaGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTAAAACCGGCg     >  1:3470571/1‑93 (MQ=255)
                                                                                     aGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTACAACCGGCGa    >  1:2354942/1‑93 (MQ=255)
                                                                                     aGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTAAAACCGGCGa    <  1:7485966/93‑1 (MQ=255)
                                                                                       aCGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTAAAACCGGCGATc  >  1:2850254/1‑93 (MQ=255)
                                                                                          |                                                                                         
ATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTAAAACCGGCGATC  >  NC_000913/3893964‑3894143

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