Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,056,373 A→G K129K (AAA→AAG)  yihK → putative GTP‑binding factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,056,3730AG98.3% 556.5 / ‑1.0 177K129K (AAA→AAG) yihKputative GTP‑binding factor
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (69/105);  total (70/107)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

Only 100 of 185 total aligned reads displayed.

TGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAA‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGG  >  NC_000913/4056283‑4056461
                                                                                          |                                                                                         
 gACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑Gtt                                                                                       >  2:4547165/1‑93 (MQ=255)
 gACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑Gtt                                                                                       <  2:2889041/93‑1 (MQ=255)
 gACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑Gtt                                                                                       <  2:3013247/93‑1 (MQ=255)
 gACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑Gtt                                                                                       <  2:3119165/93‑1 (MQ=255)
  aCGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTg                                                                                      <  1:3730007/93‑1 (MQ=255)
   cGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGa                                                                                     <  2:3914028/93‑1 (MQ=255)
    gCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGAc                                                                                    >  1:2409579/1‑93 (MQ=255)
    gCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGAc                                                                                    <  2:1338418/93‑1 (MQ=255)
    gCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGAc                                                                                    >  1:1094029/1‑93 (MQ=255)
    gCAGTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGAc                                                                                    <  1:3582153/93‑1 (MQ=255)
      aTTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCg                                                                                  <  2:3477400/93‑1 (MQ=255)
       tttGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGc                                                                                 >  2:4435986/1‑93 (MQ=255)
         tGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGccc                                                                               >  1:4411695/1‑93 (MQ=255)
          gACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCt                                                                              <  1:1725792/93‑1 (MQ=255)
           aCGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAGCAAG‑GTTGACCGCCCTg                                                                             >  2:2075788/1‑93 (MQ=255)
           aCGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTg                                                                             <  2:3578381/93‑1 (MQ=255)
           aCGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTg                                                                             <  1:2752998/93‑1 (MQ=255)
           aCGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTAT‑AACAAAGGTTGACCGCCCTg                                                                             >  1:3984918/1‑93 (MQ=255)
                ccGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGgcgcg                                                                        <  2:2607560/93‑1 (MQ=255)
                 cGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGgcgcgc                                                                       >  2:2421259/1‑93 (MQ=255)
                  gATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGcgcgcg                                                                      >  1:4436476/1‑93 (MQ=255)
                  gATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGcgcgcg                                                                      >  2:3535219/1‑93 (MQ=255)
                   aTGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGt                                                                     <  2:3949411/93‑1 (MQ=255)
                   aTGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGt                                                                     >  1:553038/1‑93 (MQ=255)
                    tGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTc                                                                    <  1:2346507/93‑1 (MQ=255)
                    tGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTc                                                                    <  1:2222749/93‑1 (MQ=255)
                        gCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGa                                                                <  2:1118065/93‑1 (MQ=255)
                         cAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGAt                                                               <  2:3911588/93‑1 (MQ=255)
                         cAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGGTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGAt                                                               >  2:1918231/1‑93 (MQ=255)
                          aaaCGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGGTATCAACAAG‑GTTGACCGCCCCGGCCCGCGCCCTGAtg                                                              >  1:1017936/1‑92 (MQ=255)
                           aaCGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTg                                                             <  1:1044905/93‑1 (MQ=255)
                           aaCGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTg                                                             <  1:627105/93‑1 (MQ=255)
                           aaCGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTg                                                             >  1:4233150/1‑93 (MQ=255)
                             cgcgCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTggg                                                           >  2:3786329/1‑93 (MQ=255)
                              gcgcTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGt                                                          >  1:3833907/1‑93 (MQ=255)
                              gcgcTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGt                                                          <  2:2025027/93‑1 (MQ=255)
                               cgcTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGtt                                                         <  2:941260/93‑1 (MQ=255)
                               cgcTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGtt                                                         <  1:1672707/93‑1 (MQ=255)
                                gcTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTtg                                                        >  2:1960842/1‑93 (MQ=255)
                                 cTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTtgt                                                       >  1:29705/1‑93 (MQ=255)
                                 cTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTtgt                                                       <  2:1179671/93‑1 (MQ=255)
                                  ttCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTtgtg                                                      >  2:3258160/1‑93 (MQ=255)
                                    cGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGTCCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGa                                                    <  1:688634/93‑1 (MQ=255)
                                        aCCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAg                                                <  1:3547220/93‑1 (MQ=255)
                                        aCCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAg                                                >  1:3319811/1‑93 (MQ=255)
                                        aCCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCACCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAg                                                <  1:3597796/93‑1 (MQ=255)
                                         ccAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATGAgg                                               >  2:2575346/1‑93 (MQ=255)
                                         ccAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAgg                                               <  1:3181992/93‑1 (MQ=255)
                                         ccAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAgg                                               <  2:3038922/93‑1 (MQ=255)
                                           aaaaaaGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGGGCGTCCTGATTGTGTTGTGGGTCAggca                                             >  1:3208419/1‑91 (MQ=255)
                                            aaaaaGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTAt                                            >  2:2511296/1‑93 (MQ=255)
                                             aaaaGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTAtt                                           <  1:3664190/93‑1 (MQ=255)
                                             aaaaGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTAtt                                           <  2:420353/93‑1 (MQ=255)
                                             aaaaGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTAtt                                           <  1:1946348/93‑1 (MQ=255)
                                             aaaaGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTAtt                                           <  1:4636171/93‑1 (MQ=255)
                                              aaaGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTc                                          <  1:50104/93‑1 (MQ=255)
                                              aaaGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTc                                          <  2:4616102/93‑1 (MQ=255)
                                               aaGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCg                                         <  2:1396904/93‑1 (MQ=255)
                                                aGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGa                                        <  1:3321642/93‑1 (MQ=255)
                                                 gCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGAt                                       >  2:2687343/1‑93 (MQ=255)
                                                 gCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGAt                                       >  2:1451794/1‑93 (MQ=255)
                                                  cGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATc                                      <  1:1262869/93‑1 (MQ=255)
                                                  cGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATc                                      <  2:3927755/93‑1 (MQ=255)
                                                   gTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCt                                     >  2:2656960/1‑93 (MQ=255)
                                                     ttGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGt                                   >  2:626663/1‑93 (MQ=255)
                                                      tGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGtt                                  <  2:4261738/93‑1 (MQ=255)
                                                        cTTACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCg                                <  1:3321737/93‑1 (MQ=255)
                                                         ttACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGt                               <  1:279995/93‑1 (MQ=255)
                                                          tACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGtt                              <  1:2712476/93‑1 (MQ=255)
                                                          tACGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGtt                              <  2:3756893/93‑1 (MQ=255)
                                                           aCGGCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTa                             >  1:3006293/1‑93 (MQ=255)
                                                             ggCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAAc                           >  1:3765736/1‑93 (MQ=255)
                                                              gCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAAcc                          <  2:3011323/93‑1 (MQ=255)
                                                              gCCTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAAcc                          >  1:2320529/1‑93 (MQ=255)
                                                               gcTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCt                         <  2:1017996/92‑1 (MQ=255)
                                                                cTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTc                        >  2:1134655/1‑93 (MQ=255)
                                                                cTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTc                        >  2:1540351/1‑93 (MQ=255)
                                                                cTGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTc                        >  2:170400/1‑93 (MQ=255)
                                                                 tGAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCg                       >  1:2205207/1‑93 (MQ=255)
                                                                  gAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGa                      <  2:1289537/93‑1 (MQ=255)
                                                                  gAAGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATATGTTCGTTAACCCCGa                      >  1:2568686/1‑93 (MQ=255)
                                                                    aGCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGAcg                    >  1:2379542/1‑93 (MQ=255)
                                                                     gCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGAcgc                   <  1:4307452/93‑1 (MQ=255)
                                                                     gCCGATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGAcgc                   <  2:4113087/93‑1 (MQ=255)
                                                                        gATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGcgac                <  2:1322037/93‑1 (MQ=255)
                                                                        gATTGTTGTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGcgac                <  1:1789816/93‑1 (MQ=255)
                                                                          ttgttgTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGcgaccg              <  1:2712473/93‑1 (MQ=255)
                                                                          ttgttgTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGcgaccg              >  1:4347032/1‑93 (MQ=255)
                                                                          ttgttgTTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGcgaccg              >  2:2922931/1‑93 (MQ=255)
                                                                           tgttgtTATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGcgaccga             <  2:5683/93‑1 (MQ=255)
                                                                            gttgttATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGcgaccgag            >  1:213312/1‑92 (MQ=255)
                                                                            gttgttATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGcgaccgac            >  1:3472508/1‑93 (MQ=255)
                                                                             ttgttATCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACcgacg           <  1:1520724/93‑1 (MQ=255)
                                                                                   tCAACAAG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGagcagc     <  2:4640166/93‑1 (MQ=255)
                                                                                    caacaaG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCt    >  2:89841/1‑93 (MQ=255)
                                                                                    caacaaG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCt    <  2:918242/93‑1 (MQ=255)
                                                                                     aacaaG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTg   <  1:1613048/93‑1 (MQ=255)
                                                                                     aacaaG‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTg   <  2:2625539/93‑1 (MQ=255)
                                                                                     aacaaA‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTg   <  1:4136924/93‑1 (MQ=255)
                                                                                      acaaG‑GTTGACCGCCCTGGCGCGCATCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTgg  <  2:3235678/93‑1 (MQ=255)
                                                                                          |                                                                                         
TGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAA‑GTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGG  >  NC_000913/4056283‑4056461

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