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

CAGAAACGCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAG  >  NC_000913/4517846‑4518041
                                                                                                      |                                                                                                 
cAGAAACGCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCt                                                                                                    >  1:721630/1‑100 (MQ=255)
cAGAAACGCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGGTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCt                                                                                                    >  1:2551971/1‑100 (MQ=255)
 aGAAACGCGTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTc                                                                                                  >  2:538641/1‑100 (MQ=255)
  gAAACGCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTcc                                                                                                 >  2:2491913/1‑100 (MQ=255)
  gAAACGCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTcc                                                                                                 <  1:2535996/100‑1 (MQ=255)
     aCGCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTGGCGTTCCTCAGGTGAAGGCGCTccccc                                                                                             <  2:481981/100‑1 (MQ=255)
       gCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTccccccc                                                                                           >  1:1689121/1‑99 (MQ=255)
        cTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTccccccct                                                                                          >  1:2410966/1‑98 (MQ=255)
          tCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTccccccctcc                                                                                        <  1:297534/100‑5 (MQ=255)
           cACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTccccccctccg                                                                                       >  2:1268457/1‑95 (MQ=255)
            aCGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCgg                                                                                       <  2:879957/100‑1 (MQ=255)
            aCGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCgg                                                                                       <  2:2793117/100‑1 (MQ=255)
            aCGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCgg                                                                                       >  2:1859903/1‑100 (MQ=255)
            aCGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCgg                                                                                       <  1:371838/100‑1 (MQ=255)
             cGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGc                                                                                      >  1:1711469/1‑100 (MQ=255)
             cGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGc                                                                                      <  2:2306059/100‑1 (MQ=255)
              gtgtTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCa                                                                                     <  1:225249/100‑1 (MQ=255)
               tgtTTTGCTGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCAt                                                                                    >  1:1546055/1‑100 (MQ=255)
                gtTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATa                                                                                   >  2:441728/1‑100 (MQ=255)
                gtTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATa                                                                                   >  1:2435660/1‑100 (MQ=255)
                     gCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAgcgc                                                                              >  2:250990/1‑100 (MQ=255)
                      ccGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGcc                                                                             <  1:924452/100‑1 (MQ=255)
                        gTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAg                                                                           <  1:2109421/100‑1 (MQ=255)
                         tttAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCatc                                                                          <  1:1950287/100‑3 (MQ=255)
                         tttAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGc                                                                          <  2:1492439/100‑1 (MQ=255)
                            aGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAgatcggaagag                                                                       >  2:66602/1‑90 (MQ=255)
                            aGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCAtc                                                                       >  2:336475/1‑100 (MQ=255)
                             gCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCAtct                                                                      >  2:1314804/1‑100 (MQ=255)
                             gCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCAtct                                                                      >  1:952448/1‑100 (MQ=255)
                              gtatctggagaTGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCAtctc                                                                     <  2:288081/90‑1 (MQ=255)
                               ccGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTcc                                                                    <  2:1595251/100‑1 (MQ=255)
                                 gTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAg                                                                  <  2:827590/100‑1 (MQ=255)
                                 gTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAg                                                                  <  1:829195/100‑1 (MQ=255)
                                    cAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATa                                                               <  1:1266661/100‑1 (MQ=255)
                                     aGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATAc                                                              >  2:2171759/1‑100 (MQ=255)
                                      gCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACg                                                             <  1:446104/100‑1 (MQ=255)
                                           cTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCtttt                                                        <  1:2187340/100‑1 (MQ=255)
                                           cTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCtttt                                                        >  2:1050017/1‑100 (MQ=255)
                                             gTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTcc                                                     >  1:2795253/1‑100 (MQ=255)
                                              tCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTcc                                                     <  2:721628/100‑1 (MQ=255)
                                               cGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCa                                                    >  2:865790/1‑100 (MQ=255)
                                               cGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCa                                                    >  2:2867099/1‑100 (MQ=255)
                                                 aggagTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCGCAGGTGAAGGCGCTCCCACCCTCCGGCATAAGCGCCAGCATCTCCAGATACCCTTTTTCCAgc                                                  >  2:1249781/1‑100 (MQ=255)
                                                  ggagTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAgcg                                                 >  2:2809841/1‑100 (MQ=255)
                                                    agTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGcgcg                                               <  2:79350/100‑1 (MQ=255)
                                                     gTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGt                                              >  2:598022/1‑100 (MQ=255)
                                                     gTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGt                                              >  2:2321937/1‑100 (MQ=255)
                                                       tGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGttt                                            >  1:44251/1‑100 (MQ=255)
                                                         tAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCg                                          <  2:1842172/100‑1 (MQ=255)
                                                             gTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGa                                      >  1:268404/1‑100 (MQ=255)
                                                             gTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGa                                      >  2:2211680/1‑100 (MQ=255)
                                                             gTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGa                                      >  2:2243123/1‑100 (MQ=255)
                                                             gTCTCGAGTTGGCGTTCGCGTTCCTCAGGTGAAGGCGCACCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGa                                      >  2:1785238/1‑100 (MQ=255)
                                                               ctcGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGaaa                                    >  2:739432/1‑100 (MQ=255)
                                                               ctcGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGaaa                                    <  2:2003544/100‑1 (MQ=255)
                                                                tcGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAAc                                   >  1:2312726/1‑100 (MQ=255)
                                                                 cGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACa                                  <  1:2243136/100‑1 (MQ=255)
                                                                    gTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGt                               >  1:506854/1‑100 (MQ=255)
                                                                     ttGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTc                              <  1:2577541/100‑1 (MQ=255)
                                                                      tGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGCCGAAAcgcgctg                             <  2:1689108/100‑7 (MQ=255)
                                                                       ggCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTcca                            >  1:612338/1‑100 (MQ=255)
                                                                       ggCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTcca                            >  1:2489227/1‑100 (MQ=255)
                                                                         cTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTccacc                          >  1:1339481/1‑100 (MQ=255)
                                                                         cTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTccacc                          <  1:250991/100‑1 (MQ=255)
                                                                         cTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTccacca                         >  1:1290696/1‑100 (MQ=255)
                                                                           ttCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCAt                        <  1:1314809/100‑1 (MQ=255)
                                                                            tCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATc                       >  1:223794/1‑100 (MQ=255)
                                                                             cgcgTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAACATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCa                      >  1:2642419/1‑100 (MQ=255)
                                                                              gcgTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCAc                     >  2:1094126/1‑100 (MQ=255)
                                                                              gcgTTCCTCAGGTGAAGGCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCAcg                    >  2:1985306/1‑100 (MQ=255)
                                                                                gTTCCTCAGGTGCCGGCGCCCCCCCCCTCCGGCATCAGCGACAGCATCTCCAGCTCCGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCa                      <  2:2410950/97‑1 (MQ=255)
                                                                                gTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCAcgc                   >  1:2185183/1‑100 (MQ=255)
                                                                                   ccTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGtt                >  1:1697069/1‑100 (MQ=255)
                                                                                    cTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGttt               >  2:2074474/1‑100 (MQ=255)
                                                                                    cTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGttt               <  2:506852/100‑1 (MQ=255)
                                                                                    cTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGtt                >  1:485343/1‑99 (MQ=255)
                                                                                     tCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTg              <  2:2185170/100‑1 (MQ=255)
                                                                                     tCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTg              >  1:547122/1‑100 (MQ=255)
                                                                                      cAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTgg             <  1:441729/100‑1 (MQ=255)
                                                                                       aGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGc            <  2:1711456/100‑1 (MQ=255)
                                                                                        ggTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCg           >  1:1207263/1‑100 (MQ=255)
                                                                                        ggTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCg           >  2:293777/1‑100 (MQ=255)
                                                                                         gTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGa          >  1:818040/1‑100 (MQ=255)
                                                                                            ttccgatcTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAg       <  1:479159/93‑1 (MQ=255)
                                                                                            aaGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAg       >  2:567130/1‑100 (MQ=255)
                                                                                            aaGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAg       <  2:1808640/100‑1 (MQ=255)
                                                                                              ggCGCTACCCCCTCCGGCATATGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGAt         <  2:2727344/95‑1 (MQ=255)
                                                                                                cgcTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCag  >  1:1081297/1‑100 (MQ=255)
                                                                                                      |                                                                                                 
CAGAAACGCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAG  >  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: