Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 635,261 A→G 47.9% V39A (GTT→GCT)  ybdM ← Spo0J family protein, ParB‑like nuclease domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913635,2610AG47.9% 22.3 / 198.2 152V39A (GTT→GCT) ybdMSpo0J family protein, ParB‑like nuclease domain
Reads supporting (aligned to +/- strand):  ref base A (49/30);  new base G (38/35);  total (87/65)
Fisher's exact test for biased strand distribution p-value = 2.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 155 total aligned reads displayed.

CGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTT  >  NC_000913/635163‑635360
                                                                                                  |                                                                                                   
cGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAAc                                                                                                    >  2:108558/1‑100 (MQ=255)
 gggTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAAcc                                                                                                   >  1:815082/1‑100 (MQ=255)
  ggtggtGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCt                                                                                                  >  2:1150725/1‑100 (MQ=255)
  ggtggtGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCt                                                                                                  >  1:720958/1‑100 (MQ=255)
     ggtgCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGAt                                                                                               <  2:1710442/100‑1 (MQ=255)
      gtgCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATg                                                                                              <  1:2052172/100‑1 (MQ=255)
       tgCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATgg                                                                                             <  2:2335025/100‑1 (MQ=255)
       tgCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATgg                                                                                             >  1:594873/1‑100 (MQ=255)
          cACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATc                                                                                          >  1:1089478/1‑100 (MQ=255)
          cACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATc                                                                                          >  1:2194275/1‑100 (MQ=255)
          cACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATc                                                                                          <  1:3686461/100‑1 (MQ=255)
            cGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGc                                                                                        <  1:545595/100‑1 (MQ=255)
            cGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGc                                                                                        <  2:2284073/100‑1 (MQ=255)
                gttTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATg                                                                                    >  2:2135818/1‑100 (MQ=255)
                 tttGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGc                                                                                   >  2:4165849/1‑100 (MQ=255)
                  ttGGGTTGTAATCATTAGGCATAAGTTGGCTTTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCg                                                                                  <  2:3092219/100‑1 (MQ=255)
                   tGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCgg                                                                                 <  2:2782383/100‑1 (MQ=255)
                     ggTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGACAGGGGCTAGCCTGATGGCTCGCCCTGCGGAt                                                                               >  1:474531/1‑100 (MQ=255)
                     ggTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGAt                                                                               <  2:1578366/100‑1 (MQ=255)
                     ggTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGAt                                                                               >  2:684463/1‑100 (MQ=255)
                     ggTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGAt                                                                               <  2:1851967/100‑1 (MQ=255)
                     ggTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGAt                                                                               >  1:1167318/1‑100 (MQ=255)
                      gTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCggtt                                                                               >  1:33488/1‑97 (MQ=255)
                      gTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGAtt                                                                              >  2:151579/1‑100 (MQ=255)
                      gTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCACCATGCGGa                                                                                >  2:69869/1‑98 (MQ=255)
                         gTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCa                                                                           <  2:3618258/100‑1 (MQ=255)
                          tAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCAt                                                                          >  1:1164069/1‑100 (MQ=255)
                          tAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCAt                                                                          >  1:2934897/1‑100 (MQ=255)
                             tCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTaa                                                                       >  1:380298/1‑100 (MQ=255)
                              cATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTa                                                                        >  1:4071103/1‑98 (MQ=255)
                                 tAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGc                                                                   >  2:3635975/1‑100 (MQ=255)
                                   ggCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCtt                                                                 <  2:1474708/100‑1 (MQ=255)
                                   ggCATAAGTTGGATGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCtt                                                                 >  2:2812289/1‑100 (MQ=255)
                                    gCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCAACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCttt                                                                >  1:1071934/1‑100 (MQ=255)
                                     cATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGGTGGATCGCCATGCGGATTTCATTAATCGCTTTg                                                               >  1:3998739/1‑100 (MQ=255)
                                     cATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTg                                                               <  2:3158267/100‑1 (MQ=255)
                                        aaGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATg                                                            <  1:3079036/100‑1 (MQ=255)
                                         aGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGc                                                           <  2:2804858/100‑1 (MQ=255)
                                         aGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGc                                                           <  1:692/100‑1 (MQ=255)
                                          gTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCg                                                          <  2:2017340/100‑1 (MQ=255)
                                          gTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCg                                                          <  2:3759163/100‑1 (MQ=255)
                                          gTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCg                                                          >  2:2066915/1‑100 (MQ=255)
                                            tGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGat                                                        <  2:2223280/100‑1 (MQ=255)
                                              gCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGatca                                                      <  2:1304705/100‑1 (MQ=255)
                                              gCTGTTTTTAACCCACAGAACGCAATACACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGatca                                                      >  1:584584/1‑100 (MQ=255)
                                                 gTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCt                                                   <  1:3138826/100‑1 (MQ=255)
                                                  tttttAACCCACAGAACGCAATCCACTGCCTCTTCACGAAAGGCGCTAGCCTGCTGGATCGCCATGCGGATTTCATTAATCGCTTTGCTGCGATCATCtt                                                  >  1:854959/1‑100 (MQ=255)
                                                   ttttAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTc                                                 >  1:4144056/1‑100 (MQ=255)
                                                      naaCCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCgg                                              >  2:3230368/2‑100 (MQ=255)
                                                       accccACAGAACGCAATCCACTGGCTCTGCACGAAAGGGGCTAGTCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGc                                             <  2:1625998/98‑1 (MQ=255)
                                                       aaCCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGc                                             >  1:2108292/1‑100 (MQ=255)
                                                        aCCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCa                                            >  2:2951649/1‑100 (MQ=255)
                                                        aCCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCa                                            >  1:3606849/1‑100 (MQ=255)
                                                        aCCAACATAACGCTATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGCATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGc                                             <  1:3162792/99‑1 (MQ=255)
                                                         cccACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCaa                                           >  2:3852658/1‑100 (MQ=255)
                                                         cccACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCaa                                           >  2:1414666/1‑100 (MQ=255)
                                                         cccACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCaa                                           >  2:1103961/1‑100 (MQ=255)
                                                         cccACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCaa                                           >  2:3300005/1‑100 (MQ=255)
                                                            acaGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCt                                        >  2:563950/1‑100 (MQ=255)
                                                            acaGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCt                                        >  2:570638/1‑100 (MQ=255)
                                                            acaGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCt                                        >  2:2971790/1‑100 (MQ=255)
                                                               gAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGc                                     <  1:4165827/100‑1 (MQ=255)
                                                                aaCGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCg                                    <  1:471227/100‑1 (MQ=255)
                                                                  cGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCgag                                  <  2:3936189/100‑1 (MQ=255)
                                                                   gCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCgaga                                 >  2:2685712/1‑100 (MQ=255)
                                                                   gCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCgagn                                 <  1:2577531/100‑2 (MQ=255)
                                                                    cAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGaa                                >  1:1310650/1‑100 (MQ=255)
                                                                      aTCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAAt                              >  1:652626/1‑100 (MQ=255)
                                                                       tCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAAtt                             <  2:2063374/100‑1 (MQ=255)
                                                                       tCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAAtt                             >  1:1195422/1‑100 (MQ=255)
                                                                        ccACTGGCTCTTCATGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTg                            <  1:2909224/100‑1 (MQ=255)
                                                                        ccACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTg                            <  1:3260327/100‑1 (MQ=255)
                                                                         cACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGc                           >  1:1640775/1‑100 (MQ=255)
                                                                          aCTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCg                          >  1:3068234/1‑100 (MQ=255)
                                                                            tGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGtt                        <  1:1586336/100‑1 (MQ=255)
                                                                            tGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGACATTGCGtt                        >  2:799743/1‑100 (MQ=255)
                                                                             ggCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTa                       <  2:1800848/100‑1 (MQ=255)
                                                                             ggCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTa                       <  2:1181277/100‑1 (MQ=255)
                                                                             ggCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTa                       >  2:3851601/1‑100 (MQ=255)
                                                                             ggCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTa                       >  1:1881364/1‑100 (MQ=255)
                                                                               ctgttCACGAAAGGGGCTAGCCTGATGGATCGCAAGGCGGATTTCATGAAGCGCTTTGATGCGATCAGCTTCCGGCAAGCTGGCGAGAAATTGCGtt                        <  1:1777380/94‑1 (MQ=255)
                                                                               ctctTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTaaa                     >  2:364346/1‑100 (MQ=255)
                                                                                tctTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAAt                    <  2:3731290/100‑1 (MQ=255)
                                                                                 ctTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATc                   >  2:908573/1‑100 (MQ=255)
                                                                                  ttCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATcc                  >  2:792885/1‑100 (MQ=255)
                                                                                  ttCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATcc                  >  2:2525741/1‑100 (MQ=255)
                                                                                   tCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATcc                  <  1:120420/99‑1 (MQ=255)
                                                                                    cACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTg                <  2:1195425/100‑1 (MQ=255)
                                                                                    cACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTg                >  1:3279426/1‑100 (MQ=255)
                                                                                     aCGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGc               <  1:3457625/100‑1 (MQ=255)
                                                                                      cGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCg              >  1:2277030/1‑100 (MQ=255)
                                                                                      cGAAAGGGGCTAACCTGATGGATTGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCg              <  1:581514/100‑1 (MQ=255)
                                                                                        aaaGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGtt            >  1:2377889/1‑100 (MQ=255)
                                                                                          aGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTaa          <  1:1971452/100‑1 (MQ=255)
                                                                                            gggCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATc        <  2:4215238/100‑1 (MQ=255)
                                                                                              gCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCgt      <  1:1324615/100‑1 (MQ=255)
                                                                                               cTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCgtt     >  2:278335/1‑100 (MQ=255)
                                                                                                tAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCgttg    <  1:2436990/100‑1 (MQ=255)
                                                                                                  aCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCgttgtt  >  2:1931124/1‑100 (MQ=255)
                                                                                                  aCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCgttgtt  <  2:2432811/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
CGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTT  >  NC_000913/635163‑635360

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