Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,008,491 2 bp→CG intergenic (+58/‑174) b4555 → / → yigM hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,008,4910GC100.0% 567.8 / NA 159intergenic (+58/‑175)b4555/yigMhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (74/85);  total (74/85)
*NC_0009134,008,4920CG99.4% 561.6 / NA 158intergenic (+59/‑174)b4555/yigMhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base G (72/85);  total (73/85)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 167 total aligned reads displayed.

AGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGAAAGCTGCATTATGGGCAGTGGCCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACT  >  NC_000913/4008409‑4008577
                                                                                        ||                                                                                       
gccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCgcacc                                                                                     >  1:791574/3‑93 (MQ=255)
   tGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCCGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGc                                                                                  >  1:4251499/1‑93 (MQ=255)
   tGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAgg                                                                                  >  1:4227871/1‑92 (MQ=255)
   tGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGc                                                                                  <  2:4328641/93‑1 (MQ=255)
   tGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGc                                                                                  >  1:3585612/1‑93 (MQ=255)
   tGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGc                                                                                  >  1:956820/1‑93 (MQ=255)
   tGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGc                                                                                  <  1:66123/93‑1 (MQ=255)
    gAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCg                                                                                 <  1:1641550/93‑1 (MQ=255)
    gAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCg                                                                                 >  2:1296504/1‑93 (MQ=255)
    gAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCg                                                                                 >  1:5755649/1‑93 (MQ=255)
       ggATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCt                                                                              <  1:5916574/93‑1 (MQ=255)
                    tttggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCAt                                                                 <  2:1727595/89‑1 (MQ=255)
                      tggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATcc                                                               <  2:704227/91‑1 (MQ=255)
                      tggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATcc                                                               <  2:704226/91‑1 (MQ=255)
                      tggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATcc                                                               <  1:1618191/91‑1 (MQ=255)
                      tggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATcc                                                               <  2:4971818/91‑1 (MQ=255)
                          gACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAg                                                           >  2:2100642/1‑93 (MQ=255)
                            cGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCt                                                         <  2:5184845/93‑1 (MQ=255)
                             gggATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTg                                                       >  2:586952/1‑93 (MQ=255)
                               gATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGcggga                                                      >  1:2704025/1‑89 (MQ=255)
                               gATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTgga                                                     >  2:2900952/1‑91 (MQ=255)
                                 tGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTggaag                                                   <  2:1949765/93‑5 (MQ=255)
                                 tGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTggaag                                                   <  1:3735156/93‑5 (MQ=255)
                                   aaCGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTg                                                 >  2:264074/1‑93 (MQ=255)
                                   aaCGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTg                                                 <  1:3115705/93‑1 (MQ=255)
                                   aaCGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTg                                                 <  2:3137003/93‑1 (MQ=255)
                                    acgacgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGa                                                >  2:3514420/1‑93 (MQ=255)
                                    acgacgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGa                                                >  1:632002/1‑93 (MQ=255)
                                    acgacgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGa                                                <  2:5551018/93‑1 (MQ=255)
                                      gacgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAtt                                              <  1:3436527/93‑1 (MQ=255)
                                      gacgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAtt                                              <  1:1384626/93‑1 (MQ=255)
                                      gacgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAtt                                              <  2:2727450/93‑1 (MQ=255)
                                       acgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTg                                             >  1:5366060/1‑93 (MQ=255)
                                       acgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTg                                             >  2:5156872/1‑93 (MQ=255)
                                       acgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTg                                             <  2:2245207/93‑1 (MQ=255)
                                        cgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCCCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTgg                                            >  1:1516147/1‑93 (MQ=255)
                                        cgCCGAAATGCTGTAGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTgg                                            >  2:3286360/1‑93 (MQ=255)
                                         gCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGt                                           <  1:4745399/93‑1 (MQ=255)
                                           cGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTAc                                         <  1:3296422/93‑1 (MQ=255)
                                            gAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGGAGTGATTGGTACt                                        >  2:2659482/1‑93 (MQ=255)
                                            gAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACt                                        >  2:5539366/1‑93 (MQ=255)
                                             aaaTGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTa                                       >  2:4710536/1‑93 (MQ=255)
                                             aaaTGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGATGGAAGTGATTGGTACTa                                       <  1:265504/93‑1 (MQ=255)
                                              aaTGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTaa                                      <  1:5666223/93‑1 (MQ=255)
                                              aaTGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTaa                                      <  2:421874/93‑1 (MQ=255)
                                              aaTGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTaa                                      <  2:421873/93‑1 (MQ=255)
                                              aaTGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTaa                                      <  2:5288156/93‑1 (MQ=255)
                                              aaTGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTaa                                      <  1:3216533/93‑1 (MQ=255)
                                                tGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATg                                    >  1:1075907/1‑93 (MQ=255)
                                                tGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATg                                    >  2:4451836/1‑93 (MQ=255)
                                                  cTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGcc                                  <  1:4431194/93‑1 (MQ=255)
                                                  cTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGcc                                  >  1:3893941/1‑93 (MQ=255)
                                                  cTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGcc                                  >  2:326993/1‑93 (MQ=255)
                                                   tGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcg                                 >  2:1670387/1‑93 (MQ=255)
                                                   tGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcg                                 >  2:4618886/1‑93 (MQ=255)
                                                    ttCGATGGTGGGGAAAGGCTGCATGATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcga                                <  2:4922589/92‑1 (MQ=255)
                                                    gTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcga                                <  2:776132/93‑1 (MQ=255)
                                                     tCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgac                               >  2:1778771/1‑93 (MQ=255)
                                                     tCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgac                               <  2:2862580/93‑1 (MQ=255)
                                                      cGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACATTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgacg                              <  1:1943809/93‑1 (MQ=255)
                                                       gATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgacga                             >  2:896219/1‑93 (MQ=255)
                                                       gATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgacga                             >  2:1246233/1‑93 (MQ=255)
                                                           ggcggggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGTAAGTGATTGGTACTAATGCCGACGAcgc                          <  1:5054425/89‑1 (MQ=255)
                                                            gcggggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGCAGTGATTGGTACTAATGCCGACGAcgcg                         >  1:437965/4‑93 (MQ=255)
                                                             cggggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCgg                        <  2:4466181/91‑1 (MQ=255)
                                                             cggggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCgg                        <  1:5861390/91‑1 (MQ=255)
                                                              ggggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGt                       <  2:2260528/92‑1 (MQ=255)
                                                               gggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTg                      >  1:3665109/1‑93 (MQ=255)
                                                               gggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTg                      >  2:3550725/1‑93 (MQ=255)
                                                               gggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTg                      >  2:4813301/1‑93 (MQ=255)
                                                                ggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGc                     <  1:2351478/93‑1 (MQ=255)
                                                                   aaggCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGc                  >  2:3725327/4‑93 (MQ=255)
                                                                   aaggCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGCGATTGGTACTAATGCCGACGACGCGGTGCCGc                  >  1:5460017/4‑93 (MQ=255)
                                                                    aggCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCa                 >  1:4346356/3‑93 (MQ=255)
                                                                    aggCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCa                 <  1:1183419/91‑1 (MQ=255)
                                                                     ggCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCAt                <  2:5841925/92‑1 (MQ=255)
                                                                      gCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCAtt               >  1:2628476/1‑93 (MQ=255)
                                                                       cTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTa              <  2:378179/93‑1 (MQ=255)
                                                                        tGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGCGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTAt             >  2:5766905/1‑93 (MQ=255)
                                                                        tGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTAt             >  1:50387/1‑93 (MQ=255)
                                                                         gCATTATGGGCAGTGCGCACCAGCTTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATc            <  1:1657154/93‑1 (MQ=255)
                                                                         gCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATc            <  2:383212/93‑1 (MQ=255)
                                                                         gCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATc            <  1:573060/93‑1 (MQ=255)
                                                                         gCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATc            >  1:1989713/1‑93 (MQ=255)
                                                                         gCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATc            <  1:1804050/93‑1 (MQ=255)
                                                                         gCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATc            <  1:1054860/93‑1 (MQ=255)
                                                                          cATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCt           <  1:4976465/93‑1 (MQ=255)
                                                                          cATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCt           >  1:940281/1‑93 (MQ=255)
                                                                          cATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCt           <  2:199396/93‑1 (MQ=255)
                                                                            ttATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGc         <  1:2562207/93‑1 (MQ=255)
                                                                            ttATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGc         <  2:1325476/93‑1 (MQ=255)
                                                                              aTGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGt       >  2:4701548/1‑93 (MQ=255)
                                                                              aTGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGt       <  2:3767761/93‑1 (MQ=255)
                                                                              aTGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGt       <  2:1459593/93‑1 (MQ=255)
                                                                              aTGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGt       >  2:2884237/1‑93 (MQ=255)
                                                                               tGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTa      >  2:6102919/1‑93 (MQ=255)
                                                                                gggCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTaa     <  1:4710266/93‑1 (MQ=255)
                                                                                gggCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTaa     <  2:3571405/93‑1 (MQ=255)
                                                                                  gCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAGAc   >  1:4288826/1‑93 (MQ=255)
                                                                                  gCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAAc   <  2:4992451/93‑1 (MQ=255)
                                                                                        ||                                                                                       
AGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGAAAGCTGCATTATGGGCAGTGGCCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACT  >  NC_000913/4008409‑4008577

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