Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 857,578 C→G 15.5% intergenic (+23/+218) ybiT → / ← ybiU ABC‑F family putative regulatory ATPase/DUF1479 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913857,5780CG15.5% 508.2 / 77.2 182intergenic (+23/+218)ybiT/ybiUABC‑F family putative regulatory ATPase/DUF1479 family protein
Reads supporting (aligned to +/- strand):  ref base C (68/85);  new base G (17/11);  total (85/97)
Fisher's exact test for biased strand distribution p-value = 1.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.28e-01

Only 100 of 187 total aligned reads displayed.

TGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGAC  >  NC_000913/857479‑857676
                                                                                                   |                                                                                                  
tGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAc                                                                                                    >  1:21315/1‑100 (MQ=255)
tGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAc                                                                                                    <  1:1414051/100‑1 (MQ=255)
tGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAc                                                                                                    >  2:4310558/1‑100 (MQ=255)
 ggAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACt                                                                                                   >  1:2245892/1‑100 (MQ=255)
 ggAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACt                                                                                                   >  2:4272281/1‑100 (MQ=255)
       cACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCgg                                                                                              <  1:3228893/99‑1 (MQ=255)
        aCCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCgggg                                                                                            <  2:1357417/100‑1 (MQ=255)
         ccccGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGt                                                                                           <  1:2211817/100‑1 (MQ=255)
         ccccGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGt                                                                                           >  1:147721/1‑100 (MQ=255)
         ccccGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCGGGGt                                                                                           <  1:1898486/100‑1 (MQ=255)
         ccccGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCCGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCGGGGt                                                                                           <  1:957315/100‑1 (MQ=255)
          cccGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTg                                                                                          >  2:1257452/1‑100 (MQ=255)
          cccGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCGGGGTg                                                                                          <  2:1882630/100‑1 (MQ=255)
           ccGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAGAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCggg                                                                                             <  1:2031394/96‑1 (MQ=255)
            cGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCGGGGTGAc                                                                                        <  2:3519941/100‑1 (MQ=39)
                   cgTGATCGACTTTAGCGGTAATTACGAAGATTACTTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCGGGGTGAcggataa                                                                                  <  2:1708204/99‑5 (MQ=37)
                   cgTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCt                                                                                 <  1:1829130/100‑1 (MQ=255)
                    gTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCtg                                                                                >  2:835454/1‑100 (MQ=255)
                     tGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCtga                                                                               <  2:3954318/100‑1 (MQ=255)
                      gATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCtgat                                                                              >  1:705834/1‑100 (MQ=255)
                       aTCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCtgatg                                                                             <  2:1416121/100‑1 (MQ=255)
                       aTCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCtgatg                                                                             <  2:1314193/100‑1 (MQ=255)
                        tCGACTTTAGCGGTAATTACGATGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCtgatga                                                                            <  2:885187/100‑1 (MQ=255)
                        tCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCtgatga                                                                            >  1:3284058/1‑100 (MQ=255)
                           aCTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACaa                                                                         <  1:1684012/100‑1 (MQ=255)
                             tttAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAAt                                                                       <  1:3286980/100‑1 (MQ=255)
                                 gCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCaa                                                                   <  1:4385692/100‑1 (MQ=255)
                                 gCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCaa                                                                   >  1:4269678/1‑100 (MQ=255)
                                 gCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCaa                                                                   >  2:684718/1‑100 (MQ=255)
                                 gCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCaa                                                                   <  1:2725619/100‑1 (MQ=255)
                                    gTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCAGTCACCCCGATACTGTCGGGGTGAGGGCGGCTGATGACAAATGCAAAAt                                                                >  1:2020975/1‑100 (MQ=255)
                                        ttACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGcc                                                            >  2:3000542/1‑100 (MQ=255)
                                          aCGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTg                                                          >  1:706843/1‑100 (MQ=255)
                                          aCGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTg                                                          >  1:3139973/1‑100 (MQ=255)
                                           cGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGa                                                         <  2:3576355/100‑1 (MQ=255)
                                             aaGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATg                                                       <  1:1948898/100‑1 (MQ=255)
                                             aaGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATg                                                       >  2:3107576/1‑100 (MQ=255)
                                             aaGATTACCTGCGTAGTAAAGGGATCGAGTAAGAGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGAt                                                        >  1:368609/1‑99 (MQ=255)
                                             aaGATTACCTGCATAGTAAAGGGATCGAGTAAGTGATTATCCGTCATCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCGAAATAGCCTGATg                                                       <  2:4385013/100‑1 (MQ=255)
                                              aGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATgc                                                      <  1:1819187/100‑1 (MQ=255)
                                              aGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATgc                                                      >  1:995664/1‑100 (MQ=255)
                                               gATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATgcg                                                     <  1:2755602/100‑1 (MQ=255)
                                                 ttACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCt                                                   <  1:1652117/100‑1 (MQ=255)
                                                 ttACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGATGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCt                                                   >  1:1631596/1‑100 (MQ=255)
                                                   aCCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTAc                                                 <  2:918773/100‑1 (MQ=255)
                                                    ccTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACg                                                >  1:3959867/1‑100 (MQ=255)
                                                      tGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCt                                              >  1:66607/1‑100 (MQ=255)
                                                      tGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCt                                              >  2:546163/1‑100 (MQ=255)
                                                          tagtagAGGGATCGAGTATGTGATTATCCGTCACCCCGATACGGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCTCTACGCTTATc                                          <  1:1329440/100‑1 (MQ=38)
                                                          tagtaAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATc                                          <  2:2501393/100‑1 (MQ=255)
                                                          tagtaAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATc                                          <  2:1839365/100‑1 (MQ=255)
                                                            gtaAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAg                                        <  2:3800490/100‑1 (MQ=255)
                                                             taAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAgg                                       <  2:250876/100‑1 (MQ=255)
                                                              aaaGGGATGGAGTAAGTGATTATCCGGCACCCCGTTACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAGATAGCCTGATGCGCTACGCGTATCAGGc                                      <  1:3491805/100‑1 (MQ=38)
                                                              aaaGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGc                                      <  2:3880218/100‑1 (MQ=255)
                                                              aaaGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGc                                      >  2:2644687/1‑100 (MQ=255)
                                                              aaaGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATACAAAATAGCCTGATTCGCTACGCTTATCAGGc                                      >  1:1084056/1‑100 (MQ=255)
                                                                  ggATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTAc                                  >  2:1146458/1‑100 (MQ=255)
                                                                    aTCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGc                                <  1:2012318/100‑1 (MQ=255)
                                                                    aTCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGc                                <  2:4129057/100‑1 (MQ=255)
                                                                    aTCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGc                                <  2:3689204/100‑1 (MQ=255)
                                                                     tCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCa                               >  1:3994497/1‑100 (MQ=255)
                                                                       gAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGAGGCGCTACGCTTATCAGGCCTACACAGt                             >  1:734887/1‑100 (MQ=255)
                                                                          tAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTAcc                          <  2:1757928/100‑1 (MQ=255)
                                                                            aGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTg                        >  2:3340845/1‑100 (MQ=255)
                                                                             gTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGc                       <  2:3159862/100‑1 (MQ=255)
                                                                              tccgatctcGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCa                      <  1:2847726/92‑1 (MQ=255)
                                                                              tGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCa                      >  2:4057004/1‑100 (MQ=255)
                                                                               gATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCaa                     >  1:2581117/1‑100 (MQ=255)
                                                                               gATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCaa                     <  1:3901443/100‑1 (MQ=255)
                                                                                 gccgccGTCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAata                   >  2:197272/5‑100 (MQ=25)
                                                                                 gccgccGTCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAata                   >  2:3197999/5‑100 (MQ=25)
                                                                                   cgccGTCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAatata                 >  1:30456/3‑100 (MQ=25)
                                                                                    tCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCTTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAatatat                >  2:4368976/1‑100 (MQ=255)
                                                                                    tCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAatatat                >  1:1756467/1‑100 (MQ=255)
                                                                                    tCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAatatat                >  1:2806041/1‑100 (MQ=255)
                                                                                    gccGTCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTCCGCAGTACCTGCAatatat                >  2:1630347/2‑100 (MQ=25)
                                                                                    gccGTCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAatatat                >  2:2340548/2‑100 (MQ=25)
                                                                                    gccGTCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAatatat                >  2:636454/2‑100 (MQ=25)
                                                                                     ccGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATAtt               >  2:1682048/1‑100 (MQ=255)
                                                                                     ccGTCACCCCGATACTGACGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAatatat                <  1:4023052/99‑1 (MQ=255)
                                                                                     ccGTCACCCCGATAATGTAGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCCTATCAGGACTACGCAGTACCTGCAatatat                <  1:4124917/99‑1 (MQ=25)
                                                                                     ccGTCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGGCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATAtt               >  1:4350645/1‑100 (MQ=25)
                                                                                     ccGTCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATAtt               >  2:2095277/1‑100 (MQ=37)
                                                                                     ccGTCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATAtt               >  1:1548498/1‑100 (MQ=25)
                                                                                     ccGTCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATAtt               >  2:3793238/1‑100 (MQ=25)
                                                                                      cGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCaagatcgga              >  2:2847742/1‑93 (MQ=255)
                                                                                      cGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTg              <  1:1374664/100‑1 (MQ=255)
                                                                                       gTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGAAAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGa             >  1:2353257/1‑100 (MQ=255)
                                                                                        tCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGaa            <  2:289493/100‑1 (MQ=255)
                                                                                         cACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAAt           >  1:1050970/1‑100 (MQ=25)
                                                                                           ccccGTTACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCt         <  2:3686222/100‑1 (MQ=255)
                                                                                           ccccGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTGCCTGCAATATATTGATTc          >  1:4244474/1‑99 (MQ=255)
                                                                                           ccccGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCt         >  2:1933108/1‑100 (MQ=255)
                                                                                            cccGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTg        >  2:2198460/1‑100 (MQ=255)
                                                                                             ccGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGc       <  2:2029010/100‑1 (MQ=255)
                                                                                              cGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCa      <  2:2373261/100‑1 (MQ=255)
                                                                                               gATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGc       <  1:3761054/98‑1 (MQ=255)
                                                                                                aTACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAgg    <  2:3757932/100‑1 (MQ=255)
                                                                                                aTACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAgg    >  1:3947287/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
TGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGAC  >  NC_000913/857479‑857676

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: