Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,517,951:1 (C)6→7 100% 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.C74.0% 240.5 / 91.3 96coding (291/522 nt)fecIRNA polymerase sigma‑19 factor, fec operon‑specific; ECF sigma factor
Reads supporting (aligned to +/- strand):  ref base . (16/9);  new base C (34/37);  total (50/46)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 106 total aligned reads displayed.

CAGAAACGCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGC  >  NC_000913/4517846‑4518039
                                                                                                      |                                                                                                
cAGAAACGCTTCACGTGTTTGGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCt                                                                                                   >  1:1579272/1‑100 (MQ=255)
 aGAAACGCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTc                                                                                                 <  2:1403168/100‑1 (MQ=255)
 aGAAACGCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTc                                                                                                 <  2:1186830/100‑1 (MQ=255)
  gAAACGCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTcc                                                                                                >  1:534616/1‑100 (MQ=255)
   aaaCGCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTccc                                                                                               <  1:1326098/100‑1 (MQ=255)
      cGCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTcccccc                                                                                          >  1:1225623/1‑100 (MQ=255)
         ttCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTccccccctc                                                                                       <  2:716079/100‑4 (MQ=255)
         ttCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCTcc                                                                                       >  2:1848957/1‑100 (MQ=255)
          tCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTccccccctcc                                                                                      <  1:2427393/100‑5 (MQ=255)
            aCGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCgg                                                                                     >  2:2592147/1‑100 (MQ=255)
            aCGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCTCCGGc                                                                                    >  2:1379328/1‑100 (MQ=255)
            aCGTGTTTTGCCGGTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCAGGTTGCCTTGCGCGTTCCTCCGGTGAAGGCGCCcccgcccgccg                                                                                     >  2:2321471/1‑94 (MQ=255)
             cGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCTCCGGCa                                                                                   <  2:203905/100‑1 (MQ=255)
             cGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCTCCGGCa                                                                                   >  2:16223/1‑100 (MQ=255)
              gtgtTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCa                                                                                   <  2:1869696/100‑1 (MQ=255)
              gtgtTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGCGTCTCCAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCCCCCCCCCTCCgg                                                                                     <  1:226174/98‑1 (MQ=255)
               tgtTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCAt                                                                                  <  1:773382/100‑1 (MQ=255)
               tgtTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCAt                                                                                  >  2:2564866/1‑100 (MQ=255)
                gtTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATa                                                                                 <  2:709363/100‑1 (MQ=255)
                gtTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCTCCGGCATaa                                                                                >  1:696102/1‑100 (MQ=255)
                  tttGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAg                                                                               >  2:860778/1‑100 (MQ=255)
                  tttGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAg                                                                               <  2:1519754/100‑1 (MQ=255)
                    tGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAgcg                                                                             <  2:1102563/100‑1 (MQ=255)
                     gCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAgcgc                                                                            >  1:729479/1‑100 (MQ=255)
                       cGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCa                                                                          <  2:2144567/100‑1 (MQ=255)
                       cGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCTCCGGCATAAGCGCCAg                                                                         <  2:1010436/100‑1 (MQ=255)
                        gTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAg                                                                         >  1:2169648/1‑100 (MQ=255)
                         tttAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCTCCGGCATAAGCGCCAGCa                                                                       <  1:1877679/100‑1 (MQ=255)
                           tAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCAt                                                                      >  2:446042/1‑100 (MQ=255)
                            aGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCTCCGGCATAAGCGCCAGCAtct                                                                    >  1:2213634/1‑100 (MQ=255)
                             gCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCCTAAGCGCCAGCAtct                                                                    >  2:1145810/1‑100 (MQ=255)
                             gCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCAtct                                                                    >  2:2184056/1‑100 (MQ=255)
                              cccGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTcc                                                                  <  2:177662/100‑1 (MQ=255)
                               ccGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTcc                                                                  <  2:563059/100‑1 (MQ=255)
                                cGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCa                                                                 <  2:1565930/100‑1 (MQ=255)
                                cGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCa                                                                 <  2:1782775/100‑1 (MQ=255)
                                 gTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGa                                                               <  2:1225627/100‑1 (MQ=255)
                                  tCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGAt                                                              >  2:1375827/1‑100 (MQ=255)
                                   ccAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGAt                                                              <  2:2584261/100‑1 (MQ=255)
                                    cAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATa                                                             <  1:342982/100‑1 (MQ=255)
                                    cAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATa                                                             >  1:1828109/1‑100 (MQ=255)
                                    cAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATa                                                             >  1:1760203/1‑100 (MQ=255)
                                         tGCTGTCGAGGAGTTGTCGGGTCTCGAGTTGGCTTTCGCGTTACTCAGGTGAGGGCGCTCCCCCCCTTCGGCCTAAGCGCCAGCCTCCCCAgtt                                                              >  1:1987809/1‑92 (MQ=255)
                                            tGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCttttt                                                     <  2:1858063/100‑1 (MQ=255)
                                             gTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTc                                                    >  1:265112/1‑100 (MQ=255)
                                              tCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCTCCAGCATCTCCAGATACGCTTTTTcc                                                   <  2:1579281/100‑1 (MQ=255)
                                              tCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTcc                                                   <  1:742478/100‑1 (MQ=255)
                                              tCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTcc                                                   <  1:2002317/100‑1 (MQ=255)
                                               cGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCa                                                  >  1:1032420/1‑100 (MQ=255)
                                               cGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCa                                                  >  1:1987179/1‑99 (MQ=255)
                                               cGAGGAGTTGGAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTcc                                                   <  1:2093621/99‑1 (MQ=255)
                                                gaggagTTGTAGGGTCTCGAGTTGGCTTTCGCGTTGCTCAGGTGAAGGCGCTCCCCACCTCCGGCATACGCGCCCGCATCTCCCGATACGCTTTTTcc                                                   >  1:1725498/1‑98 (MQ=255)
                                                gaggagTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAgc                                                >  2:45283/1‑100 (MQ=255)
                                                       tgttaggTCTCGGGGTGGCTTTCGCCGTCCCCAAGTGAAGGCCCCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGtt                                           <  2:1429300/95‑1 (MQ=255)
                                                        gTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCg                                        >  1:468963/1‑100 (MQ=255)
                                                         tcgggTCTTGAGTGGTTTGTCGCTTTCCTCTGGTAAAGGCGCTCCCCCGCTCCGGCATAAGCGTCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCg                                        <  1:847239/98‑1 (MQ=255)
                                                         tAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCg                                        <  2:573573/100‑1 (MQ=255)
                                                          aGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCgg                                       <  2:2295169/100‑1 (MQ=255)
                                                          aGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCgg                                       <  1:1525835/100‑1 (MQ=255)
                                                          aGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAgcgagatcggaa                                      >  1:902733/1‑91 (MQ=255)
                                                            ggTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCTGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCg                                     <  1:1381525/100‑1 (MQ=255)
                                                              tctcGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGaa                                   >  2:1778430/1‑100 (MQ=255)
                                                               ctcGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGaaa                                  >  1:127407/1‑100 (MQ=255)
                                                               cgcGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGaa                                   <  2:44426/97‑1 (MQ=255)
                                                                tcGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAAc                                 >  1:2181414/1‑100 (MQ=255)
                                                                tcGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAAc                                 <  1:939553/100‑1 (MQ=255)
                                                                 cGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACa                                >  2:2310328/1‑100 (MQ=255)
                                                                 cGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACa                                <  2:1672489/100‑1 (MQ=255)
                                                                 cGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACa                                >  1:2167169/1‑100 (MQ=255)
                                                                   aGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAgg                              <  2:2062010/100‑1 (MQ=255)
                                                                   aGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAgg                              <  1:516828/100‑1 (MQ=255)
                                                                    gTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGt                             >  2:897220/1‑100 (MQ=255)
                                                                    gTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGt                             <  1:1207963/100‑1 (MQ=255)
                                                                     ttGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTc                            <  2:448490/100‑1 (MQ=255)
                                                                       ggCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTccac                         >  1:2590072/1‑100 (MQ=255)
                                                                        gCTTTCGCGTTCCTCAGGTGAAGTCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGt                             <  1:1270408/96‑1 (MQ=255)
                                                                         cTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCTGCTCGTTTCCGCGCAACAGGTccacc                        >  1:43997/1‑100 (MQ=255)
                                                                          tttCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTccacca                       <  2:1139001/100‑1 (MQ=255)
                                                                          tttCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTccacca                       <  1:2111424/100‑1 (MQ=255)
                                                                            tCGCGTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATc                     >  2:360568/1‑100 (MQ=255)
                                                                             cgcgTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCa                    >  2:1582614/1‑100 (MQ=255)
                                                                              gcgTTCCTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCAc                   >  1:1770019/1‑100 (MQ=255)
                                                                                   ccTCAGGTGACGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGGTTCGGCGAAACAGGTCCACCATCCcgcgat              >  1:1574897/1‑98 (MQ=255)
                                                                                   ccTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGtt              >  2:1890439/1‑100 (MQ=255)
                                                                                    cTCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGttt             >  2:1695877/1‑100 (MQ=255)
                                                                                    cTCAGGTGAAGGCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTg            >  1:445365/1‑100 (MQ=255)
                                                                                    cTCAGGTGAAGGCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTg            >  1:1850159/1‑100 (MQ=255)
                                                                                     tCAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTg            <  1:1453216/100‑1 (MQ=255)
                                                                                      cAGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACACTTTTTCCAGCGAGTTTGGGCGAACAAGGTCCACCATCACGCGttt             >  1:1776580/1‑98 (MQ=255)
                                                                                       aGGTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGc          >  2:2103722/1‑100 (MQ=255)
                                                                                        ggTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCg         >  2:1794362/1‑100 (MQ=255)
                                                                                        ggTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCg         >  2:1039138/1‑100 (MQ=255)
                                                                                         gTGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGa        <  1:202827/100‑1 (MQ=255)
                                                                                          tGAAGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGAt       <  1:1454492/100‑1 (MQ=255)
                                                                                          tGAAGGCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATa      >  1:720549/1‑100 (MQ=255)
                                                                                            aaGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAg     <  2:857773/100‑1 (MQ=255)
                                                                                             aGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGt    >  2:105025/1‑100 (MQ=255)
                                                                                             aGGCGCTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGt    >  2:108532/1‑100 (MQ=255)
                                                                                              ggCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGc  >  2:2578728/1‑100 (MQ=255)
                                                                                               gcgcTCCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGc  <  2:2005364/100‑1 (MQ=255)
                                                                                                      |                                                                                                
CAGAAACGCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGC  >  NC_000913/4517846‑4518039

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: