Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 635,261 A→G 73.8% 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,2610AG73.8% 244.1 / 92.2 145V39A (GTT→GCT) ybdMSpo0J family protein, ParB‑like nuclease domain
Reads supporting (aligned to +/- strand):  ref base A (12/26);  new base G (54/53);  total (66/79)
Fisher's exact test for biased strand distribution p-value = 5.78e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.33e-01

Only 100 of 146 total aligned reads displayed.

GGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGT  >  NC_000913/635164‑635359
                                                                                                 |                                                                                                  
   tggtgCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGa                                                                                               <  2:467527/100‑1 (MQ=255)
   tggtgCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGa                                                                                               <  2:3768200/100‑1 (MQ=255)
    ggtgCCACGTTGTTTGGGTTGTAATCATTTGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGAt                                                                                              <  2:2499682/100‑1 (MQ=255)
    ggtgCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGAt                                                                                              >  2:382762/1‑100 (MQ=255)
     gtgCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATg                                                                                             <  2:1209801/100‑1 (MQ=255)
     gtgCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATg                                                                                             >  1:1781563/1‑100 (MQ=255)
     gtgCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATg                                                                                             <  1:2485875/100‑1 (MQ=255)
       gCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGa                                                                                           >  2:3664252/1‑100 (MQ=255)
          aCGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCg                                                                                        <  1:158442/100‑1 (MQ=255)
           cGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGc                                                                                       <  1:2860235/100‑1 (MQ=255)
             ttgttTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCa                                                                                     <  1:840316/100‑1 (MQ=255)
               gttTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATg                                                                                   <  1:3828468/100‑1 (MQ=255)
               gttTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATg                                                                                   <  2:1533736/100‑1 (MQ=255)
                tttGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGc                                                                                  >  1:2421477/1‑100 (MQ=255)
                 ttGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCTCGAAAGGGGCTAACCTGATGGATCGCCATGCg                                                                                 >  2:3045320/1‑100 (MQ=255)
                 ttGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCg                                                                                 <  1:1021489/100‑1 (MQ=255)
                 ttGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCg                                                                                 <  1:605378/100‑1 (MQ=255)
                 ttGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCg                                                                                 <  2:182045/100‑1 (MQ=255)
                    ggTTGTAATCATTAGGCATAAGTTGGCTGTTTTTACCCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGAt                                                                              <  1:835221/100‑1 (MQ=255)
                     gTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGAtt                                                                             <  1:4034871/100‑1 (MQ=255)
                     gTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGAtt                                                                             <  1:1792902/100‑1 (MQ=255)
                       tGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTc                                                                           <  1:3256497/100‑1 (MQ=255)
                       tGTAATCAGTAGGCATGAGTTGGCTGTTTTTAGCCCACAGAGCGCAAGACTCTGGTTCTTCACGAAAGGAGCTAGCTTGATGGATCGCCATGCGGATTTc                                                                           <  2:846149/100‑1 (MQ=255)
                         tAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCAt                                                                         <  1:2433120/100‑1 (MQ=255)
                          aaTCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCAtt                                                                        >  2:241352/1‑100 (MQ=255)
                          aaTCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCAtt                                                                        <  1:2950338/100‑1 (MQ=255)
                             cATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAAt                                                                     >  1:1998248/1‑100 (MQ=255)
                              aTTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAAt                                                                     <  2:2060220/99‑1 (MQ=255)
                               ttAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCg                                                                   <  2:3095801/100‑1 (MQ=255)
                                tAGGCATAAGTTGGCTGTTTGTAATCCAATGAGTGCGATCCGCTGGCTCTTCGCGAAAGGGGCTAGGCTGATGGATCGCCATGCGGATTTCAGTAATCg                                                                   <  2:3502061/99‑1 (MQ=255)
                                  ggCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCtt                                                                >  1:3501967/1‑100 (MQ=255)
                                  ggCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCtt                                                                >  2:136942/1‑100 (MQ=255)
                                  ggCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCtt                                                                >  1:3607097/1‑100 (MQ=255)
                                  ggCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCtt                                                                >  2:566335/1‑100 (MQ=255)
                                      tAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGAt                                                            <  1:2252518/100‑1 (MQ=255)
                                       aaGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATg                                                           <  2:3890459/100‑1 (MQ=255)
                                        aGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGc                                                          <  1:818814/100‑1 (MQ=255)
                                        aGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGc                                                          >  2:1638253/1‑100 (MQ=255)
                                         gTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCg                                                         >  1:307286/1‑100 (MQ=255)
                                         gTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCg                                                         >  1:4088704/1‑100 (MQ=255)
                                           tGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGat                                                       <  1:3529451/100‑1 (MQ=255)
                                           tGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGat                                                       <  1:3708630/100‑1 (MQ=255)
                                           tGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTTATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGat                                                       <  2:3659482/100‑1 (MQ=255)
                                            ggCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGatc                                                      >  2:1444393/1‑100 (MQ=255)
                                            ggCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGatc                                                      >  2:95778/1‑100 (MQ=255)
                                             gCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGatca                                                     <  1:3451279/100‑1 (MQ=255)
                                             gCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGatca                                                     >  1:26368/1‑100 (MQ=255)
                                             gCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGatca                                                     <  2:1440822/100‑1 (MQ=255)
                                              cTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGatcat                                                    <  1:2097723/100‑1 (MQ=255)
                                               tGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGatcatc                                                   >  1:1350630/1‑100 (MQ=255)
                                               tGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGatcatc                                                   >  2:747648/1‑100 (MQ=255)
                                                  ttttAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTc                                                >  2:4004796/1‑100 (MQ=255)
                                                   tttAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTcc                                               <  1:2642871/100‑1 (MQ=255)
                                                   tttAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTcc                                               <  1:2904534/100‑1 (MQ=255)
                                                      aaCCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGTGCTAGCCTGATGGATCGCCATGCGGATTTCATGAATCGCTTTGATGCGATCATCGTCCGGc                                            >  1:190111/1‑100 (MQ=255)
                                                       aCCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCa                                           >  1:2137028/1‑100 (MQ=255)
                                                       aCCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCa                                           <  1:1016599/100‑1 (MQ=255)
                                                        cccACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCaa                                          >  2:1477610/1‑100 (MQ=255)
                                                         ccACAGCACGCAATCCACTGGCTCTTCACGAAAAGGGCTAGCCTGATGGATCGCCATGCGGATTGCATTAATCGCTTTGATGCGAGCATTTTCCGGCAAg                                         >  2:1603525/1‑100 (MQ=255)
                                                         ccACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAg                                         <  1:3271500/100‑1 (MQ=255)
                                                         ccACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAg                                         >  2:2615807/1‑100 (MQ=255)
                                                         ccACAGAACGAAATACACTGGCGCTTCACGAAAGGGGGTAGCCTGATGGATCGCCATGCGGATTTCATTAATTGCTTTGATGCGATCATCTTCCGGCAAg                                         >  2:965440/1‑100 (MQ=255)
                                                              gAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGc                                    >  2:496796/1‑100 (MQ=255)
                                                              gAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGc                                    >  1:2216339/1‑100 (MQ=255)
                                                               aaCGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCg                                   <  1:2840397/100‑1 (MQ=255)
                                                               aaCGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCg                                   <  1:341283/100‑1 (MQ=255)
                                                                aCGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCga                                  >  2:2378877/1‑100 (MQ=255)
                                                                aCGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCCTTGATGCGATCATCTTCCGGCAAGCTGGCga                                  >  1:1216477/1‑100 (MQ=255)
                                                                 cGCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCgag                                 <  1:1021972/100‑1 (MQ=255)
                                                                  gCAATCCACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCgaga                                <  2:3506605/100‑1 (MQ=255)
                                                                     aTCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAAt                             >  2:18218/1‑100 (MQ=255)
                                                                       ccACTGGCTCTTCACGAAAGGGGCTGGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCgaggagtag                           >  2:620327/1‑94 (MQ=255)
                                                                       ccACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTg                           >  1:701921/1‑100 (MQ=255)
                                                                       ccACTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTg                           >  2:3204604/1‑100 (MQ=255)
                                                                       ccACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTg                           <  2:385592/100‑1 (MQ=255)
                                                                          cttgCTCTTCACTAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAAtt                            <  1:78663/93‑1 (MQ=255)
                                                                          cTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGt                        <  1:1886285/100‑1 (MQ=255)
                                                                          cTGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGt                        <  1:1228188/100‑1 (MQ=255)
                                                                           tGGCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGtt                       <  2:920103/100‑1 (MQ=255)
                                                                            ggCTCTTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTa                      <  1:713748/100‑1 (MQ=255)
                                                                               tctTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAAt                   >  2:2092091/1‑100 (MQ=255)
                                                                               tctTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAAt                   <  1:432264/100‑1 (MQ=255)
                                                                               tctTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAAt                   <  1:1924610/100‑1 (MQ=255)
                                                                                ctTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATc                  >  1:485768/1‑100 (MQ=255)
                                                                                ctTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATc                  >  1:725020/1‑100 (MQ=255)
                                                                                ctTCACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATc                  >  1:93979/1‑100 (MQ=255)
                                                                                  tCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCt                >  1:987595/1‑100 (MQ=255)
                                                                                   cACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCGTCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTg               >  1:1235456/1‑100 (MQ=255)
                                                                                   cACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTg               >  2:2860398/1‑100 (MQ=255)
                                                                                   cACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTg               >  2:905087/1‑100 (MQ=255)
                                                                                   cACGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGAGCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATcc                 <  2:2564042/98‑1 (MQ=255)
                                                                                   cACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTg               <  1:1798235/100‑1 (MQ=255)
                                                                                    aCGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGc              >  1:3128600/1‑100 (MQ=255)
                                                                                    aCGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGc              >  1:3477930/1‑100 (MQ=255)
                                                                                     cGAAAGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCg             <  1:60483/100‑1 (MQ=255)
                                                                                     cGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCg             >  1:3619407/1‑100 (MQ=255)
                                                                                       aaaGGGGCTAGCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGtt           <  1:1181832/100‑1 (MQ=255)
                                                                                              cTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCgtt    <  2:1243911/100‑1 (MQ=255)
                                                                                              cTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCAGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCgtt    >  1:1363273/1‑100 (MQ=255)
                                                                                                aaCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCgttgt  <  1:1757851/100‑1 (MQ=255)
                                                                                                 |                                                                                                  
GGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGT  >  NC_000913/635164‑635359

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