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% 807.7 / NA 253E113A (GAG→GCG) yidWregulator protein for dgo operon
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base G (118/134);  minor base C (0/1);  total (118/135)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 271 total aligned reads displayed.

TCGTTAAACGCTTCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTCATTTTCCAGCACC  >  NC_000913/3872358‑3872536
                                                                                             |                                                                                      
  gTTAAACGCTTCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCtgc                                                                                       <  2:2848843/93‑3 (MQ=255)
    tAAACGCTTCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTgccc                                                                                    >  2:2968507/1‑92 (MQ=255)
     aaaCGCTTCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTgccca                                                                                   >  2:2468339/1‑91 (MQ=255)
     aaaCGCTTCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTgccca                                                                                   >  2:4381739/1‑91 (MQ=255)
      aaCGCTTCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCaa                                                                                   >  1:3059575/1‑93 (MQ=255)
       aCGCTTCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAAc                                                                                  >  1:207568/1‑93 (MQ=255)
       aCGCTTCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGACGCGCGCTCTGCTGCCCAAc                                                                                  <  2:4511793/93‑1 (MQ=255)
          cTTCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCCACggg                                                                               >  1:946269/1‑93 (MQ=255)
            tCGCGGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCa                                                                             <  1:1350871/93‑1 (MQ=255)
             cgcgGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCaa                                                                            <  2:3759773/93‑1 (MQ=255)
             cgcgGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCaa                                                                            <  2:60137/93‑1 (MQ=255)
             cgcgGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCaa                                                                            <  2:3370611/93‑1 (MQ=255)
               cgGTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATc                                                                          <  2:1114638/93‑1 (MQ=255)
               cgGTCCTGATTGTTGGCAATCATCTCGTGCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATc                                                                          <  1:1245814/93‑1 (MQ=255)
                 gTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGc                                                                        <  1:708158/93‑1 (MQ=255)
                 gTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGc                                                                        <  2:3880041/93‑1 (MQ=255)
                   ccTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCg                                                                      >  2:2065690/1‑93 (MQ=255)
                    cTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCgg                                                                     >  2:3620532/1‑93 (MQ=255)
                     tGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGt                                                                    >  1:3648473/1‑93 (MQ=255)
                       aTTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTc                                                                  >  2:2025390/1‑93 (MQ=255)
                       aTTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTc                                                                  >  1:1405693/1‑93 (MQ=255)
                        ttgttgGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTcc                                                                 <  2:4169042/93‑1 (MQ=255)
                        ttgttgGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTcc                                                                 <  1:1532954/93‑1 (MQ=255)
                           ttgGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTccacc                                                              <  2:2943125/93‑1 (MQ=255)
                           ttgGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTccacc                                                              >  2:214615/1‑93 (MQ=255)
                            tgGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTccacca                                                             <  1:3662464/93‑1 (MQ=255)
                              gCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGa                                                           <  2:1290265/93‑1 (MQ=255)
                              gCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGa                                                           <  1:3956967/93‑1 (MQ=255)
                              gCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGa                                                           <  2:1976488/93‑1 (MQ=255)
                                aaTCATCTTGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGAtt                                                         >  1:3738943/1‑93 (MQ=255)
                                 atcatcTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGAttt                                                        >  2:4663171/1‑93 (MQ=255)
                                 atcatcTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGAttt                                                        <  1:27996/93‑1 (MQ=255)
                                  tcatcTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTc                                                       >  2:1627198/1‑93 (MQ=255)
                                  tcatcTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTc                                                       <  1:4573383/93‑1 (MQ=255)
                                     tctcGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCa                                                    <  2:4331036/93‑1 (MQ=255)
                                     tctcGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCa                                                    <  1:3510178/93‑1 (MQ=255)
                                     tctcGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCa                                                    >  2:4605770/1‑93 (MQ=255)
                                     tctcGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCa                                                    <  2:3509836/93‑1 (MQ=255)
                                     tctcGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCa                                                    <  1:2203423/93‑1 (MQ=255)
                                      ctcGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCAc                                                   <  1:4359129/93‑1 (MQ=255)
                                       tcGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACt                                                  >  2:2170129/1‑93 (MQ=255)
                                        cGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACtt                                                 >  2:4235049/1‑93 (MQ=255)
                                        cGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACtt                                                 <  2:1316876/93‑1 (MQ=255)
                                         gTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTc                                                <  2:4139714/93‑1 (MQ=255)
                                         gTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTc                                                <  1:3445447/93‑1 (MQ=255)
                                         gTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTCCCCAACGGTCAATCGCCGGTTCCACCAGATTTCGCACTTc                                                <  1:4304787/93‑1 (MQ=255)
                                          ttCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCg                                               >  2:3411518/1‑93 (MQ=255)
                                           tCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGc                                              >  1:4606893/1‑93 (MQ=255)
                                              gcgcCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCa                                           >  1:869655/1‑93 (MQ=255)
                                               cgcCGATTCAATCTGCGCCAGCTCGCTGGACGTCGCGCGCTCTTCTGCCTACCGGGCAATCGCCGGTTCCACCAAATTTCGCGCTTCGCTCAt                                          <  2:4366519/93‑1 (MQ=255)
                                               cgcCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCAt                                          >  1:2157558/1‑93 (MQ=255)
                                                gcCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATg                                         >  1:99476/1‑93 (MQ=255)
                                                gcCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATg                                         <  1:3015892/93‑1 (MQ=255)
                                                  cGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGc                                       >  1:1031621/1‑93 (MQ=255)
                                                   gATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCa                                      >  2:3100066/1‑93 (MQ=255)
                                                   gATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCa                                      >  2:3795626/1‑93 (MQ=255)
                                                    aTTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCAc                                     <  1:1971299/93‑1 (MQ=255)
                                                     ttCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACt                                    >  2:1712968/1‑93 (MQ=255)
                                                     ttCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACt                                    >  2:3966784/1‑93 (MQ=255)
                                                     ttCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACt                                    <  2:706776/93‑1 (MQ=255)
                                                       cAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGa                                  <  1:3819454/93‑1 (MQ=255)
                                                        aaTCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGAt                                 >  2:2287752/1‑93 (MQ=255)
                                                        aaTCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGAt                                 >  2:2551335/1‑93 (MQ=255)
                                                         aTCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATa                                >  1:1887557/1‑93 (MQ=255)
                                                          tCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATaa                               >  2:2686678/1‑93 (MQ=255)
                                                          tCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATaa                               <  1:1459943/93‑1 (MQ=255)
                                                            tGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGc                             <  2:2407402/93‑1 (MQ=255)
                                                            tGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGc                             <  2:3020526/93‑1 (MQ=255)
                                                             gcgcCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGcc                            <  1:3418028/93‑1 (MQ=255)
                                                             gcgcCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGcc                            <  1:2393941/93‑1 (MQ=255)
                                                               gcCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGt                          >  2:1454702/1‑93 (MQ=255)
                                                                ccAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTg                         <  1:2052385/93‑1 (MQ=255)
                                                                 cAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTgg                        <  2:1518089/93‑1 (MQ=255)
                                                                  aGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTggg                       >  1:3890051/1‑93 (MQ=255)
                                                                  aGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTggg                       >  1:219766/1‑93 (MQ=255)
                                                                   gATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGt                      >  1:3822656/1‑93 (MQ=255)
                                                                   gATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGt                      <  2:3137794/93‑1 (MQ=255)
                                                                    aTCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTc                     >  1:1655001/1‑93 (MQ=255)
                                                                    aTCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTc                     <  1:2563426/93‑1 (MQ=255)
                                                                     tCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCg                    >  2:18609/1‑93 (MQ=255)
                                                                      cGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGt                   >  2:4036292/1‑93 (MQ=255)
                                                                       gCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTa                  <  1:359502/93‑1 (MQ=255)
                                                                       gCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGAAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTa                  >  2:1183888/1‑93 (MQ=255)
                                                                        cTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAg                 >  1:768702/1‑93 (MQ=255)
                                                                         tGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGt                <  1:2922160/93‑1 (MQ=255)
                                                                         tGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGt                >  1:2792835/1‑93 (MQ=255)
                                                                          ggAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTc               >  2:537134/1‑93 (MQ=255)
                                                                             aGTCGCGCTCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTCAtt            <  2:2804040/93‑1 (MQ=255)
                                                                             aGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTCAtt            >  1:3264809/1‑93 (MQ=255)
                                                                             aGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTCAtt            <  2:625101/93‑1 (MQ=255)
                                                                             aGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTCAtt            >  2:1947643/1‑93 (MQ=255)
                                                                              ggcgcgcgCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTCAttt           <  2:330369/91‑1 (MQ=255)
                                                                              gTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTCAttt           <  1:1107055/93‑1 (MQ=255)
                                                                              gTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTCAttt           <  2:3261575/93‑1 (MQ=255)
                                                                                cgcgcgCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTCATTTTc         <  2:968012/93‑1 (MQ=255)
                                                                                 gcgcgcTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTCATTTTcc        >  2:782135/1‑93 (MQ=255)
                                                                                     gcTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTCATTTTCCAGCa    >  2:891804/1‑93 (MQ=255)
                                                                                      ctctGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTCATTTTCCAGCAc   >  2:1474241/1‑93 (MQ=255)
                                                                                       tctGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTCATTTTCCAGCAcc  <  1:4057200/93‑1 (MQ=255)
                                                                                       tctGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGTAGTAATTTTCCAGCAcc  >  2:742316/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 < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: