Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,872,450 +G coding (338/387 nt) yidW ← regulator protein for dgo operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,872,4501.G99.2% 1222.2 / 1.2 356coding (338/387 nt)yidWregulator protein for dgo operon
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base G (196/157);  total (197/159)
Fisher's exact test for biased strand distribution p-value = 5.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.32e-01

Only 100 of 394 total aligned reads displayed.

TCGTTAAACGCTTCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTCATTTTCCAGCACC  >  NC_000913/3872358‑3872536
                                                                                              |                                                                                       
 cGTTAAACGCTTCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCtg                                                                                         >  2:3793811/1‑92 (MQ=255)
  gTTAAACGCTTCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCtgc                                                                                        >  2:7982217/1‑91 (MQ=255)
  gTTAAACGCTTCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCtgc                                                                                        <  2:8372968/93‑3 (MQ=255)
    tAAACGCTTCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTgccc                                                                                     <  2:3391244/93‑2 (MQ=255)
     aaaCGCTTCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTgccca                                                                                    >  1:6161923/1‑91 (MQ=255)
      aaCGCTTCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAc                                                                                  >  2:2490733/1‑93 (MQ=255)
       aCGCTTCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTgcccgac                                                                                  <  1:6870690/93‑5 (MQ=255)
        cGCTTCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACg                                                                                 <  1:7060968/93‑1 (MQ=255)
         gCTTCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAAAgg                                                                                >  1:6499437/1‑93 (MQ=255)
           ttCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGc                                                                              >  2:2848888/1‑93 (MQ=255)
           ttCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGc                                                                              >  1:8286268/1‑93 (MQ=255)
            tCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCa                                                                             >  2:2903307/1‑93 (MQ=255)
            tCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGAGCTCTGCTGCCCAACGGGCa                                                                             >  2:8056614/1‑93 (MQ=255)
                 gTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGc                                                                        >  1:6374156/1‑93 (MQ=255)
                 gTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGc                                                                        >  2:8430691/1‑93 (MQ=255)
                 gTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGc                                                                        <  1:2499356/93‑1 (MQ=255)
                 gTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGc                                                                        <  2:6150995/93‑1 (MQ=255)
                  tCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGcc                                                                       >  2:8489720/1‑93 (MQ=255)
                     tGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGt                                                                    <  2:5048605/93‑1 (MQ=255)
                        ttgttgGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTcc                                                                 >  1:1087739/1‑93 (MQ=255)
                         tgttgGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTcca                                                                <  1:5990898/93‑1 (MQ=255)
                              gCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGa                                                           >  1:7580231/1‑93 (MQ=255)
                               cAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGAt                                                          <  1:4378862/93‑1 (MQ=255)
                               cAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGAt                                                          <  2:7119166/93‑1 (MQ=255)
                                aaTCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGAtt                                                         >  1:6569463/1‑93 (MQ=255)
                                 atcatcTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGAttt                                                        >  1:8639804/1‑93 (MQ=255)
                                 atcatcTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGAttt                                                        <  1:4546610/93‑1 (MQ=255)
                                 atcatcTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGAttt                                                        >  2:2016493/1‑93 (MQ=255)
                                 accatcTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGAttt                                                        >  1:8424834/3‑93 (MQ=255)
                                      ctcGTTCAGCGCCGATTCACTCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCACCGGGCAATCGCCGGTTCCACCCGATTTCGCAc                                                   <  1:7837754/93‑1 (MQ=255)
                                      ctcGTTCAGCGCCGATTAAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCAc                                                   <  2:8803472/93‑1 (MQ=255)
                                       tcGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACt                                                  >  1:3091302/1‑93 (MQ=255)
                                       tcGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACt                                                  <  1:8810372/93‑1 (MQ=255)
                                       tcGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACt                                                  >  2:1050293/1‑93 (MQ=255)
                                       tcGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACt                                                  <  1:7750922/93‑1 (MQ=255)
                                       tcGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACt                                                  >  2:1487362/1‑93 (MQ=255)
                                       tcGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACt                                                  >  1:3687740/1‑93 (MQ=255)
                                       tcGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACt                                                  >  2:4349448/1‑93 (MQ=255)
                                        cGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACtt                                                 <  1:6965913/93‑1 (MQ=255)
                                        cGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACtt                                                 <  1:5159303/93‑1 (MQ=255)
                                         gtacAGCGCCTAGTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTc                                                <  2:4602088/90‑1 (MQ=255)
                                         gTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTc                                                >  2:3637274/1‑93 (MQ=255)
                                         gTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTc                                                >  2:6174673/1‑93 (MQ=255)
                                          ttCAGCGCCGATTCAATCTGCGCCAGATCGCTGTAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGCTCCACCAGATTTCGCACTTCg                                               >  1:1788407/1‑93 (MQ=255)
                                          ttCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCg                                               >  1:5062773/1‑93 (MQ=255)
                                          ttCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCg                                               >  2:8116443/1‑93 (MQ=255)
                                           tCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGc                                              <  2:6691719/93‑1 (MQ=255)
                                           tCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGc                                              >  1:1296175/1‑93 (MQ=255)
                                             aGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTc                                            >  1:2677545/1‑93 (MQ=255)
                                               cgcCGATTCAATCTGCGCCAGATCGCTAGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCAt                                          <  1:3954110/93‑1 (MQ=255)
                                                gcCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATg                                         <  2:4678368/93‑1 (MQ=255)
                                                  cGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGc                                       <  1:752838/93‑1 (MQ=255)
                                                   gATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCa                                      >  1:7862486/1‑93 (MQ=255)
                                                   gATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCa                                      >  1:5268295/1‑93 (MQ=255)
                                                   gATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCa                                      >  2:3069818/1‑93 (MQ=255)
                                                    aTTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCAc                                     >  2:5682012/1‑93 (MQ=255)
                                                    aTTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCAc                                     >  1:8706758/1‑93 (MQ=255)
                                                    aTTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCAc                                     <  2:3240993/93‑1 (MQ=255)
                                                    aTTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCAc                                     <  2:1271115/93‑1 (MQ=255)
                                                     ttCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACt                                    <  2:2405801/93‑1 (MQ=255)
                                                     ttCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACt                                    >  1:4856640/1‑93 (MQ=255)
                                                     ttCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACt                                    >  2:3868805/1‑93 (MQ=255)
                                                      tCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTg                                   <  1:6779664/93‑1 (MQ=255)
                                                      tCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCACTGCAGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCACTCATGGCACTg                                   <  2:7834932/93‑1 (MQ=255)
                                                       cAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGa                                  <  2:6530669/93‑1 (MQ=255)
                                                       cAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGa                                  <  1:4810601/93‑1 (MQ=255)
                                                         agcTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATa                                <  2:7679884/91‑1 (MQ=255)
                                                            tGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGc                             <  2:6929155/93‑1 (MQ=255)
                                                            tGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGc                             <  2:1621574/93‑1 (MQ=255)
                                                            cgcgcCAGCTCGCTGGAAGTCGCGCACTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGc                             <  1:192717/92‑1 (MQ=255)
                                                             gcgcCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGcc                            <  1:53212/93‑1 (MQ=255)
                                                                ccAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTg                         <  2:2739179/93‑1 (MQ=255)
                                                                ccAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTg                         >  2:592464/1‑93 (MQ=255)
                                                                  aGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTggg                       <  2:6322399/93‑1 (MQ=255)
                                                                  aGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTggg                       <  2:2166437/93‑1 (MQ=255)
                                                                  aGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTggg                       >  2:4981641/1‑93 (MQ=255)
                                                                    aTCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTc                     >  2:3984281/1‑93 (MQ=255)
                                                                    aTCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTc                     >  1:7304883/1‑93 (MQ=255)
                                                                      cGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGt                   >  1:1180485/1‑93 (MQ=255)
                                                                       gCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTa                  >  1:2931941/1‑93 (MQ=255)
                                                                         tGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGt                <  2:1719201/93‑1 (MQ=255)
                                                                         tGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGt                <  2:8730377/93‑1 (MQ=255)
                                                                         tGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGt                >  1:3896733/1‑93 (MQ=255)
                                                                         tGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGt                >  1:4085380/1‑93 (MQ=255)
                                                                         tGGAAGTCGCGCGCTCTGCTGCACAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGt                <  2:5271049/93‑1 (MQ=255)
                                                                             attCGCACTCTCTGCTACGATACGGGCAATCGCCGGTTCCACCAGATTTCGCTCTTCGCTCATGGCACTGATAAGCCGTGGGTAGTAGTCAtt            <  2:4906222/91‑1 (MQ=255)
                                                                             aGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTCAtt            >  2:5763623/1‑93 (MQ=255)
                                                                             aGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTCAtt            <  1:4629903/93‑1 (MQ=255)
                                                                               tCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTCAtttt          <  2:1156392/93‑1 (MQ=255)
                                                                                 gcgcgcTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTCATTTTcc        >  1:1659352/1‑93 (MQ=255)
                                                                                   gcgcTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTCATTTTCCAg      >  2:5562537/1‑93 (MQ=255)
                                                                                   gcgcTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTCATTTTCCAg      >  2:7293208/1‑93 (MQ=255)
                                                                                    cgcTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTCATTTTCCAGc     >  1:5475503/1‑93 (MQ=255)
                                                                                    cgcTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTCATTTTCCAGc     <  2:19089/93‑1 (MQ=255)
                                                                                     gcTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTCATTTTCCAGCa    >  2:7533176/1‑93 (MQ=255)
                                                                                     gcTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTCATTTTCCAGCa    >  1:3315371/1‑93 (MQ=255)
                                                                                     cctctGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTCATTTTCCAGCa    <  1:8095393/92‑1 (MQ=255)
                                                                                      ctctGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTCATTTTCCAGCAc   >  1:3498073/1‑93 (MQ=255)
                                                                                       tctGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTCATTTTCCAGCAcc  >  1:4235184/1‑93 (MQ=255)
                                                                                       tctGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTCATTTTCCAGCAcc  <  1:740093/93‑1 (MQ=255)
                                                                                              |                                                                                       
TCGTTAAACGCTTCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTCATTTTCCAGCACC  >  NC_000913/3872358‑3872536

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