Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,059,023 G→C 6.2% P54A (CCG→GCG)  fecI ← RNA polymerase sigma 19 factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,059,0230GC 6.2% 249.2 / 4.3 77P54A (CCG→GCG) fecIRNA polymerase sigma 19 factor
Reads supporting (aligned to +/- strand):  ref base G (39/33);  new base C (3/2);  total (42/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGCGACGGCACGTGTTGCTCGGGCAACTGCGCCAACGCATCAAGGTAGGCCTGCTCCAGGCTACGTCGCTGGTGGAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGCCCATACGCCGGCGCAGCCAGGCTTCAAGCCAGCCGCGATGCTCCCGGTACAGGGCAGGCAGGCCAGGCTCCGGTGGCG  >  NC_002947/4058880‑4059166
                                                                                                                                               |                                                                                                                                               
tGGCGACGGCACGTGTTGCTCGGGCAACTGCGCCAACGCATCAAGGTAGGCCTGCTCCAGGCTACGTCGCTGGTGGAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTggg                                                                                                                                            >  2:263802/1‑149 (MQ=255)
  gCGACGGCACGTGTTGCTCGGGCAACTGCGCCAACGCATCAAGGTAGGCCTGCTCCAGGCTACGTCGCTGGTGGAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGgcg                                                                                                                                          >  5:392451/1‑149 (MQ=255)
     aCGGCACGTGTTACTCGGGCAACTGCGCCAACGCATCAAGGTAGGCCTGCTCCAGGCTACGTCGCTGGTGGAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTg                                                                                                                                       <  6:76083/149‑1 (MQ=255)
       ggCACGTGTTGCTCGGGCAACTGCGCCAACGCATCAAGGTAGGCCTGCTCCAGGCTACGTCGCTGGTGGAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTGGCGTATGTCTGCCAGCGGCTGGGCGCTGGc                                                                                                                                     <  1:32225/149‑1 (MQ=255)
        gCACGTGTTGCTCGGGCAACTGCGCCAACGCATCAAGGTAGGCCTGCTCCAGGCTACGTCGCTGGTGGAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCTGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGcc                                                                                                                                    <  1:213083/149‑1 (MQ=255)
            gtgtTGCTCGGGCAACTGCGCCAACGCATCAAGGTAGGCCTGCTCCAGGCTACGTCGCTGGTGGAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGCCCAGCCGCTGGc                                                                                                                                     <  4:196141/144‑1 (MQ=255)
            gtgtTGCTCGGGCAACTGCGCCAACGCATCAAGGTAGGCCTGCTCCAGGCTACGTCGCTGGTGGAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGCCCAGCCGCTGGc                                                                                                                                     >  3:196141/1‑144 (MQ=255)
            gtgtTGCTCGGGCAACTGCGCCAACGCATCAAGGTAGGCCTGCTCCAGGCTACGTCGCTGGTGGAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCACCGCCTGGGCGCTGGCCAGaa                                                                                                                                >  7:145117/1‑149 (MQ=255)
             tgtTGCTCGGGCAACTGCGCCAACGCATCAAGGTAGGCCTGCTCCAGGCTACGTCGCTGGTGGAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAAt                                                                                                                               <  6:348160/149‑1 (MQ=255)
             tgtTGCTCGGGCAACTGCGCCAACGCATCAAGGTAGGCCTGCTCCAGGCTACGTCGCTGGTGGAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAAt                                                                                                                               <  7:52671/149‑1 (MQ=255)
               ttGCTCGGGCAACTGCGCCAACGCATCAAGGTAGGCCTGCTCCAGGCTACGTCGCTGGTGGAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATgc                                                                                                                             <  1:162812/149‑1 (MQ=255)
                  cTCGGGCAACTGCGCCAACGCATCAAGGTAGGCCTGCTCCAGGCTACGTCGCTGGTGGAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATgcgc                                                                                                                           <  1:127807/148‑1 (MQ=255)
                  cTCGGGCAACTGCGCCAACGCATCAAGGTAGGCCTGCTCCAGGCTACGTCGCTGGTGGAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCa                                                                                                                          <  6:392451/149‑1 (MQ=255)
                  cTCGGGCAACTGCGCCAACGCATCAAGGTAGGCCTGCTCCAGGCTACGTCGCTGGTGGAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGCAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATgcgc                                                                                                                           <  1:281796/148‑1 (MQ=255)
                    cGGGCAACTGCGCCAACGCATCAAGGTAGGCCTGCTCCAGGCTACGTCGCTGGTGGAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAg                                                                                                                         >  5:26191/1‑148 (MQ=255)
                     gggCAACTGCGCCAACGCATCAAGGTAGGCCTGCTCCAGGCTACGTCGCTGGTGGAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGa                                                                                                                        >  4:19103/1‑148 (MQ=255)
                       gCAACTGCGCCAACGCATCAAGGTAGGCCTGCTCCAGGCTACGTCGCTGGTGGAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGaaaa                                                                                                                     <  3:194611/149‑1 (MQ=255)
                       gCAACTGCGCCAACGCATCAAGGTAGGCCTGCTCCAGGCTACGTCGCTGGTGGAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGaaaa                                                                                                                     <  3:349062/149‑1 (MQ=255)
                       gCAACTGCGCCAACGCATCAAGGTAGGCCTGCTCCAGGCTACGTCGCTGGTGGAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGaa                                                                                                                       >  5:17239/1‑147 (MQ=255)
                         aaCTGCGCCAACGCATCAAGGTAGGCCTGCTCCAGGCTACGTCGCTGGTGGAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGt                                                                                                                   <  5:37037/149‑1 (MQ=255)
                             gcgcCAACGCATCAAGGTAGGCCTGCTCCAGGCTACGTCGCTGGTGGAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATcc                                                                                                               >  5:85208/1‑149 (MQ=255)
                              cgcCAACGCATCAAGGTAGGCCTGCTCCAGGCTACGTCGCTGGTGGAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATcc                                                                                                               <  6:358622/148‑1 (MQ=255)
                                ccAACGCATCAAGGTAGGCCTGCTCCAGGCTACGTCGCTGGTGGAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTgg                                                                                                            <  3:380122/149‑1 (MQ=255)
                                   aCGCATCAAGGTAGGCCTGCTCCAGGCTACGTCGCTGGTGGAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGCCCAGCCGCTGGc                                                                                                                                     >  7:336431/1‑121 (MQ=255)
                                   aCGCATCAAGGTAGGCCTGCTCCAGGCTACGTCGCTGGTGGAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGCCCAGCCGCTGGc                                                                                                                                     <  8:336431/121‑1 (MQ=255)
                                       aTCAAGGTAGGCCTGCTCCAGGCTACGTCGCTGGTGGAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGt                                                                                                     >  3:109518/1‑149 (MQ=255)
                                       aTCAAGGTAGGCCTGCTCCAGGCTACGTCGCTGGTGGAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGt                                                                                                     >  8:294666/1‑149 (MQ=255)
                                          aaGGTAGGCCTGCTCCAGGCTACGTCGCTGGTGGAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTcgg                                                                                                  <  6:17239/149‑1 (MQ=255)
                                           aGGTAGGCCTGCTCCAGGCTACGTCGCTGGTGGAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTcggc                                                                                                 >  4:66727/1‑149 (MQ=255)
                                            ggTAGGCCTGCTCCAGGCTACGTCGCTGGTGGAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTcggcg                                                                                                <  6:26191/149‑1 (MQ=255)
                                            ggTAGGCCTGCTCCAGGCTACGTCGCTGGTGGAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTcggcg                                                                                                >  8:2162/1‑149 (MQ=255)
                                                ggCCTGCTCCAGGCTACGTCGCTGGTGGAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCAt                                                                                            <  5:204976/149‑1 (MQ=255)
                                                   cTGCTCCAGGCTACGTCGCTGGTGGAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATTCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATccc                                                                                         <  1:263802/149‑1 (MQ=255)
                                                   cTGCTCCAGGCTACGTCGCTGGTGGAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATcc                                                                                          >  4:131709/1‑148 (MQ=255)
                                                      cTCCAGGCTACGTCGCTGGTGGAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAgg                                                                                      >  5:326268/1‑149 (MQ=255)
                                                       tCCAGGCTACGTCGCTGGTGGAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAgg                                                                                      >  4:321890/1‑148 (MQ=255)
                                                              tACGTCGCTGGTGGAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGcc                                                                               >  5:243047/1‑148 (MQ=255)
                                                                      tggtggAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGCCCATACGCCg                                                                      >  8:359062/1‑149 (MQ=255)
                                                                        gtggAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGCCCATACGCCGgc                                                                    >  4:147258/1‑149 (MQ=255)
                                                                         tggAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGCCCATACGCCGgcg                                                                   >  6:274054/1‑149 (MQ=255)
                                                                              gTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGCCCATACGCCGGCGCAGcc                                                              >  8:124089/1‑149 (MQ=255)
                                                                              gTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGCCCATACGCCGGCGCAGcc                                                              <  3:344678/149‑1 (MQ=255)
                                                                                 gCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGCCCATACGCCGGCGCAGCCAgg                                                           >  6:47134/1‑149 (MQ=255)
                                                                                  cTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGCCCATACGCCGGCGCAGCCAGGc                                                          <  7:359062/149‑1 (MQ=255)
                                                                                  cTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGCCCATACGCCGGCGCAGCCAGGc                                                          >  4:132950/1‑149 (MQ=255)
                                                                                   tGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGCCCATACGCCGGCGCAGCCAGGCt                                                         >  8:276131/1‑149 (MQ=255)
                                                                                   tGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATATCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGCCCATACGCCGGCGCAGCCAGGCt                                                         >  1:286210/1‑149 (MQ=255)
                                                                                   tGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGAACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGCCCATACGCCGGCGCAGCCAGGCt                                                         >  3:375425/1‑149 (MQ=255)
                                                                                       cAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGCCCATACGCCGGCGCAGCCAGGCTTCaa                                                     >  5:245144/1‑149 (MQ=255)
                                                                                         aaCGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGCCCATACGCCGGCGCAGCCAGGCTTCAagc                                                   <  5:274054/149‑1 (MQ=255)
                                                                                                  ccACAGTCAGCAGGTAGGTACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGCCCATACGCCGGCGCAGCCAGGCTTCAAGCCAGCCGCGa                                          >  6:223440/1‑149 (MQ=255)
                                                                                                  ccACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGCCCATACGCCGGCGCAGCCAGGCTTCAAGCCAGCCGCGa                                          >  6:277407/1‑149 (MQ=255)
                                                                                                   cacaGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGCCCATACGCCGGCGCAGCCAGGCTTCAAGCCAGCCGCGAt                                         >  6:277102/1‑149 (MQ=255)
                                                                                                   cacaGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGCCCATACGCCGGCGCAGCCAGGCTTCAAGCCAGCCGCGAt                                         >  6:74775/1‑149 (MQ=255)
                                                                                                    acaGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGCCCATACGCCGGCGCAGCCAGGCTTCAAGCCAGCCGCGATg                                        >  7:98819/1‑149 (MQ=255)
                                                                                                          agcagGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGCCCATACGCCGGCGCAGCCAGGCTTCAAGCCAGCCGCGATGCTccc                                   <  5:277407/148‑1 (MQ=255)
                                                                                                           gcagGTAGGCACGCGGTTCGCGTATGTCTGCCTGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGCCCATACGCCGGCGCAGCCAGGCTTCAAGCCAGCCGCGATGCTccctg                                 <  8:145117/149‑3 (MQ=255)
                                                                                                               gTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGCCCATACGCCGGCGCAGCCAGGCTTCAAGCCAGCCGCGATGCTCCCGGTa                               >  8:338412/1‑147 (MQ=255)
                                                                                                                tAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCACAGGCATTGCCCATACGCCGGAGCAGCCAGGCTTCAAGCCAGCCGCGATGCTCCCGGTACAg                            >  4:146923/1‑149 (MQ=255)
                                                                                                                  ggCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGCCCATACGCCGGCGCAGCCAGGCTTCAAGCCAGCCGCGATGCTCCCGGTACAgg                           <  3:321890/148‑1 (MQ=255)
                                                                                                                    cACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGCCCATACGCCGGCGCAGCCAGGCTTCAAGCCAGCCGCGATGCTCCCGGTACAggg                          >  8:47782/1‑147 (MQ=255)
                                                                                                                     aCGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATcc                                                                                          >  2:274664/1‑82 (MQ=255)
                                                                                                                     aCGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATcc                                                                                          <  1:274664/82‑1 (MQ=255)
                                                                                                                      cgcgGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGCCCATACGCCGGCGCAGCCAGGCTTCAAGCCAGCCGCGATGCTCCCGGTACAGggcagg                      <  7:47782/149‑1 (MQ=255)
                                                                                                                       gcgGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGCCCATACGCCGGCGCAGCCAGGCTTCAAGCCAGCCGCGATGCTCCCGGTACAGggcaggc                     <  7:294666/149‑1 (MQ=255)
                                                                                                                         ggTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGCCCATACGCCGGCGCAGCCAGGCTTCAAGCCAGCCGCGATGCTCCCGGTACAGGgcaggcag                   <  5:223440/149‑1 (MQ=255)
                                                                                                                              gcgTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGCCCATACGCCGGCGCAGCCAGGCTTCAAGCCAGCCGCGATGCTCCCGGTACAGGGCAGGcaggccag              >  8:21870/1‑149 (MQ=255)
                                                                                                                                 tATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGCCCATACGCCGGCGCAGCCAGGCTTCAAGCCAGCCGCGATGCTCCCGGTACAGGGCAGGCAGGCCAGGCt           <  5:47134/149‑1 (MQ=255)
                                                                                                                                       tGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGCCCATACGCCGGCGCAGCCAGGCTTCAAGCCAGCCGCGATGCTCCCGGTACAGGGCAGGCAGGCCAGGCTCCGGTg     >  1:366666/1‑149 (MQ=255)
                                                                                                                                       tGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGCCCATACGCCGGCGCAGCCAGGCTTCAAGCCAGCCGCGATGCTCCCGGTACAGGGCAGGCAGGCCAGGCTCCGGTg     >  4:354039/1‑149 (MQ=255)
                                                                                                                                        gCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGCCCATACGCCgg                                                                     <  7:177580/84‑1 (MQ=255)
                                                                                                                                        gCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGCCCATACGCCgg                                                                     >  8:177580/1‑84 (MQ=255)
                                                                                                                                        gCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGCCCATACGCCGGCGCAGCCAGGCTTCAAGCCAGCCGCGATGCTCCCGGTACAGGGCAGGCAGGCCAGGCTCCGGTgg    >  2:19556/1‑149 (MQ=255)
                                                                                                                                         ccAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGCCCATACGCCGGCGCAGCCAGGCTTCAAGCCAGCCGCGATGCTCCCGGTACAGGGCAGGCAGGCCAGGCTCCGGt      <  8:127676/146‑1 (MQ=255)
                                                                                                                                         ccAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGCCCATACGCCGGCGCAGCCAGGCTTCAAGCCAGCCGCGATGCTCCCGGTACAGGGCAGGCAGGCCAGGCTCCGGt      >  7:127676/1‑146 (MQ=255)
                                                                                                                                         ccAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGCCCATACGCCGGCGCAGCCAGGCTTCAAGCCAGCCGCGATGCTCCCGGTACAGGGCAGGCAGGCCAGGCTCCGGt      <  7:124089/146‑1 (MQ=255)
                                                                                                                                          cAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGCCCATACGCCGGCGCAGCCAGGCTTCAAGCCAGCCGCGATGCTCCCGGTACAGGGCAGGCAGGCCAGGCTCCGGTggcg  <  3:66727/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                               
TGGCGACGGCACGTGTTGCTCGGGCAACTGCGCCAACGCATCAAGGTAGGCCTGCTCCAGGCTACGTCGCTGGTGGAAGTTGCTGAGCAAACGCTTGCCCACAGTCAGCAGGTAGGCACGCGGTTCGCGTATGTCTGCCAGCGGCTGGGCGCTGGCCAGAATGCGCAGAAAAGTATCCTGGCTGAGGTCGGCGGCATCCCAGGCATTGCCCATACGCCGGCGCAGCCAGGCTTCAAGCCAGCCGCGATGCTCCCGGTACAGGGCAGGCAGGCCAGGCTCCGGTGGCG  >  NC_002947/4058880‑4059166

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