Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,819,789 Δ1 bp 7.3% intergenic (+26/‑449) rpoS → / → PP_1624 RNA polymerase sigma S factor/group II intron‑encoding maturase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,819,7890G. 7.3% 254.7 / 13.6 82intergenic (+26/‑449)rpoS/PP_1624RNA polymerase sigma S factor/group II intron‑encoding maturase
Reads supporting (aligned to +/- strand):  ref base G (43/33);  new base . (4/2);  total (47/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

TTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGG  >  NC_002947/1819645‑1819933
                                                                                                                                                 |                                                                                                                                                 
ttGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTg                                                                                                                                                <  8:285867/148‑1 (MQ=255)
     cTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCg                                                                                                                                         >  3:349599/1‑149 (MQ=255)
      tGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCgg                                                                                                                                        >  6:381672/1‑149 (MQ=255)
        gagaTCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCggg                                                                                                                                       <  5:29096/148‑1 (MQ=255)
        gagaTCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCggg                                                                                                                                       <  4:53309/148‑1 (MQ=255)
          gaTCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCg                                                                                                                                    >  3:325713/1‑149 (MQ=255)
              ggCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGtttt                                                                                                                                >  3:3791/1‑149 (MQ=255)
              ggCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGtatt                                                                                                                                >  8:128334/1‑149 (MQ=255)
               gCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGttttt                                                                                                                               >  4:21036/1‑149 (MQ=255)
                      ccGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCGGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTgcg                                                                                                                        >  2:315804/1‑149 (MQ=255)
                      ccGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTgcg                                                                                                                        >  2:6000/1‑149 (MQ=255)
                                     gCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGTGGCGTTTTTGTGTGCGCCAGGCATGGCGCGt                                                                                                         >  6:50192/1‑149 (MQ=255)
                                      cAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGtt                                                                                                        <  6:52138/149‑1 (MQ=255)
                                        aaTCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTgc                                                                                                      <  4:349599/149‑1 (MQ=255)
                                         aTCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTgcg                                                                                                     >  4:259147/1‑149 (MQ=255)
                                            cAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGcgcg                                                                                                   >  8:327436/1‑148 (MQ=255)
                                               gTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTtgtg                                                                                                                            <  1:7375/120‑1 (MQ=255)
                                               gTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTtgtg                                                                                                                            >  2:7375/1‑120 (MQ=255)
                                                 ggAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGc                                                                                                                  >  2:105640/1‑128 (MQ=255)
                                                 ggAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGc                                                                                                                  <  1:105640/128‑1 (MQ=255)
                                                 ggAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGCAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGTGCCAGGCATGGCGCGTTGCGCGTATgcg                                                                                             >  8:334004/1‑149 (MQ=255)
                                                    aGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTgcgc                                                                                                                       >  4:67151/1‑120 (MQ=255)
                                                    aGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTgcgc                                                                                                                       <  3:67151/120‑1 (MQ=255)
                                                        cTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCgg                                                                                                                                        <  5:321136/99‑1 (MQ=255)
                                                        cTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCgg                                                                                                                                        >  6:321136/1‑99 (MQ=255)
                                                        cTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAAccc                                                                                       >  6:260491/1‑148 (MQ=255)
                                                            aGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGt                                                                                                  <  8:321236/133‑1 (MQ=255)
                                                            aGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGt                                                                                                  >  7:321236/1‑133 (MQ=255)
                                                             gcgcTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATgcg                                                                                             <  4:301736/137‑1 (MQ=255)
                                                             gcgcTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATgcg                                                                                             >  3:301736/1‑137 (MQ=255)
                                                                 ttGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTtgt                                                                             >  7:150607/1‑149 (MQ=255)
                                                                     gcgAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTgtggtg                                                                         <  4:325713/149‑1 (MQ=255)
                                                                       gagaTCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGCGCGCCATGCCTGgcgc                                                                                                           >  7:149757/1‑113 (MQ=255)
                                                                       gagaTCCTTGAAAAGAATGCTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGCGCGCCATGCCTGgcgc                                                                                                           <  8:149757/113‑1 (MQ=255)
                                                                          aTCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCACCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGAGTATGCGCAACCCGCTTGGTTGTGGTGGcc                                                                      >  5:305509/1‑147 (MQ=255)
                                                                                 aaaaGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCGAAGTGGGTg                                                             >  1:331327/1‑149 (MQ=255)
                                                                                       aTGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTgg                                                       <  5:50192/149‑1 (MQ=255)
                                                                                          gTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCtt                                                                                   <  7:204755/118‑1 (MQ=255)
                                                                                          gTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCtt                                                                                   >  8:204755/1‑118 (MQ=255)
                                                                                           tCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGCCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGt                                                   >  4:126785/1‑149 (MQ=255)
                                                                                             tGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGa                                                 <  5:260491/149‑1 (MQ=255)
                                                                                              gTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTtgtg                                                                            <  6:22556/121‑1 (MQ=255)
                                                                                              gTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTtgtg                                                                            >  5:22556/1‑121 (MQ=255)
                                                                                               tCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGaaa                                               <  3:21036/149‑1 (MQ=255)
                                                                                                 cAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCa                                                                                                                     >  3:64058/1‑77 (MQ=255)
                                                                                                 cAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCa                                                                                                                     <  4:64058/77‑1 (MQ=255)
                                                                                                       gTCATTGTTCCGGTAGCTAGCGCTGTAACAAAAAACGCCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCctct                                        <  1:315804/149‑1 (MQ=255)
                                                                                                       gTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGGGAAAGTCCTCTa                                       >  8:352461/1‑149 (MQ=255)
                                                                                                          aTTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTaca                                     >  5:60767/1‑148 (MQ=255)
                                                                                                          aTTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTaca                                     <  4:281085/148‑1 (MQ=255)
                                                                                                          aTTGTTCCAGCAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGTGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTacac                                    <  7:334004/149‑1 (MQ=255)
                                                                                                           ttGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAg                                              >  8:219779/1‑138 (MQ=255)
                                                                                                           ttGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAg                                              <  7:219779/138‑1 (MQ=255)
                                                                                                              ttCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACAcc                                 <  8:150607/148‑1 (MQ=255)
                                                                                                              ttCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACAcc                                 <  3:126785/148‑1 (MQ=255)
                                                                                                                 cAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCgg                              >  8:257796/1‑148 (MQ=255)
                                                                                                                      gctagcGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAggg                        >  6:20921/1‑149 (MQ=255)
                                                                                                                      gctagcGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAgg                         <  4:3791/148‑1 (MQ=255)
                                                                                                                         agcGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAgggg                       >  5:170357/1‑147 (MQ=255)
                                                                                                                            gcTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAgtg                  >  4:122962/1‑149 (MQ=255)
                                                                                                                             cTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGa                                                      <  1:140499/112‑1 (MQ=255)
                                                                                                                             cTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGa                                                      >  2:140499/1‑112 (MQ=255)
                                                                                                                              tGTAACAAAAAACGCCCCGGTGATGTCGGGGCGTTTTTGTGTGGACCAGGCATGGCCCGTTGCGCGTATGCGCAACCCCCTTGGTTGTGGTGGCCAGGTGGGTGACTTGGAGGGGAAGGTCGTCTAAACACCCGGCAAGGGGAAGTGtt                >  1:268433/1‑149 (MQ=17)
                                                                                                                               gcacacaaaaaCGCCCCGAGATCACCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAgtgt                 >  1:216510/7‑147 (MQ=11)
                                                                                                                                tAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGtt                <  7:257796/147‑1 (MQ=255)
                                                                                                                                tAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAg              >  7:80126/1‑149 (MQ=255)
                                                                                                                                  aCAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGcc            >  8:187273/1‑149 (MQ=255)
                                                                                                                                   acaaaaaCGCCCCGAGATCACCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCg           >  2:99097/3‑149 (MQ=11)
                                                                                                                                    caaaaaCGCCCCGAGATCACCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCg           >  1:21828/2‑148 (MQ=11)
                                                                                                                                    aaaaaaCGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTa               >  6:242099/1‑144 (MQ=255)
                                                                                                                                    aaaaaaCGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTa               <  5:242099/144‑1 (MQ=255)
                                                                                                                                    aaaaaaCGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCgg          <  3:122962/149‑1 (MQ=255)
                                                                                                                                    aaaaaaCGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCg           >  4:203044/1‑148 (MQ=255)
                                                                                                                                    aaaaaaCGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTCGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCg           >  4:129774/1‑148 (MQ=255)
                                                                                                                                    aaaaaaCGCCCCGAGATCACCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGaa                     <  7:39758/138‑1 (MQ=11)
                                                                                                                                     aaaaaCGCCCCGAGATCACCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTa               <  2:186063/143‑1 (MQ=11)
                                                                                                                                     aaaaaCGCCCCGAGATCACCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAACTGTTAgaga           >  1:186063/1‑144 (MQ=11)
                                                                                                                                             ccccGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAgg  <  5:381672/149‑1 (MQ=255)
                                                                                                                                               ccGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGc             <  3:349836/136‑1 (MQ=255)
                                                                                                                                               ccGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGc             >  4:349836/1‑136 (MQ=255)
                                                                                                                                               ccGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGc             <  3:344218/136‑1 (MQ=255)
                                                                                                                                               ccGGTGAGCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGc             >  4:344218/1‑136 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGG  >  NC_002947/1819645‑1819933

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: