Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,927,062 +G coding (2/1284 nt)
coding (169/183 nt)
yieM ←
b4564 ←
hypothetical protein
hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,927,0621.G99.1% 787.6 / ‑1.3 222coding (2/1284 nt)
coding (169/183 nt)
yieM
b4564
hypothetical protein
hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base G (100/120);  total (100/122)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

Only 100 of 252 total aligned reads displayed.

ACTGTGGCGTGGAGAGGAGCTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCACGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAAAGAAGACCGCCAGCTGCGGTGAGG  >  NC_000913/3926970‑3927148
                                                                                              |                                                                                      
aCTGTGGCGTGGAGAGGAGCTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCa                                                                                         >  1:3728565/1‑93 (MQ=255)
aCTGTGGCGTGGAGAGGAGCTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCa                                                                                         >  1:414481/1‑93 (MQ=255)
aCTGTGGCGTGGAGAGGAGCTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCa                                                                                         <  2:144042/93‑1 (MQ=255)
aCTGTGGCGTGGAGAGGAGCTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCa                                                                                         <  2:3974229/93‑1 (MQ=255)
  tgtgGCGTGGAGAGGAGCTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCagc                                                                                       >  2:739687/1‑91 (MQ=255)
  tgtgGCGTGGAGAGGAGCTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCagc                                                                                       <  1:3247750/93‑3 (MQ=255)
  tgtgGCGTGGAGAGGAGCTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCagc                                                                                       <  1:1454769/93‑3 (MQ=255)
  tgtgGCGTGGAGAGGAGCTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCagc                                                                                       <  1:1403796/93‑3 (MQ=255)
  tgtgGCGTGGAGAGGAGCTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCagc                                                                                       >  1:1391799/1‑91 (MQ=255)
   gtgGCGTGGAGAGGTGCTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCagcg                                                                                      >  1:4027102/1‑90 (MQ=255)
   gtgGCGTGGAGAGGAGCTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCagcg                                                                                      >  1:399477/1‑90 (MQ=255)
       ggTGGAGAGGAGCTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTc                                                                                  <  1:5284343/92‑1 (MQ=255)
       cGTGGAGAGGAGCTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTc                                                                                  >  2:1521527/1‑93 (MQ=255)
        gTGGAGAGGAGCTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCa                                                                                 <  2:5318438/93‑1 (MQ=255)
        gTGGAGAGGAGCTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCa                                                                                 >  2:5441249/1‑93 (MQ=255)
        gTGGAGAGGAGCTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCa                                                                                 <  1:3378611/93‑1 (MQ=255)
         tGGAGAGGAGCTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCAc                                                                                >  2:3751788/1‑93 (MQ=255)
         tGGAGAGGAGCTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCAc                                                                                >  2:203768/1‑93 (MQ=255)
         tGGAGAGGAGCTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCAc                                                                                >  1:5164383/1‑93 (MQ=255)
           gagaGGAGCTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGc                                                                              >  1:2740446/1‑93 (MQ=255)
             gaggagCTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCa                                                                            <  2:1324219/93‑1 (MQ=255)
                gagCTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCg                                                                         <  1:5002912/93‑1 (MQ=255)
                 agCTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCgg                                                                        <  2:4480676/93‑1 (MQ=255)
                 agCTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCgg                                                                        <  1:3564492/93‑1 (MQ=255)
                  gCTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCggg                                                                       <  1:2108027/93‑1 (MQ=255)
                   cTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCGAGCgggg                                                                      >  2:3115570/1‑93 (MQ=255)
                   cTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCgggg                                                                      <  2:3141987/93‑1 (MQ=255)
                    tGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGa                                                                     >  2:5554585/1‑93 (MQ=255)
                    tGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGa                                                                     >  2:3203769/1‑93 (MQ=255)
                      gCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAAc                                                                   <  1:232479/93‑1 (MQ=255)
                       gtttGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAAca                                                                  <  2:4275938/92‑1 (MQ=255)
                       cTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAAca                                                                  >  2:2792158/1‑93 (MQ=255)
                       cTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAAca                                                                  <  2:4868975/93‑1 (MQ=255)
                         ttGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAAcatc                                                                >  2:3784023/1‑93 (MQ=255)
                          tGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAAcatca                                                               <  2:182772/93‑1 (MQ=255)
                            cTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACatcatc                                                             <  1:4183529/93‑1 (MQ=255)
                             taaTAGCAACTCCCCTCTTCGGTGCGTTCCGGCGGGCCTCAGGCATCATAAAGCCGACTGCGCAGCGATTCACGCCAGCGGGGAACAtcatca                                                            <  1:5453631/91‑1 (MQ=255)
                             tGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACAtcatca                                                            >  2:2852906/1‑93 (MQ=255)
                             tGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACAtcatca                                                            <  1:4867566/93‑1 (MQ=255)
                                tAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTg                                                         >  1:6120149/1‑93 (MQ=255)
                                 aGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGa                                                        <  2:599440/93‑1 (MQ=255)
                                  gCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGAt                                                       <  1:5219361/93‑1 (MQ=255)
                                     caTCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTgc                                                    >  2:1140869/1‑93 (MQ=255)
                                     caTCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTgc                                                    >  2:5921614/1‑93 (MQ=255)
                                     caTCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACCGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTgc                                                    >  2:650688/1‑93 (MQ=255)
                                      aTCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTgcc                                                   <  2:5505301/93‑1 (MQ=255)
                                      aTCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTgcc                                                   <  2:181253/93‑1 (MQ=255)
                                       tCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTgccg                                                  >  2:5168975/1‑93 (MQ=255)
                                          ccTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCt                                               >  2:5208015/1‑93 (MQ=255)
                                           ctctTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCtt                                              <  2:2426199/93‑1 (MQ=255)
                                            tctACGGTGAGTTCCGGCGTGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTc                                             <  1:5960441/93‑1 (MQ=255)
                                             ctTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCa                                            >  1:1743877/1‑93 (MQ=255)
                                             ctTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCa                                            <  1:1689101/93‑1 (MQ=255)
                                             ctTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCa                                            <  1:1114653/93‑1 (MQ=255)
                                              ttCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAg                                           <  2:5290602/93‑1 (MQ=255)
                                                 ggTGAGTTCCGGCGGGTCTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCg                                        <  1:920614/93‑1 (MQ=255)
                                                  gTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCgg                                       >  2:5339003/1‑92 (MQ=255)
                                                  gTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGt                                       <  2:167269/93‑1 (MQ=255)
                                                    gAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTgg                                     <  2:4353384/93‑1 (MQ=255)
                                                      gTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGa                                   >  1:5837835/1‑93 (MQ=255)
                                                       ttCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGaa                                  >  1:2213397/1‑93 (MQ=255)
                                                       ttCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGaa                                  <  1:14446/93‑1 (MQ=255)
                                                       ttCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGaa                                  >  1:5714959/1‑93 (MQ=255)
                                                            gcggGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAAttttt                             >  1:2382293/1‑93 (MQ=255)
                                                             cggGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTc                            >  1:3178983/1‑93 (MQ=255)
                                                             cggGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTc                            <  2:572161/93‑1 (MQ=255)
                                                              gggACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCa                           >  1:3807726/1‑93 (MQ=255)
                                                               ggACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCaa                          <  2:2168107/93‑1 (MQ=255)
                                                                gACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCaaa                         >  2:4959879/1‑93 (MQ=255)
                                                                   tCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCGGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAaagaa                      >  1:5526039/1‑93 (MQ=255)
                                                                   tCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAaagaa                      <  2:5408193/93‑1 (MQ=255)
                                                                    cGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAaagaag                     <  2:3924582/93‑1 (MQ=255)
                                                                     gggCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAAagaaga                    <  2:3561102/93‑1 (MQ=255)
                                                                     gggCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAAagaaga                    <  2:5665573/93‑1 (MQ=255)
                                                                     gggCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAAagaaga                    <  2:264723/93‑1 (MQ=255)
                                                                      ggCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAAAGAAGAc                   >  1:3445437/1‑93 (MQ=255)
                                                                      ggCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAAAGAAGAc                   <  1:689877/93‑1 (MQ=255)
                                                                       gCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAAAGAAGAcc                  <  1:2806956/93‑1 (MQ=255)
                                                                       gCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAAAGAAGAcc                  <  2:4257858/93‑1 (MQ=255)
                                                                        cATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAAAGAAGACCg                 <  1:3275067/93‑1 (MQ=255)
                                                                         aTCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAAAGAAGACCGc                >  1:1064324/1‑93 (MQ=255)
                                                                         aTCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAAAGAAGACCGc                <  2:3337935/93‑1 (MQ=255)
                                                                         aTCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAAAGAAGACCGc                >  2:4797031/1‑93 (MQ=255)
                                                                          tCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAAAGAAGACCGcc               >  1:2715090/1‑93 (MQ=255)
                                                                           cTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAAAGAAGACCGCCa              <  2:204772/93‑1 (MQ=255)
                                                                            tttCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAAAGAAGACCGCCAg             <  2:4011012/93‑1 (MQ=255)
                                                                                aGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAAAGAAGACCGCCAGCTGc         >  2:3421260/1‑93 (MQ=255)
                                                                                 gCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAAAGAAGACCGCCAGCTGCg        >  1:4768838/1‑93 (MQ=255)
                                                                                 gCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAAAGAAGACCGCCAGCTGCg        >  1:5636426/1‑93 (MQ=255)
                                                                                   cGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAAAGAAGACCGCCAGCTGCggg      >  1:4358156/1‑92 (MQ=255)
                                                                                   cGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAAAGAAGACCGCCAGCTGCGGt      <  2:5603616/93‑1 (MQ=255)
                                                                                    gACTGCGCGGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAAAGAAGACCGCCAGCTGCGGTg     <  1:5474097/93‑1 (MQ=255)
                                                                                     aCTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAAAGAAGACCGCCAGCTGCGGTGa    <  2:2137323/93‑1 (MQ=255)
                                                                                     aCTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAAAGAAGACCGCCAGCTGCGGTGa    <  1:2470524/93‑1 (MQ=255)
                                                                                     aCTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAAAGAAGACCGCCAGCTGCGGTGa    >  1:1884926/1‑93 (MQ=255)
                                                                                      cTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAAAGAAGACCGCCAGCTGCGGTGAg   >  1:2422094/1‑93 (MQ=255)
                                                                                      cTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAAAGAAGACCGCCAGCTGCGGTGAg   <  1:199863/93‑1 (MQ=255)
                                                                                      cTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAAAGAAGACCGCCAGCTGCGGTGAg   <  1:4622141/93‑1 (MQ=255)
                                                                                       tGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAAAGAAGACCGCCAGCTGCGGTGAgg  >  2:4703661/1‑93 (MQ=255)
                                                                                       tGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAAAGAAGACCGCCAGCTGCGGTGAgc  >  2:679938/1‑92 (MQ=255)
                                                                                              |                                                                                      
ACTGTGGCGTGGAGAGGAGCTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCACGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAAAGAAGACCGCCAGCTGCGGTGAGG  >  NC_000913/3926970‑3927148

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