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.6% 935.5 / ‑4.2 259coding (338/387 nt)yidWregulator protein for dgo operon
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (143/115);  total (143/116)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 287 total aligned reads displayed.

TCGTTAAACGCTTCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTCATTTTCCAGCACC  >  NC_000913/3872358‑3872536
                                                                                              |                                                                                      
tCGTTAAACGCTTCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCt                                                                                         <  2:2692331/93‑1 (MQ=255)
 cGTTAAACGCTTCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCtg                                                                                        >  1:3018303/1‑92 (MQ=255)
 cGTTAAACGCTTCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCtg                                                                                        >  2:1687170/1‑92 (MQ=255)
  gTTAAACGCTTCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCtgc                                                                                       >  2:2891607/1‑91 (MQ=255)
   ttAAACGCTTCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTgcc                                                                                     >  2:5773987/1‑93 (MQ=255)
    tAAACGCTTCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTgccc                                                                                    <  2:3701462/93‑2 (MQ=255)
     aaaCGCTTCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTgccca                                                                                   >  2:3184763/1‑91 (MQ=255)
      aaCGCTTCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCaa                                                                                   >  1:868405/1‑93 (MQ=255)
       aCGCTTCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCGGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAAc                                                                                  >  2:1200302/1‑93 (MQ=255)
       aCGCTTCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAAc                                                                                  <  1:542225/93‑1 (MQ=255)
        cGCTTCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACg                                                                                 >  2:3875899/1‑93 (MQ=255)
        cGCTTCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACg                                                                                 <  1:1340900/93‑1 (MQ=255)
         gCTTCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACgg                                                                                <  2:6069150/93‑1 (MQ=255)
          cTTCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACggg                                                                               >  2:1124135/1‑93 (MQ=255)
           ttCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGc                                                                              >  2:1596236/1‑93 (MQ=255)
             cgcgGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCaa                                                                            <  1:4898769/93‑1 (MQ=255)
              gcgGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAAt                                                                           <  1:4379610/93‑1 (MQ=255)
                ggTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCg                                                                         >  2:115599/1‑93 (MQ=255)
                ggTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCg                                                                         >  1:1988196/1‑93 (MQ=255)
                 gTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGc                                                                        >  2:2046297/1‑93 (MQ=255)
                    cTGATTGTTGGCAATCATCTCGTTCATCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCgg                                                                     <  1:127315/93‑1 (MQ=255)
                    cTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCgg                                                                     <  2:492975/93‑1 (MQ=255)
                    cTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCgg                                                                     <  1:4761563/93‑1 (MQ=255)
                    cTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCgg                                                                     >  1:2519082/1‑93 (MQ=255)
                    cTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCgg                                                                     >  2:3256285/1‑93 (MQ=255)
                     tGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGt                                                                    >  1:2847703/1‑93 (MQ=255)
                     tGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGt                                                                    >  2:2595539/1‑93 (MQ=255)
                     tGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGt                                                                    <  1:2645746/93‑1 (MQ=255)
                     tGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGt                                                                    >  1:3922075/1‑93 (MQ=255)
                      gATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGCCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGtt                                                                   <  2:142881/93‑1 (MQ=255)
                        ttgttgGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTcc                                                                 >  2:2374296/1‑93 (MQ=255)
                        ttgttgGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTcc                                                                 >  2:4555549/1‑93 (MQ=255)
                          gttgGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTccac                                                               >  1:2744846/1‑93 (MQ=255)
                          gttgGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTccac                                                               <  1:5899976/93‑1 (MQ=255)
                            tgGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTccacca                                                             >  2:5403623/1‑93 (MQ=255)
                              gCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGa                                                           >  1:4112355/1‑93 (MQ=255)
                              gCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGa                                                           <  2:3290402/93‑1 (MQ=255)
                               cAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGAt                                                          >  1:6121778/1‑93 (MQ=255)
                                aaTCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGAtt                                                         >  1:644255/1‑93 (MQ=255)
                                 atcatcTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGCttt                                                        >  2:901523/1‑93 (MQ=255)
                                 atcatcTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGAttt                                                        >  1:3507872/1‑93 (MQ=255)
                                   catcTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTTCCCAACGGGCAATCGCCGGTTCCACCAGATTTCg                                                      <  1:3971423/93‑1 (MQ=255)
                                     tctcGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCa                                                    <  2:3773650/93‑1 (MQ=255)
                                     tctcGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGAAGTCGCa                                                    >  2:5466561/1‑93 (MQ=255)
                                      ctcGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCAc                                                   <  2:5954932/93‑1 (MQ=255)
                                        cGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACtt                                                 >  2:4320921/1‑93 (MQ=255)
                                          ttCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCg                                               <  2:1986927/93‑1 (MQ=255)
                                          ttCAGCGCCGATTCAATCTGCGCAAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCCCCAGATTTCGCACTTCg                                               <  1:2285149/93‑1 (MQ=255)
                                            cAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCt                                             <  2:5191207/93‑1 (MQ=255)
                                             aGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTc                                            <  1:5974901/93‑1 (MQ=255)
                                              gcgcCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCa                                           <  1:2767576/93‑1 (MQ=255)
                                                gcCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATg                                         <  2:4429249/93‑1 (MQ=255)
                                                gcCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATg                                         >  2:3915174/1‑93 (MQ=255)
                                                gcCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATg                                         <  1:4496293/93‑1 (MQ=255)
                                                gcCGATTCAATCTGCGCCAGAGCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATg                                         <  2:1587600/93‑1 (MQ=255)
                                                 ccGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATgg                                        <  2:3226615/93‑1 (MQ=255)
                                                  cGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGc                                       <  1:21827/93‑1 (MQ=255)
                                                  cGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGc                                       >  1:4632315/1‑93 (MQ=255)
                                                   gATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCa                                      >  1:3536758/1‑93 (MQ=255)
                                                    aTTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCAc                                     >  2:688624/1‑93 (MQ=255)
                                                    aTTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCAc                                     >  2:180990/1‑93 (MQ=255)
                                                    aTTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCAc                                     >  1:3846583/1‑93 (MQ=255)
                                                     ttCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACt                                    <  1:5838119/93‑1 (MQ=255)
                                                         aTCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCTGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATa                                <  1:5665827/93‑1 (MQ=255)
                                                          tCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATaa                               >  2:4644339/1‑93 (MQ=255)
                                                          tCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATaa                               <  1:3296490/93‑1 (MQ=255)
                                                           cTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAg                              <  2:3492082/93‑1 (MQ=255)
                                                           cTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAg                              >  1:515857/1‑93 (MQ=255)
                                                            tGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGc                             <  2:4006134/93‑1 (MQ=255)
                                                            tGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGc                             <  2:5232722/93‑1 (MQ=255)
                                                             gcgcCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGcc                            <  1:5757006/93‑1 (MQ=255)
                                                               gcCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGt                          >  2:2817691/1‑93 (MQ=255)
                                                               gcCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGt                          >  1:3151290/1‑93 (MQ=255)
                                                               gcCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGt                          <  2:4822994/93‑1 (MQ=255)
                                                               gcCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGt                          >  2:4380273/1‑93 (MQ=255)
                                                                ccAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTg                         <  2:2483785/93‑1 (MQ=255)
                                                                 cAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAATGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTgg                        <  2:5029258/93‑1 (MQ=255)
                                                                 cAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTgg                        <  1:4296491/93‑1 (MQ=255)
                                                                 cAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTgg                        <  2:3650636/93‑1 (MQ=255)
                                                                  aGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTggg                       >  2:1279844/1‑93 (MQ=255)
                                                                  aGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTggg                       <  1:2725250/93‑1 (MQ=255)
                                                                    aTCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTc                     >  1:2040326/1‑93 (MQ=255)
                                                                    aTCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTc                     >  1:1030087/1‑93 (MQ=255)
                                                                     tCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCg                    <  1:1143126/93‑1 (MQ=255)
                                                                     tCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCg                    >  2:1074218/1‑93 (MQ=255)
                                                                      cGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGt                   <  2:1329445/93‑1 (MQ=255)
                                                                      cGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGt                   <  1:1481622/93‑1 (MQ=255)
                                                                        cTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAg                 >  2:778098/1‑93 (MQ=255)
                                                                          ggAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTc               >  2:1766116/1‑93 (MQ=255)
                                                                          ggAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTc               >  2:6051238/1‑93 (MQ=255)
                                                                           gAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTCa              <  1:1763117/93‑1 (MQ=255)
                                                                             aGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTCAtt            <  1:5806235/93‑1 (MQ=255)
                                                                                cgcgcgCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTCATTTTc         <  1:1679138/93‑1 (MQ=255)
                                                                                 gcgcgcTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTCATTTTcc        >  2:4449066/1‑93 (MQ=255)
                                                                                  cgcgcTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTCATTTTCCa       >  1:4206576/1‑93 (MQ=255)
                                                                                  cgcgcTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTCATTTTCCa       <  2:5691234/93‑1 (MQ=255)
                                                                                  cgcgcTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGGCGTAGTCATTTTCCa       >  2:4967807/1‑93 (MQ=255)
                                                                                   gcgcTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTCATTTTCCat      <  1:2362611/93‑2 (MQ=255)
                                                                                      ctctGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTCATTTTCCAGCAc   >  2:2113311/1‑93 (MQ=255)
                                                                                       tctGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTCATTTTCCAGCAcc  >  1:819684/1‑93 (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 < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: