Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 524,866 A→G 47.4% T536A (ACC→GCC)  rhsD → Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913524,8660AG47.4% ‑4.9 / 235.8 176T536A (ACC→GCC) rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base A (48/45);  new base G (40/43);  total (88/88)
Fisher's exact test for biased strand distribution p-value = 7.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.27e-01

Only 100 of 177 total aligned reads displayed.

CATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGTCCA  >  NC_000913/524768‑524963
                                                                                                  |                                                                                                 
cATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGAc                                                                                                  >  2:1951650/1‑100 (MQ=255)
cATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGAc                                                                                                  >  1:3236751/1‑100 (MQ=255)
cATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGAc                                                                                                  >  2:2831652/1‑100 (MQ=255)
 aTCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGcc                                                                                                 <  1:733867/100‑1 (MQ=255)
   cgcgCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTg                                                                                               <  1:2103814/100‑1 (MQ=255)
      gcAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAg                                                                                            <  2:2422742/100‑1 (MQ=255)
       cAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGGGTTCCCGGCGAATACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTg                                                                                               <  2:1769855/96‑1 (MQ=255)
          cGGGCAGAAAGTGCGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTgg                                                                                              <  2:3896803/94‑1 (MQ=255)
            gggcggcAGTGCGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGc                                                                                        <  2:3116079/92‑1 (MQ=255)
             ggAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTAc                                                                                     >  2:518583/1‑100 (MQ=255)
              gagaCAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCGGCACCCGGCAGATGACCTGGGGCCGCTACg                                                                                    >  1:1495975/1‑100 (MQ=255)
              gagaCAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACg                                                                                    <  1:3729087/100‑1 (MQ=255)
               agaCAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACgg                                                                                   >  1:4115745/1‑100 (MQ=255)
                gaCAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGCGCCGCTACggg                                                                                  >  1:2936206/1‑100 (MQ=255)
                gaCAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACggg                                                                                  <  2:549465/100‑1 (MQ=255)
                 aCAGTACGCTACCGCTACGATGACGCGCACAGGGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGc                                                                                 <  1:1652570/100‑1 (MQ=255)
                  cAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGc                                                                                 >  2:452002/1‑99 (MQ=255)
                  cAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCa                                                                                >  1:2996830/1‑100 (MQ=255)
                   agtacccggccgctacGATGACGCGAACTGTCGGTTACCGGCGTCGACCACGGAGTCGACGGGCAGCACCCGGCAGATGACCTGGTGCCGCt                                                                                       <  2:1585910/83‑1 (MQ=255)
                    gTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGc                                                                                 <  2:2213537/97‑1 (MQ=255)
                    gTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGt                                                                              >  2:466757/1‑100 (MQ=255)
                     tACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGtt                                                                             <  2:1108773/100‑1 (MQ=255)
                       cgctaccgctacGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGc                                                                           >  2:302435/1‑100 (MQ=255)
                       cgctaccgctacGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGc                                                                           >  1:2643557/1‑100 (MQ=255)
                       cgctaccgctacGATGACGCGCACAATGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTCCGGGCAGTTGc                                                                           >  1:2023140/1‑100 (MQ=255)
                          taccgctacGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTgg                                                                        <  2:2629599/100‑1 (MQ=255)
                           accgctacGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTgg                                                                        <  2:89437/99‑1 (MQ=255)
                            ccgctacGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCg                                                                      <  1:2521590/100‑1 (MQ=255)
                               ctacGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTc                                                                   <  1:1604239/100‑1 (MQ=255)
                                  cGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCAcc                                                                >  1:2469544/1‑100 (MQ=255)
                                  cGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCAcc                                                                <  1:4076235/100‑1 (MQ=255)
                                   gATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCg                                                               >  2:3541094/1‑100 (MQ=255)
                                      gACGCGCACAGTGAGTTACCGGCGACGACAACGCATGCGACGGGCAGCACCCGGCACATGACCTGGAGGCGCTACGGGCAGGTCCTGGCGTTCCCCGAc                                                             >  2:2618697/1‑99 (MQ=255)
                                       aCGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTg                                                           <  1:345125/100‑1 (MQ=255)
                                        cgcgCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGc                                                          >  1:181721/1‑100 (MQ=255)
                                        cgcgCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGc                                                          >  1:2994910/1‑100 (MQ=255)
                                            cacaGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCgg                                                      <  2:1701208/100‑1 (MQ=255)
                                              caGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGc                                                    <  1:1133446/100‑1 (MQ=255)
                                              caGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGc                                                    <  2:1307947/100‑1 (MQ=255)
                                                gTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTa                                                  >  1:3916153/1‑100 (MQ=255)
                                                  gAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCTGGCTAcc                                                <  1:2171600/100‑1 (MQ=255)
                                                  gAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTAcc                                                >  2:3887144/1‑100 (MQ=255)
                                                  gAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTAcc                                                <  2:1323035/100‑1 (MQ=255)
                                                   aGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCa                                               <  1:3524512/100‑1 (MQ=255)
                                                    gTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAg                                              >  1:3250340/1‑100 (MQ=255)
                                                    gTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAg                                              >  1:560653/1‑100 (MQ=255)
                                                     ttACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGa                                             >  1:1209024/1‑100 (MQ=255)
                                                      tACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCCGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGAc                                            <  1:2920464/100‑1 (MQ=255)
                                                      tACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGAc                                            >  1:810426/1‑100 (MQ=255)
                                                      tACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGAc                                            <  2:781197/100‑1 (MQ=255)
                                                        ccGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGAccc                                          >  2:1478204/1‑100 (MQ=255)
                                                         cGGCGACGACAACGGATGCGACGGGCAGCACTCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCg                                         <  2:734526/100‑1 (MQ=255)
                                                         cGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCg                                         >  1:2542167/1‑100 (MQ=255)
                                                           gCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGtt                                       <  2:1557546/100‑1 (MQ=255)
                                                            cgacgaCAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTa                                      <  1:2333610/100‑1 (MQ=255)
                                                            cgacgaCAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTa                                      <  1:3184227/100‑1 (MQ=255)
                                                              acgacAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATg                                    <  2:468467/100‑1 (MQ=255)
                                                              acgacAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATg                                    <  2:2845709/100‑1 (MQ=255)
                                                              acgacAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATg                                    <  2:3434656/100‑1 (MQ=255)
                                                              acgacAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATg                                    <  2:3457268/100‑1 (MQ=255)
                                                               cgacAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGa                                   >  1:2488895/1‑100 (MQ=255)
                                                               cgacAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGa                                   >  1:1099715/1‑100 (MQ=255)
                                                                 acAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAAt                                 >  1:1281002/1‑100 (MQ=255)
                                                                 acAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAAt                                 <  1:2049166/100‑1 (MQ=255)
                                                                 acAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAAt                                 >  2:2526242/1‑100 (MQ=255)
                                                                  cAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCCGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATa                                >  2:2930038/1‑100 (MQ=255)
                                                                  cAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATa                                <  1:3308897/100‑1 (MQ=255)
                                                                  cAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATa                                <  1:1044128/100‑1 (MQ=255)
                                                                  cAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATa                                <  2:640333/100‑1 (MQ=255)
                                                                  cAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATa                                >  1:3610448/1‑100 (MQ=255)
                                                                      ggATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGAc                            <  1:1547786/100‑1 (MQ=255)
                                                                        aTGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCg                          <  2:3905197/100‑1 (MQ=255)
                                                                         tGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGc                         >  2:265057/1‑100 (MQ=255)
                                                                         tGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCg                          <  2:1270224/99‑1 (MQ=255)
                                                                          gCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCt                        <  1:2724921/100‑1 (MQ=255)
                                                                           cGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTGCCAGACCCGTTATGAATACGACCGCtt                       >  2:2351779/1‑100 (MQ=255)
                                                                           cGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCtt                       >  2:1976203/1‑100 (MQ=255)
                                                                            gACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTc                      >  1:698824/1‑100 (MQ=255)
                                                                            gACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTc                      <  2:1582026/100‑1 (MQ=255)
                                                                             aCGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCg                     >  1:1244574/1‑100 (MQ=255)
                                                                                ggCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGcc                  >  1:3275419/1‑100 (MQ=255)
                                                                                 gcagcaCCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCa                 >  1:4020020/1‑100 (MQ=255)
                                                                                  cagcaCCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAg                <  1:896357/100‑1 (MQ=255)
                                                                                   agcaCCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGa               <  1:3980563/100‑1 (MQ=255)
                                                                                   agcaCCCGGCAGATGACCTGGAGCCGCTACGGGCATTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGa               <  2:2293253/100‑1 (MQ=255)
                                                                                    gcaCCCGGCAGCTGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGAt              >  2:3077596/1‑100 (MQ=255)
                                                                                    gcaCCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGAt              <  2:2812041/100‑1 (MQ=255)
                                                                                    gcaCCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGa               <  2:181720/99‑1 (MQ=255)
                                                                                    gcaCCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGAt              >  2:2755959/1‑100 (MQ=255)
                                                                                      aCCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGa            >  2:3950982/1‑100 (MQ=255)
                                                                                       cccGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGGTTCGGCCAGCTGAc           >  2:4076005/1‑100 (MQ=255)
                                                                                       cccGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGAc           >  1:1948958/1‑100 (MQ=255)
                                                                                       cccGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGAc           >  1:52897/1‑100 (MQ=255)
                                                                                        ccGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGAcg          >  1:3045297/1‑100 (MQ=255)
                                                                                        ccGGCAGATCACCTGGAGCCGCCACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTAGGAATACGACCGCTACGGCCAGATGAcg          >  2:1317712/1‑100 (MQ=255)
                                                                                          ggCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGAcggc        <  2:1923059/100‑1 (MQ=255)
                                                                                            cAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGAcggcgg      <  1:172915/100‑1 (MQ=255)
                                                                                             aGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGt     <  2:3164493/100‑1 (MQ=255)
                                                                                              gctGACCTGGAGACGCTACGGGCAGTTGATGGCGTTCCCCGAATGCTCGTGCTACCAGACTCTTTATGACTACGACCGCTTCGACCAGATGACGGCGGt     <  2:1846198/97‑1 (MQ=37)
                                                                                               aTGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGTcc   <  2:623220/100‑1 (MQ=255)
                                                                                                  |                                                                                                 
CATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGTCCA  >  NC_000913/524768‑524963

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