Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 339,557 T→C 74.6% G411G (GGT→GGC yahG → DUF1116 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913339,5570TC74.6% 242.1 / 74.6 133G411G (GGT→GGCyahGDUF1116 family protein
Reads supporting (aligned to +/- strand):  ref base T (23/11);  new base C (42/57);  total (65/68)
Fisher's exact test for biased strand distribution p-value = 1.64e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 136 total aligned reads displayed.

GCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCG  >  NC_000913/339458‑339654
                                                                                                   |                                                                                                 
gcTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGt                                                                                                   >  1:2887115/1‑100 (MQ=255)
  tggtggGCGGCACGGTGGAAGAAGCTATTGATTTCTCACGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGctt                                                                                                 <  2:1668534/100‑1 (MQ=255)
    gtggGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTa                                                                                               >  1:2266217/1‑100 (MQ=255)
    gtggGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTa                                                                                               >  2:420644/1‑100 (MQ=255)
     tggGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTAt                                                                                              <  2:246968/100‑1 (MQ=255)
        gcggcACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATggg                                                                                           <  2:689622/100‑1 (MQ=255)
          ggcACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCg                                                                                         <  1:1109533/100‑1 (MQ=255)
          ggcACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCg                                                                                         <  1:189069/100‑1 (MQ=255)
           gcACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGt                                                                                        >  1:898527/1‑100 (MQ=255)
           gcACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGt                                                                                        >  2:2250782/1‑100 (MQ=255)
            cACGGTGGAAGAAGCTATTGATCTCTCCCGTCAGATGCGCGAAATCACCTTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTg                                                                                       <  1:1957189/100‑1 (MQ=255)
                  ggAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCg                                                                                 <  1:2377348/100‑1 (MQ=255)
                  ggAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCg                                                                                 <  1:1497225/100‑1 (MQ=255)
                  ggAAGAAGCTATTCATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGCGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCg                                                                                 >  2:1976833/1‑100 (MQ=255)
                   gaagaaGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCgg                                                                                <  1:2606083/100‑1 (MQ=255)
                    aagaagCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGc                                                                               <  1:1163814/100‑1 (MQ=255)
                     agaagCTATTGATTTCTCCCGTAAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCa                                                                              >  1:435100/1‑100 (MQ=255)
                      gaagCTATTGATTTCTCCCTTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCATCGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCaa                                                                             <  1:274697/100‑1 (MQ=255)
                         gCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCg                                                                          <  2:1081227/100‑1 (MQ=255)
                          cTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGa                                                                         <  1:2072794/100‑1 (MQ=255)
                          cTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGa                                                                         <  1:2145674/100‑1 (MQ=255)
                           tATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGAc                                                                        >  2:317245/1‑100 (MQ=255)
                           tATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGAc                                                                        >  1:1755732/1‑100 (MQ=255)
                            aTTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACa                                                                       <  1:2288829/100‑1 (MQ=255)
                               gATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGACAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATc                                                                     >  1:1786898/1‑99 (MQ=255)
                               gATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCa                                                                    >  1:1049342/1‑100 (MQ=255)
                                 tttCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCAcc                                                                  >  2:2440729/1‑100 (MQ=255)
                                 tttCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCAcc                                                                  >  1:2652164/1‑100 (MQ=255)
                                  ttCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCAccc                                                                 >  1:2683709/1‑100 (MQ=255)
                                  ttCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCAccc                                                                 <  2:426556/100‑1 (MQ=255)
                                    ctcCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCcgc                                                               >  2:130258/1‑100 (MQ=255)
                                    ctcCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCcgc                                                               >  2:1845681/1‑100 (MQ=255)
                                     tcCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCcgcg                                                              >  2:641512/1‑100 (MQ=255)
                                      cccGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCcgcg                                                              >  2:473832/1‑99 (MQ=255)
                                      cccGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGt                                                             >  2:2502508/1‑100 (MQ=255)
                                      cccGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGt                                                             >  2:377821/1‑100 (MQ=255)
                                         gTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTggg                                                          <  1:791207/100‑1 (MQ=255)
                                         gTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTggg                                                          <  1:2851281/100‑1 (MQ=255)
                                         gTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTggg                                                          <  1:641511/100‑1 (MQ=255)
                                         gTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTggg                                                          <  2:477561/100‑1 (MQ=255)
                                         gTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTggg                                                          <  2:2542970/100‑1 (MQ=255)
                                           cAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTa                                                        <  1:1574713/100‑1 (MQ=255)
                                           cAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTa                                                        <  1:1887417/100‑1 (MQ=255)
                                             gATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCGCCATTCCGCTGCTCGGCGTTATGGGCGGGCCGGCGGCAAGCGACATCACCCGCGTGGGTagc                                                      >  1:1299812/1‑100 (MQ=255)
                                             gATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTagc                                                      <  1:1316079/100‑1 (MQ=255)
                                              aTGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTagca                                                     <  2:1251513/100‑1 (MQ=255)
                                                 cgcgAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCg                                                  >  2:2868919/1‑100 (MQ=255)
                                                  gcgAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCgg                                                 >  2:2809220/1‑100 (MQ=255)
                                                       aTCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGATAGCAGCGGCATTc                                            >  1:1979782/1‑100 (MQ=255)
                                                          aCCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGc                                         <  1:2239913/100‑1 (MQ=255)
                                                           cccTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGcc                                        >  1:621177/1‑100 (MQ=255)
                                                           cccTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGcc                                        <  1:2242720/100‑1 (MQ=255)
                                                            ccTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCg                                       >  1:1099009/1‑100 (MQ=255)
                                                            ccTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCg                                       >  1:1531567/1‑100 (MQ=255)
                                                             cTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCg                                       >  1:2882151/1‑99 (MQ=255)
                                                               cGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCTGCGTGGGTAGCAGCGGCATTCTGCCgg                                      <  1:636046/98‑1 (MQ=255)
                                                               cGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTg                                    >  1:2163349/1‑100 (MQ=255)
                                                               cGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTg                                    <  2:1777346/100‑1 (MQ=255)
                                                               cGGTGAAAACCCCAACGTCACCATTCCGCTGCCCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGCCATTCTGCCgg                                      >  2:131724/1‑98 (MQ=255)
                                                                 gTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGAt                                  >  2:407798/1‑100 (MQ=255)
                                                                 gTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGAt                                  >  2:1458728/1‑100 (MQ=255)
                                                                  tGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATc                                 <  2:955231/100‑1 (MQ=255)
                                                                  tGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATc                                 <  2:950604/100‑1 (MQ=255)
                                                                  tGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATc                                 >  2:2249101/1‑100 (MQ=255)
                                                                  tGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATc                                 <  2:1053008/100‑1 (MQ=255)
                                                                   gAAAGCCCCGCCGTCCCCCGTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCa                                <  1:729311/100‑1 (MQ=255)
                                                                   gAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCa                                <  1:1627889/100‑1 (MQ=255)
                                                                     aaaCCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAac                              <  2:1001905/100‑1 (MQ=255)
                                                                     aaaCCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAac                              >  1:130275/1‑100 (MQ=255)
                                                                     aaaCCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAac                              >  1:1862162/1‑100 (MQ=255)
                                                                       aCCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAacac                            <  1:344583/100‑1 (MQ=255)
                                                                       aCCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAacac                            <  2:2853311/100‑1 (MQ=255)
                                                                        ccccAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACAcc                           >  2:1986439/1‑100 (MQ=255)
                                                                        ccccAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACAcc                           >  2:2177873/1‑100 (MQ=255)
                                                                        ccccAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACAcc                           >  2:2868776/1‑100 (MQ=255)
                                                                             aCGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCAt                      >  1:2664661/1‑100 (MQ=255)
                                                                              cGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATc                     >  1:405074/1‑100 (MQ=255)
                                                                                 cACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGcc                  <  2:1328556/100‑1 (MQ=255)
                                                                                 cACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGcc                  <  2:2466028/100‑1 (MQ=255)
                                                                                   ccATTCCTCTGCCCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCa                <  1:2473441/100‑1 (MQ=255)
                                                                                   ccATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCa                >  2:2011169/1‑100 (MQ=255)
                                                                                   ccATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCa                >  2:1572880/1‑100 (MQ=255)
                                                                                   ccATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCa                <  2:585373/100‑1 (MQ=255)
                                                                                   ccATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCa                >  2:478652/1‑100 (MQ=255)
                                                                                    cATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATGCTGCCGGTGATCAACACCGCCATCGCCCAt               <  1:2280870/100‑1 (MQ=255)
                                                                                     aTTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATa              <  1:2868908/100‑1 (MQ=255)
                                                                                     aTTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATa              >  1:740189/1‑100 (MQ=255)
                                                                                     aTTCCGCTGCTCGGCTTAAAAGGCGTGCCGTCGGCAATCGACATCACGCGCGTGGGTAGCAGCGGAAATATGCCGGTGATCAACACCGCCATTGCCCAt               >  1:1266980/1‑99 (MQ=255)
                                                                                      ttCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATaa             >  2:1206690/1‑100 (MQ=255)
                                                                                       tCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATaaa            <  1:2806825/100‑1 (MQ=255)
                                                                                       tCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATaaa            <  2:1186412/100‑1 (MQ=255)
                                                                                         cGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGa          >  1:1388171/1‑100 (MQ=255)
                                                                                          gctgctCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGAt         >  1:1496938/1‑100 (MQ=255)
                                                                                           ctgctCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGCGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATg        >  2:2225262/1‑100 (MQ=255)
                                                                                             gctCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAgccgag      >  2:1711589/1‑95 (MQ=255)
                                                                                               tCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCggg    <  2:1901685/100‑1 (MQ=255)
                                                                                               tCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCggg    >  1:1230651/1‑100 (MQ=255)
                                                                                                cGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGc   >  2:1867086/1‑100 (MQ=255)
                                                                                                 ggTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCg  >  1:855332/1‑100 (MQ=255)
                                                                                                 ggTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCg  <  2:2321564/100‑1 (MQ=255)
                                                                                                   |                                                                                                 
GCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCG  >  NC_000913/339458‑339654

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: