Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,866,015 G→A R87C (CGC→TGC)  yidE ← putative transport protein YidE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,866,0150GA99.7% 1370.9 / ‑4.9 321R87C (CGC→TGC) yidEputative transport protein YidE
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (149/171);  total (150/171)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 330 total aligned reads displayed.

CAGTGGAATATCAAACAGTTTATGCAGGATGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCGCAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACAGGATCAGGCCAAATTCCTGAATAACATGCAGCATATCGCTACTT  >  NC_000913/3865917‑3866112
                                                                                                  |                                                                                                 
 aGTGGAATATCAAACAGTTTATGCAGGATGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACaa                                                                                                <  1:6584519/101‑1 (MQ=255)
 aGTGGAATATCAAACAGTTTATGCAGGATGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACaa                                                                                                >  2:6585068/1‑101 (MQ=255)
 aGTGGAATATCAAACAGTTTATGCAGGATGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACaa                                                                                                <  2:6584513/101‑1 (MQ=255)
  gTGGAATATCAAACAGTTTATGCAGGATGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAAt                                                                                               >  1:6585069/1‑101 (MQ=255)
   tGGAATATCAAACAGTTTATGCAGGATGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATg                                                                                              >  2:6585072/1‑101 (MQ=255)
     gAATATCAAACAGTTTATGCAGGATGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAg                                                                                            <  1:6584831/101‑1 (MQ=255)
     gAATATCAAACAGTTTATGCAGGATGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAg                                                                                            <  1:6584691/101‑1 (MQ=255)
      aaTATCAAACAGTTTATGCAGGATGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAgg                                                                                           >  1:6585077/1‑101 (MQ=255)
       atatCAAACAGTTTATGCAGGATGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGc                                                                                          <  2:6584727/101‑1 (MQ=255)
        tatCAAACAGTTTATGCAGGATGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCa                                                                                         >  1:6585080/1‑101 (MQ=255)
         atCAAACAGTTTATGCAGGATGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCaa                                                                                        <  2:6584584/101‑1 (MQ=255)
         atCAAACAGTTTATGCAGGATGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCaa                                                                                        <  1:6584628/101‑1 (MQ=255)
          tCAAACAGTTTATGCAGGATGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCaaa                                                                                       >  2:6585084/1‑101 (MQ=255)
           cAAACAGTTTATGCAGGATGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAaag                                                                                      <  1:6584744/101‑1 (MQ=255)
           cAAACAGTTTATGCAGGATGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAaag                                                                                      >  1:6585085/1‑101 (MQ=255)
           cAAACAGTTTATGCAGGATGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAaag                                                                                      <  2:6584474/101‑1 (MQ=255)
            aaaCAGTTTATGCAGGATGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAaaga                                                                                     <  2:6584695/101‑1 (MQ=255)
             aaCAGTTTATGCAGGATGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAaagaa                                                                                    <  1:6584571/101‑1 (MQ=255)
                  tttATGCAGGATGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCg                                                                               >  1:6585091/1‑101 (MQ=255)
                   ttATGCAGGATGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCgg                                                                              >  2:6585095/1‑101 (MQ=255)
                     aTGCAGGATGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGGGACGCACAATGAGGCAAAGAAGCCCGGcc                                                                            <  1:6584455/101‑1 (MQ=255)
                     aTGCAGGATGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGcc                                                                            >  1:6585096/1‑101 (MQ=255)
                     aTGCAGGATGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGcc                                                                            <  2:6584717/101‑1 (MQ=255)
                       gCAGGATGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCt                                                                          >  1:6585101/1‑101 (MQ=255)
                       gCAGGATGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCt                                                                          >  2:6585103/1‑101 (MQ=255)
                         aGGATGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTAc                                                                        >  1:6585105/1‑101 (MQ=255)
                           gATGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCt                                                                      <  1:6584615/101‑1 (MQ=255)
                           gATGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCt                                                                      <  1:6584582/101‑1 (MQ=255)
                           gATGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCt                                                                      <  2:6584567/101‑1 (MQ=255)
                           gATGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCt                                                                      <  2:6584686/101‑1 (MQ=255)
                            aTGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGCGGCAAAGAAGCCCGGCCCTACCTg                                                                     >  2:6585109/1‑101 (MQ=255)
                            aTGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTg                                                                     <  2:6584597/101‑1 (MQ=255)
                            aTGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTg                                                                     <  2:6584423/101‑1 (MQ=255)
                              ggCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGaa                                                                   >  1:6585112/1‑101 (MQ=255)
                                 ggTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATcc                                                                >  1:6585118/1‑101 (MQ=255)
                                 ggTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATcc                                                                <  1:6584331/101‑1 (MQ=255)
                                 ggTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATcc                                                                >  2:6585116/1‑101 (MQ=255)
                                  gTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGaa                                                                   <  2:6584779/97‑1 (MQ=255)
                                   tAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCg                                                              <  2:6584701/101‑1 (MQ=255)
                                    aaCCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGa                                                             <  2:6584758/101‑1 (MQ=255)
                                    aaCCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGa                                                             >  2:6585123/1‑101 (MQ=255)
                                    aaCCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGa                                                             >  1:6585124/1‑101 (MQ=255)
                                     aCCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGAt                                                            >  2:6585128/1‑101 (MQ=255)
                                     aCCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGAt                                                            <  1:6584801/101‑1 (MQ=255)
                                        aGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGt                                                         >  1:6585132/1‑101 (MQ=255)
                                        aGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGt                                                         <  1:6584752/101‑1 (MQ=255)
                                          accaccGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAAt                                                                                               <  2:6584606/61‑1 (MQ=255)
                                          accaccGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGtat                                                       >  1:6585136/1‑101 (MQ=255)
                                          accaccGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGtat                                                       <  1:6584699/101‑1 (MQ=255)
                                          accaccGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGtat                                                       >  2:6585137/1‑101 (MQ=255)
                                           ccaccGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGtata                                                      <  1:6584490/101‑1 (MQ=255)
                                            caccGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATaa                                                     <  1:6584609/101‑1 (MQ=255)
                                            caccGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATaa                                                     <  2:6584644/101‑1 (MQ=255)
                                             accGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATaaa                                                    >  2:6585143/1‑101 (MQ=255)
                                              ccGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAAc                                                   <  1:6584765/101‑1 (MQ=255)
                                              ccGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAAc                                                   >  1:6585144/1‑101 (MQ=255)
                                                  tgatgaCGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAAc                                               <  1:6584438/101‑1 (MQ=255)
                                                  tgatgaCGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAAc                                               >  2:6585148/1‑101 (MQ=255)
                                                   gatgaCGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACa                                              >  1:6585151/1‑101 (MQ=255)
                                                    atgaCGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACAg                                             <  2:6584243/101‑1 (MQ=255)
                                                     tgaCGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACAgg                                            >  2:6585155/1‑101 (MQ=255)
                                                      gaCGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACAgg                                            <  2:6584661/100‑1 (MQ=255)
                                                       aCGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACAGGAt                                          >  1:6585158/1‑101 (MQ=255)
                                                        cGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACAg                                             >  2:6585163/1‑97 (MQ=255)
                                                        cGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACAGGATc                                         <  1:6584531/101‑1 (MQ=255)
                                                        cGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACAGGATc                                         >  1:6585161/1‑101 (MQ=255)
                                                         gATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACAGGATCa                                        <  2:6584677/101‑1 (MQ=255)
                                                            cAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGATTCCCGATAGTATAAACGAACAGGATCAGGc                                     <  1:6584667/101‑1 (MQ=255)
                                                            cAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACAGGATCAGGc                                     <  1:6584401/101‑1 (MQ=255)
                                                            cAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACAGGATCAGGc                                     >  1:6585168/1‑101 (MQ=255)
                                                              gAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACAGGATCAGGCCa                                   >  2:6585171/1‑101 (MQ=255)
                                                               aaCAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACAGGATCAGGCCaa                                  <  1:6584794/101‑1 (MQ=255)
                                                               aaCAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACAGGATCAGGCCaa                                  <  2:6584805/101‑1 (MQ=255)
                                                               aaCAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACAGGATCAGGCCaa                                  >  1:6585173/1‑101 (MQ=255)
                                                               aaCAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACAGGATCAGGCCaa                                  <  2:6584931/101‑1 (MQ=255)
                                                                aCAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACCaa                                                >  2:6585176/1‑86 (MQ=255)
                                                                 cAGCAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACAGGATCAGGCCAAAt                                <  1:6584817/101‑1 (MQ=255)
                                                                    cAAACAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACAGGATCAGGCCAAATTcc                             <  1:6584646/101‑1 (MQ=255)
                                                                      aaCAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTTCCTGAATCCCGATAGTATAAACGAACAGGAt                                          <  2:6584500/86‑1 (MQ=255)
                                                                        cAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACAGGATCAGGCCAAATTCCTGaa                         >  2:6585181/1‑101 (MQ=255)
                                                                        cAGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACAGGATCAGGCCAAATTCCTGaa                         >  1:6585182/1‑101 (MQ=255)
                                                                         aGGTTGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACAGGATCAGGCCAAATTCCTGAAt                        <  2:6584877/101‑1 (MQ=255)
                                                                             tGAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACAGGATCAGGCCAAATTCCTGAATAACa                    >  2:6585188/1‑101 (MQ=255)
                                                                              gAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACAGGATCAGGCCAAATTCCTGAATAACAt                   <  2:6584896/101‑1 (MQ=255)
                                                                              gAGGCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACAGGATCAGGCCAAATTCCTGAATAACAt                   <  1:6584912/101‑1 (MQ=255)
                                                                                ggCGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACAGGATCAGGCCAAATTCCTGAATAACATgc                 >  1:6585190/1‑101 (MQ=255)
                                                                                  cGTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACAGGATCAGGCCAAATTCCTGAATAACATgcag               <  2:6584829/101‑1 (MQ=255)
                                                                                   gTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACAGGATCAGGCCAAATTCCTGAATAACATgcagc              <  2:6584545/101‑1 (MQ=255)
                                                                                   gTAATCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACAGGATCAGGCCAAATTCCTGAATAACATgcagc              >  2:6585194/1‑101 (MQ=255)
                                                                                     aaTCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACAGGATCAGGCCAAATTCCTGAATAACATGCAGCat            >  1:6585198/1‑101 (MQ=255)
                                                                                     aaTCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACAGGATCAGGCCAAATTCCTGAATAACATGCAGCat            >  2:6585200/1‑101 (MQ=255)
                                                                                      aTCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACAGGATCAGGcc                                    >  2:6585204/1‑76 (MQ=255)
                                                                                       tCCGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACAGGATCAGGCCAAATTCCTGAATAACATGCAGCatat          <  1:6584971/101‑1 (MQ=255)
                                                                                        ccGGAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACAGGATCAGGCCAAATTCCTGAATAACATGCAGCATATc         >  1:6585207/1‑101 (MQ=255)
                                                                                          ggAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACAGGATCAGGCCAAATTCCTGAATAACATGCAGCATATCGc       <  2:6584818/101‑1 (MQ=255)
                                                                                          ggAGACGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACAGGATCAGGCCAAATTCCTGAATAACATGCAGCATATCGc       >  2:6585210/1‑101 (MQ=255)
                                                                                           gagaCGCAAAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAAc                                               <  1:6584548/60‑1 (MQ=255)
                                                                                            agaCGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACAGGATCAGGCCAAATTCCTGAATAACATGCAGCATATCGCTa     >  2:6585214/1‑101 (MQ=255)
                                                                                             gaCGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACAGGATCAGGCCAAATTCCTGAATAACATGCAGCATATCGCTAc    >  1:6585215/1‑101 (MQ=255)
                                                                                              aCGCACAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACAGGATCAGGCCAAATTCCTGAATAACATGCAGCATATCGCTACt   <  2:6584625/101‑1 (MQ=255)
                                                                                                  |                                                                                                 
CAGTGGAATATCAAACAGTTTATGCAGGATGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCGTAATCCGGAGACGCGCAATGAGGCAAAGAAGCCCGGCCCTACCTGAATCCCGATAGTATAAACGAACAGGATCAGGCCAAATTCCTGAATAACATGCAGCATATCGCTACTT  >  NC_000913/3865917‑3866112

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