Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,517,951 (C)6→7 coding (285/522 nt) fecI ← RNA polymerase sigma‑19 factor, fec operon‑specific; ECF sigma factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,517,9451.C88.5% 248.4 / 25.2 78coding (291/522 nt)fecIRNA polymerase sigma‑19 factor, fec operon‑specific; ECF sigma factor
Reads supporting (aligned to +/- strand):  ref base . (6/3);  new base C (39/30);  total (45/33)
Fisher's exact test for biased strand distribution p-value = 7.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGAAACGCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑T‑CCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAG  >  NC_000913/4517846‑4518041
                                                                                                      |                                                                                                 
cAGAAACGCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑t                                                                                                    >  1:721630/1‑100 (MQ=255)
cAGAAACGCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGGTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑t                                                                                                    >  1:2551971/1‑100 (MQ=255)
 aGAAACGCGTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑T‑c                                                                                                  >  2:538641/1‑100 (MQ=255)
  gAAACGCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑T‑cc                                                                                                 >  2:2491913/1‑100 (MQ=255)
  gAAACGCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑T‑cc                                                                                                 <  1:2535996/100‑1 (MQ=255)
     aCGCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTGGCGTTCCTCAGGTGAAGGCG‑C‑T‑ccc‑cc                                                                                             <  2:481981/100‑1 (MQ=255)
       gCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑T‑ccc‑cccc                                                                                           >  1:1689121/1‑99 (MQ=255)
        cTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑T‑ccc‑cccct                                                                                          >  1:2410966/1‑98 (MQ=255)
          tCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑T‑ccc‑cccctcc                                                                                        <  1:297534/100‑5 (MQ=255)
           cACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑T‑ccc‑cccctccg                                                                                       >  2:1268457/1‑95 (MQ=255)
            aCGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCgg                                                                                       <  2:879957/100‑1 (MQ=255)
            aCGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCgg                                                                                       <  2:2793117/100‑1 (MQ=255)
            aCGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCgg                                                                                       >  2:1859903/1‑100 (MQ=255)
            aCGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCgg                                                                                       <  1:371838/100‑1 (MQ=255)
             cGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGc                                                                                      >  1:1711469/1‑100 (MQ=255)
             cGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGc                                                                                      <  2:2306059/100‑1 (MQ=255)
              gtgtTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCa                                                                                     <  1:225249/100‑1 (MQ=255)
               tgtTTTGCTGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCAt                                                                                    >  1:1546055/1‑100 (MQ=255)
                gtTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATa                                                                                   >  2:441728/1‑100 (MQ=255)
                gtTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATa                                                                                   >  1:2435660/1‑100 (MQ=255)
                     gCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAgcgc                                                                              >  2:250990/1‑100 (MQ=255)
                      ccGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGcc                                                                             <  1:924452/100‑1 (MQ=255)
                        gTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAg                                                                           <  1:2109421/100‑1 (MQ=255)
                         tttAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCatc                                                                          <  1:1950287/100‑3 (MQ=255)
                         tttAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGc                                                                          <  2:1492439/100‑1 (MQ=255)
                            aGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAgatcggaagag                                                                       >  2:66602/1‑90 (MQ=255)
                            aGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCAtc                                                                       >  2:336475/1‑100 (MQ=255)
                             gCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCAtct                                                                      >  2:1314804/1‑100 (MQ=255)
                             gCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCAtct                                                                      >  1:952448/1‑100 (MQ=255)
                              gtatctggagaTGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCAtctc                                                                     <  2:288081/90‑1 (MQ=255)
                               ccGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTcc                                                                    <  2:1595251/100‑1 (MQ=255)
                                 gTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAg                                                                  <  2:827590/100‑1 (MQ=255)
                                 gTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAg                                                                  <  1:829195/100‑1 (MQ=255)
                                    cAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATa                                                               <  1:1266661/100‑1 (MQ=255)
                                     aGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATAc                                                              >  2:2171759/1‑100 (MQ=255)
                                      gCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACg                                                             <  1:446104/100‑1 (MQ=255)
                                           cTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCtttt                                                        <  1:2187340/100‑1 (MQ=255)
                                           cTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCtttt                                                        >  2:1050017/1‑100 (MQ=255)
                                             gTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑T‑CCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTcc                                                     >  1:2795253/1‑100 (MQ=255)
                                              tCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTcc                                                     <  2:721628/100‑1 (MQ=255)
                                               cGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCa                                                    >  2:865790/1‑100 (MQ=255)
                                               cGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCa                                                    >  2:2867099/1‑100 (MQ=255)
                                                 aggagTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCGCAGGTGAAGGCG‑C‑T‑CCCACCCTCCGGCATAAGCGCCAGCATCTCCAGATACCCTTTTTCCAgc                                                  >  2:1249781/1‑100 (MQ=255)
                                                  ggagTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAgcg                                                 >  2:2809841/1‑100 (MQ=255)
                                                    agTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGcgcg                                               <  2:79350/100‑1 (MQ=255)
                                                     gTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGt                                              >  2:598022/1‑100 (MQ=255)
                                                     gTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGt                                              >  2:2321937/1‑100 (MQ=255)
                                                       tGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGttt                                            >  1:44251/1‑100 (MQ=255)
                                                         tAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCg                                          <  2:1842172/100‑1 (MQ=255)
                                                             gTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGa                                      >  1:268404/1‑100 (MQ=255)
                                                             gTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGa                                      >  2:2211680/1‑100 (MQ=255)
                                                             gTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGa                                      >  2:2243123/1‑100 (MQ=255)
                                                             gTCTCGAGTTGGCGTTCGCGTTCCTCAGGTGAAGGCG‑CAC‑CCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGa                                      >  2:1785238/1‑100 (MQ=255)
                                                               ctcGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGaaa                                    >  2:739432/1‑100 (MQ=255)
                                                               ctcGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGaaa                                    <  2:2003544/100‑1 (MQ=255)
                                                                tcGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAAc                                   >  1:2312726/1‑100 (MQ=255)
                                                                 cGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACa                                  <  1:2243136/100‑1 (MQ=255)
                                                                    gTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGt                               >  1:506854/1‑100 (MQ=255)
                                                                     ttGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTc                              <  1:2577541/100‑1 (MQ=255)
                                                                      tGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGCCGAAAcgcgctg                             <  2:1689108/100‑7 (MQ=255)
                                                                       ggCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTcca                            >  1:612338/1‑100 (MQ=255)
                                                                       ggCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTcca                            >  1:2489227/1‑100 (MQ=255)
                                                                         cTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTccacc                          >  1:1339481/1‑100 (MQ=255)
                                                                         cTTTCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTccacc                          <  1:250991/100‑1 (MQ=255)
                                                                         cTTTCGCGTTCCTCAGGTGAAGGCG‑C‑T‑CCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTccacca                         >  1:1290696/1‑100 (MQ=255)
                                                                           ttCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCAt                        <  1:1314809/100‑1 (MQ=255)
                                                                            tCGCGTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATc                       >  1:223794/1‑100 (MQ=255)
                                                                             cgcgTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAACATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCa                      >  1:2642419/1‑100 (MQ=255)
                                                                              gcgTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCAc                     >  2:1094126/1‑100 (MQ=255)
                                                                              gcgTTCCTCAGGTGAAGGCG‑C‑T‑CCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCAcg                    >  2:1985306/1‑100 (MQ=255)
                                                                                gTTCCTCAGGTGCCGGCGCC‑C‑CCC‑CCCTCCGGCATCAGCGACAGCATCTCCAGCTCCGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCa                      <  2:2410950/97‑1 (MQ=255)
                                                                                gTTCCTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCAcgc                   >  1:2185183/1‑100 (MQ=255)
                                                                                   ccTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGtt                >  1:1697069/1‑100 (MQ=255)
                                                                                    cTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGttt               >  2:2074474/1‑100 (MQ=255)
                                                                                    cTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGttt               <  2:506852/100‑1 (MQ=255)
                                                                                    cTCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGtt                >  1:485343/1‑99 (MQ=255)
                                                                                     tCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTg              <  2:2185170/100‑1 (MQ=255)
                                                                                     tCAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTg              >  1:547122/1‑100 (MQ=255)
                                                                                      cAGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTgg             <  1:441729/100‑1 (MQ=255)
                                                                                       aGGTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGc            <  2:1711456/100‑1 (MQ=255)
                                                                                        ggTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCg           >  1:1207263/1‑100 (MQ=255)
                                                                                        ggTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCg           >  2:293777/1‑100 (MQ=255)
                                                                                         gTGAAGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGa          >  1:818040/1‑100 (MQ=255)
                                                                                            ttccgatc‑T‑CCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAg       <  1:479159/93‑1 (MQ=255)
                                                                                            aaGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAg       >  2:567130/1‑100 (MQ=255)
                                                                                            aaGGCG‑C‑TCCCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAg       <  2:1808640/100‑1 (MQ=255)
                                                                                              ggCG‑C‑T‑ACC‑CCCTCCGGCATATGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGAt         <  2:2727344/95‑1 (MQ=255)
                                                                                                cg‑c‑T‑CCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCag  >  1:1081297/1‑100 (MQ=255)
                                                                                                      |                                                                                                 
CAGAAACGCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCG‑C‑T‑CCC‑CCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAG  >  NC_000913/4517846‑4518041

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