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.

CAGAAACGCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CT‑CCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGC  >  NC_000913/4517846‑4518039
                                                                                                      |                                                                                                
cAGAAACGCTTCACGTGTTTGGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑Ct                                                                                                   >  1:1579272/1‑100 (MQ=255)
 aGAAACGCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CT‑c                                                                                                 <  2:1403168/100‑1 (MQ=255)
 aGAAACGCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CT‑c                                                                                                 <  2:1186830/100‑1 (MQ=255)
  gAAACGCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CT‑cc                                                                                                >  1:534616/1‑100 (MQ=255)
   aaaCGCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CT‑ccc                                                                                               <  1:1326098/100‑1 (MQ=255)
      cGCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CT‑cccc‑c‑c                                                                                          >  1:1225623/1‑100 (MQ=255)
         ttCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CT‑cccc‑c‑cctc                                                                                       <  2:716079/100‑4 (MQ=255)
         ttCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CT‑CCCC‑C‑CTcc                                                                                       >  2:1848957/1‑100 (MQ=255)
          tCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CT‑cccc‑c‑cctcc                                                                                      <  1:2427393/100‑5 (MQ=255)
            aCGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCgg                                                                                     >  2:2592147/1‑100 (MQ=255)
            aCGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CT‑CCCC‑C‑CTCCGGc                                                                                    >  2:1379328/1‑100 (MQ=255)
            aCGTGTTTTGCCGGTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCAGGTTGCCTTGCGCGTTCCTCCGGTGAAGG‑CG‑CC‑cccg‑c‑ccgccg                                                                                     >  2:2321471/1‑94 (MQ=255)
             cGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CT‑CCCC‑C‑CTCCGGCa                                                                                   <  2:203905/100‑1 (MQ=255)
             cGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CT‑CCCC‑C‑CTCCGGCa                                                                                   >  2:16223/1‑100 (MQ=255)
              gtgtTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCa                                                                                   <  2:1869696/100‑1 (MQ=255)
              gtgtTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGCGTCTCCAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CGCCC‑CCCC‑C‑CTCCgg                                                                                     <  1:226174/98‑1 (MQ=255)
               tgtTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCAt                                                                                  <  1:773382/100‑1 (MQ=255)
               tgtTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCAt                                                                                  >  2:2564866/1‑100 (MQ=255)
                gtTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATa                                                                                 <  2:709363/100‑1 (MQ=255)
                gtTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CT‑CCCC‑C‑CTCCGGCATaa                                                                                >  1:696102/1‑100 (MQ=255)
                  tttGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAg                                                                               >  2:860778/1‑100 (MQ=255)
                  tttGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAg                                                                               <  2:1519754/100‑1 (MQ=255)
                    tGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAgcg                                                                             <  2:1102563/100‑1 (MQ=255)
                     gCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAgcgc                                                                            >  1:729479/1‑100 (MQ=255)
                       cGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCa                                                                          <  2:2144567/100‑1 (MQ=255)
                       cGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CT‑CCCC‑C‑CTCCGGCATAAGCGCCAg                                                                         <  2:1010436/100‑1 (MQ=255)
                        gTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAg                                                                         >  1:2169648/1‑100 (MQ=255)
                         tttAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CT‑CCCC‑C‑CTCCGGCATAAGCGCCAGCa                                                                       <  1:1877679/100‑1 (MQ=255)
                           tAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCAt                                                                      >  2:446042/1‑100 (MQ=255)
                            aGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CT‑CCCC‑C‑CTCCGGCATAAGCGCCAGCAtct                                                                    >  1:2213634/1‑100 (MQ=255)
                             gCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCCTAAGCGCCAGCAtct                                                                    >  2:1145810/1‑100 (MQ=255)
                             gCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCAtct                                                                    >  2:2184056/1‑100 (MQ=255)
                              cccGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CT‑CCCC‑C‑CTCCGGCATAAGCGCCAGCATCTcc                                                                  <  2:177662/100‑1 (MQ=255)
                               ccGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTcc                                                                  <  2:563059/100‑1 (MQ=255)
                                cGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCa                                                                 <  2:1565930/100‑1 (MQ=255)
                                cGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCa                                                                 <  2:1782775/100‑1 (MQ=255)
                                 gTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CT‑CCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGa                                                               <  2:1225627/100‑1 (MQ=255)
                                  tCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CT‑CCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGAt                                                              >  2:1375827/1‑100 (MQ=255)
                                   ccAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGAt                                                              <  2:2584261/100‑1 (MQ=255)
                                    cAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATa                                                             <  1:342982/100‑1 (MQ=255)
                                    cAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATa                                                             >  1:1828109/1‑100 (MQ=255)
                                    cAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATa                                                             >  1:1760203/1‑100 (MQ=255)
                                         tGCTGTCGAGGAGTTGTCGGGTCTCGAGTTGGCTTTCGCGTTACTCAGGTGAGGG‑CG‑CTCCCCC‑C‑CTTCGGCCTAAGCGCCAGCCTCCCCAgtt                                                              >  1:1987809/1‑92 (MQ=255)
                                            tGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCttttt                                                     <  2:1858063/100‑1 (MQ=255)
                                             gTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTc                                                    >  1:265112/1‑100 (MQ=255)
                                              tCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCTCCAGCATCTCCAGATACGCTTTTTcc                                                   <  2:1579281/100‑1 (MQ=255)
                                              tCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTcc                                                   <  1:742478/100‑1 (MQ=255)
                                              tCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTcc                                                   <  1:2002317/100‑1 (MQ=255)
                                               cGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCa                                                  >  1:1032420/1‑100 (MQ=255)
                                               cGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CT‑CCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCa                                                  >  1:1987179/1‑99 (MQ=255)
                                               cGAGGAGTTGGAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTcc                                                   <  1:2093621/99‑1 (MQ=255)
                                                gaggagTTGTAGGGTCTCGAGTTGGCTTTCGCGTTGCTCAGGTGAAGG‑CG‑CT‑CCCCAC‑CTCCGGCATACGCGCCCGCATCTCCCGATACGCTTTTTcc                                                   >  1:1725498/1‑98 (MQ=255)
                                                gaggagTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CT‑CCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAgc                                                >  2:45283/1‑100 (MQ=255)
                                                       tgttaggTCTCGGGGTGGCTTTCGCCGTCCCCAAGTGAAGGCCC‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGtt                                           <  2:1429300/95‑1 (MQ=255)
                                                        gTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CT‑CCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCg                                        >  1:468963/1‑100 (MQ=255)
                                                         tcgggTCTTGAGTGGTTTGTCGCTTTCCTCTGGTAAAGG‑CG‑CT‑CCCC‑CGCTCCGGCATAAGCGTCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCg                                        <  1:847239/98‑1 (MQ=255)
                                                         tAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCg                                        <  2:573573/100‑1 (MQ=255)
                                                          aGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCgg                                       <  2:2295169/100‑1 (MQ=255)
                                                          aGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCgg                                       <  1:1525835/100‑1 (MQ=255)
                                                          aGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CT‑CCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAgcgagatcggaa                                      >  1:902733/1‑91 (MQ=255)
                                                            ggTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCTGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCg                                     <  1:1381525/100‑1 (MQ=255)
                                                              tctcGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGaa                                   >  2:1778430/1‑100 (MQ=255)
                                                               ctcGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGaaa                                  >  1:127407/1‑100 (MQ=255)
                                                               cgcGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGaa                                   <  2:44426/97‑1 (MQ=255)
                                                                tcGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAAc                                 >  1:2181414/1‑100 (MQ=255)
                                                                tcGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAAc                                 <  1:939553/100‑1 (MQ=255)
                                                                 cGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACa                                >  2:2310328/1‑100 (MQ=255)
                                                                 cGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACa                                <  2:1672489/100‑1 (MQ=255)
                                                                 cGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACa                                >  1:2167169/1‑100 (MQ=255)
                                                                   aGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAgg                              <  2:2062010/100‑1 (MQ=255)
                                                                   aGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAgg                              <  1:516828/100‑1 (MQ=255)
                                                                    gTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGt                             >  2:897220/1‑100 (MQ=255)
                                                                    gTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGt                             <  1:1207963/100‑1 (MQ=255)
                                                                     ttGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTc                            <  2:448490/100‑1 (MQ=255)
                                                                       ggCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CT‑CCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTccac                         >  1:2590072/1‑100 (MQ=255)
                                                                        gCTTTCGCGTTCCTCAGGTGAAGT‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGt                             <  1:1270408/96‑1 (MQ=255)
                                                                         cTTTCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCTGCTCGTTTCCGCGCAACAGGTccacc                        >  1:43997/1‑100 (MQ=255)
                                                                          tttCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTccacca                       <  2:1139001/100‑1 (MQ=255)
                                                                          tttCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTccacca                       <  1:2111424/100‑1 (MQ=255)
                                                                            tCGCGTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATc                     >  2:360568/1‑100 (MQ=255)
                                                                             cgcgTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCa                    >  2:1582614/1‑100 (MQ=255)
                                                                              gcgTTCCTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCAc                   >  1:1770019/1‑100 (MQ=255)
                                                                                   ccTCAGGTGACGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGGTTCGGCGAAACAGGTCCACCATCCcgcgat              >  1:1574897/1‑98 (MQ=255)
                                                                                   ccTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGtt              >  2:1890439/1‑100 (MQ=255)
                                                                                    cTCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGttt             >  2:1695877/1‑100 (MQ=255)
                                                                                    cTCAGGTGAAGG‑CG‑CT‑CCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTg            >  1:445365/1‑100 (MQ=255)
                                                                                    cTCAGGTGAAGG‑CG‑CT‑CCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTg            >  1:1850159/1‑100 (MQ=255)
                                                                                     tCAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTg            <  1:1453216/100‑1 (MQ=255)
                                                                                      cAGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACACTTTTTCCAGCGAGTTTGGGCGAACAAGGTCCACCATCACGCGttt             >  1:1776580/1‑98 (MQ=255)
                                                                                       aGGTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGc          >  2:2103722/1‑100 (MQ=255)
                                                                                        ggTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCg         >  2:1794362/1‑100 (MQ=255)
                                                                                        ggTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCg         >  2:1039138/1‑100 (MQ=255)
                                                                                         gTGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGa        <  1:202827/100‑1 (MQ=255)
                                                                                          tGAAGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGAt       <  1:1454492/100‑1 (MQ=255)
                                                                                          tGAAGG‑CG‑CT‑CCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATa      >  1:720549/1‑100 (MQ=255)
                                                                                            aaGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAg     <  2:857773/100‑1 (MQ=255)
                                                                                             aGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGt    >  2:105025/1‑100 (MQ=255)
                                                                                             aGG‑CG‑CTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGt    >  2:108532/1‑100 (MQ=255)
                                                                                              gg‑CG‑CT‑CCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGc  >  2:2578728/1‑100 (MQ=255)
                                                                                               g‑cg‑cTCCCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGc  <  2:2005364/100‑1 (MQ=255)
                                                                                                      |                                                                                                
CAGAAACGCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGG‑CG‑CT‑CCCC‑C‑CTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGC  >  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: ‑