Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,580,714 G→A S26L (TCG→TTG)  ydeM ← putative YdeN‑specific sulfatase‑maturating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,580,7140GA97.8% 288.4 / ‑2.1 89S26L (TCG→TTG) ydeMputative YdeN‑specific sulfatase‑maturating enzyme
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (48/39);  total (50/39)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

CCAGGTAAAATAGACCTGATTGCCAGACGCTGCGATATATTGTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGC  >  NC_000913/1580617‑1580811
                                                                                                 |                                                                                                 
ccAGGTAAAATAGACCTGATTGCCAGACGCTGCGATATATTGTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAAc                                                                                                 <  1:191655/100‑1 (MQ=255)
 cAGGTAAAATAGACCTGATTGCCAGACGCTGCGATATATTGTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAAct                                                                                                >  1:876791/1‑100 (MQ=255)
    gTAAAATAGACCTGATTGCCAGACGCTGCGATATATTGTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCtt                                                                                             >  2:150009/1‑100 (MQ=255)
    gTAAAATAGACCTGATTGCCAGACGCTGCGATATATTGTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCtt                                                                                             >  2:342382/1‑100 (MQ=255)
       aaaTAGACCTGATTGCCAGACGCTGCGATATATTGTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCttttt                                                                                          >  2:31136/1‑100 (MQ=255)
        aaTAGACCTGATTGCCAGACGCTGCGATATATTGTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTc                                                                                         <  1:1564411/100‑1 (MQ=255)
          tAGACCTGATTGCCAGACGCTGCGATATATTGTTTGATGAACTCTTCCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCaa                                                                                       <  2:1142123/100‑1 (MQ=255)
           aGACCTGATTGCCAGACGCTGCGATATATTGTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAg                                                                                      <  1:1097881/100‑1 (MQ=255)
            gACCTGATTGCCAGACGCTGCGATATATTGTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAgg                                                                                     <  2:974557/100‑1 (MQ=255)
             aCCTGATTGCCAGACGCTGCGATATATTGTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGt                                                                                    >  1:688290/1‑100 (MQ=255)
              cttccgatctaGACGCTGCGATATATTGTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTa                                                                                   <  1:1493035/90‑1 (MQ=255)
              ccTGATTGCCAGACGCTGCGATATATTGTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCTTTTTTTTATTGCAACGCTTTTTCAAGGTa                                                                                   >  1:1314037/1‑100 (MQ=255)
                tGATTGCCAGACGCTGCGATATATTGTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTaaa                                                                                 <  1:593306/100‑1 (MQ=255)
                tGATTGCCAGACGCTGCGATATATTGTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTaaa                                                                                 >  1:1273734/1‑100 (MQ=255)
                  aTTGCCAGACGCTGCGATATATTGTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTaaaa                                                                                <  2:1432261/99‑1 (MQ=255)
                    tGCCAGACGCTGCGATATATTGTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACa                                                                             <  1:471589/100‑1 (MQ=255)
                     gCCAGACGCTGCGATATATTGTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAg                                                                            >  1:58911/1‑100 (MQ=255)
                     gCCAGACGCTGCGATATATTGTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAg                                                                            >  1:725777/1‑100 (MQ=255)
                     gCCAGACGCTGCGATATATTGTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAg                                                                            >  1:272533/1‑100 (MQ=255)
                     gCCAGACGCTGCGATATATTGTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAg                                                                            >  2:1597702/1‑100 (MQ=255)
                       cAGACGCTGCGATATATTGTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTa                                                                          <  1:1555454/100‑1 (MQ=255)
                        aGACGCTGCGATATATTGTTTGATGAATTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAATCTGCAACTCTTTTTCAAGGTAAAAACAGTaa                                                                         <  1:414962/100‑1 (MQ=255)
                        aGACGCTGCGATATATTGTTTGATGAACTCTTTCAGATTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTa                                                                          <  1:278281/99‑1 (MQ=255)
                        aGACGCTGCGATATATTGTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTaagatcggaag                                                                         >  2:1493046/1‑91 (MQ=255)
                        aGACGCTGCGATATATTGTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTaa                                                                         >  1:1176733/1‑100 (MQ=255)
                         gACGCTGCGATATATTGTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGATAAAAACAGTAAt                                                                        >  1:348462/1‑100 (MQ=255)
                           cGCTGCGATATATTGTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATca                                                                      >  2:843962/1‑100 (MQ=255)
                           cGCTGCGATATATTGTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTATCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATca                                                                      <  2:1046768/100‑1 (MQ=255)
                             cTGCGATATATTGTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATcaca                                                                    <  1:342380/100‑1 (MQ=255)
                               gCGATATATTGTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACAtt                                                                  >  2:1417779/1‑100 (MQ=255)
                               gCGATATATTGTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATAACAtt                                                                  >  2:892984/1‑100 (MQ=255)
                                  atatatTGTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAAGTGCAACTCTTTTTCAAGGTAAAAAAAGTAATCACATTTga                                                               >  1:1343192/1‑100 (MQ=255)
                                  atatatTGTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTga                                                               >  1:1219865/1‑100 (MQ=255)
                                  atatatTGTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTg                                                                <  2:926385/99‑1 (MQ=255)
                                   tatatTGTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTgag                                                              <  1:1453251/100‑1 (MQ=255)
                                   tatatTGTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTgag                                                              <  1:214796/100‑1 (MQ=255)
                                    atatTGTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTgaga                                                             <  2:1533519/100‑1 (MQ=255)
                                    atatTGTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTgaga                                                             <  2:899354/100‑1 (MQ=255)
                                        tGTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTAc                                                         >  1:1143226/1‑100 (MQ=255)
                                         gTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACa                                                        >  2:23581/1‑100 (MQ=255)
                                         gTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACa                                                        >  2:23648/1‑100 (MQ=255)
                                           ctgatgaACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACAtt                                                      <  2:1262401/99‑1 (MQ=255)
                                            tgatgaACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTg                                                     >  2:917185/1‑100 (MQ=255)
                                            tgatgaACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTg                                                     >  1:1388493/1‑100 (MQ=255)
                                             gatgaACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACAttag                                                    >  1:1262390/1‑98 (MQ=255)
                                               tgaACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGaaa                                                  <  2:31120/100‑1 (MQ=255)
                                                gcactctgaaagagTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGaaaa                                                 >  2:1787502/11‑100 (MQ=255)
                                                 aaCTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAAc                                                >  2:1595414/1‑100 (MQ=255)
                                                   ctctTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTg                                              >  2:1676954/1‑100 (MQ=255)
                                                    tctTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTgg                                             <  2:1141731/100‑1 (MQ=255)
                                                     ctTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGa                                            >  1:158440/1‑100 (MQ=255)
                                                     ctTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGa                                            >  2:769501/1‑100 (MQ=255)
                                                      tttCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAg                                           <  2:700504/100‑1 (MQ=255)
                                                      tttCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAg                                           <  1:928428/100‑1 (MQ=255)
                                                       ttCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAgg                                          <  2:688294/100‑1 (MQ=255)
                                                       ttCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAgg                                          >  2:1236550/1‑100 (MQ=255)
                                                         cAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGc                                        >  1:1763943/1‑100 (MQ=255)
                                                           gagTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTTAGATTACATTGAAAACTGGAGGGCtt                                      >  2:1725679/1‑100 (MQ=255)
                                                            agTGCTGTCATCCATCCATTTTTTATGAGTAAACTTCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTg                                     <  2:77376/100‑1 (MQ=255)
                                                              tGCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGc                                   <  2:524270/100‑1 (MQ=255)
                                                               gCTGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGc                                   <  1:866147/99‑1 (MQ=255)
                                                                 tGTCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGt                                >  1:1807429/1‑100 (MQ=255)
                                                                  gTCATCCATCCATTTTTCATGAGTAAACTTCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTa                               <  2:1492369/100‑1 (MQ=255)
                                                                   tCATCCATCCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTaa                              <  2:160192/100‑1 (MQ=255)
                                                                    catccatcCATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAAc                             >  1:828839/1‑100 (MQ=255)
                                                                     atccatccATTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACg                            >  2:380433/1‑100 (MQ=255)
                                                                       ccatccatTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACg                            >  1:1105123/1‑98 (MQ=255)
                                                                        catccatTTTTCATGAGTAAACTGTAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGc                         <  2:114630/100‑1 (MQ=255)
                                                                        catccatTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGc                         >  1:43292/1‑100 (MQ=255)
                                                                        catccatTTTTCATGAGTAAACTACAACTCTTTTTCAAAGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGc                         >  1:43135/1‑100 (MQ=255)
                                                                          tccatTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCAt                       >  1:714743/1‑100 (MQ=255)
                                                                           ccatTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATc                      >  2:139655/1‑100 (MQ=255)
                                                                           ccatTTTTCATGAGTAAACTGCAACTCTTTTTCAAGCTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATc                      <  2:1273745/100‑1 (MQ=255)
                                                                            catTTTTCATGAGTGAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCg                     <  1:1236539/100‑1 (MQ=255)
                                                                            catTTTTCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCg                     >  2:1229109/1‑100 (MQ=255)
                                                                               ttttCATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTa                  <  1:398525/100‑1 (MQ=255)
                                                                                 ttCATGAGTAAACTGCAACTCTTTTTCAAGTTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTAtc                <  2:404209/100‑1 (MQ=255)
                                                                                   cATGAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTAtctc              >  1:752053/1‑100 (MQ=255)
                                                                                      gAGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCt           >  2:1163477/1‑100 (MQ=255)
                                                                                       aGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTc          >  2:1002367/1‑100 (MQ=255)
                                                                                       aGTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTc          >  1:368406/1‑100 (MQ=255)
                                                                                        gttatcTGCAACTCTTTTTCAAGGTAGAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTc          <  2:1247300/94‑1 (MQ=255)
                                                                                        gTAAACTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCa         <  2:1353304/100‑1 (MQ=255)
                                                                                          aaaCTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAAATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAAt       <  1:843956/100‑1 (MQ=255)
                                                                                           aaCTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATa      >  1:537936/1‑100 (MQ=255)
                                                                                           aaCTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATa      <  1:278626/100‑1 (MQ=255)
                                                                                            aCTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATaa     >  2:535442/1‑100 (MQ=255)
                                                                                             cTGCAACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAg    >  1:847968/1‑100 (MQ=255)
                                                                                               gCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAggc  >  2:797786/1‑100 (MQ=255)
                                                                                                 |                                                                                                 
CCAGGTAAAATAGACCTGATTGCCAGACGCTGCGATATATTGTTTGATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGC  >  NC_000913/1580617‑1580811

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: