Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,819,795:1 +A 7.4% intergenic (+32/‑443) 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,7951.A 7.4% 336.6 / 15.4 81intergenic (+32/‑443)rpoS/PP_1624RNA polymerase sigma S factor/group II intron‑encoding maturase
Reads supporting (aligned to +/- strand):  ref base . (43/32);  new base A (4/2);  total (47/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

CTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCGG  >  NC_002947/1819650‑1819942
                                                                                                                                                   |                                                                                                                                                   
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                                                                                                                           <  1:105640/128‑1 (MQ=255)
                                            ggAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGc                                                                                                                           >  2:105640/1‑128 (MQ=255)
                                            ggAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGCAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGTGCCAGGCATGGCGCGTTGCGCGTATgcg                                                                                                      >  8:334004/1‑149 (MQ=255)
                                               aGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTgcgc                                                                                                                                <  3:67151/120‑1 (MQ=255)
                                               aGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTgcgc                                                                                                                                >  4:67151/1‑120 (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                                                                                                      >  3:301736/1‑137 (MQ=255)
                                                        gcgcTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATgcg                                                                                                      <  4:301736/137‑1 (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                                                                                            >  8:204755/1‑118 (MQ=255)
                                                                                     gTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCtt                                                                                            <  7:204755/118‑1 (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                                                                                                                              <  4:64058/77‑1 (MQ=255)
                                                                                            cAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCa                                                                                                                              >  3:64058/1‑77 (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                        <  5:242099/144‑1 (MQ=255)
                                                                                                                               aaaaaaCGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTa                        >  6:242099/1‑144 (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                      >  4:349836/1‑136 (MQ=255)
                                                                                                                                          ccGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGc                      <  3:349836/136‑1 (MQ=255)
                                                                                                                                          ccGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGc                      <  3:344218/136‑1 (MQ=255)
                                                                                                                                          ccGGTGAGCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGc                      >  4:344218/1‑136 (MQ=255)
                                                                                                                                                  ctcGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCgg  <  6:60767/148‑1 (MQ=34)
                                                                                                                                                   |                                                                                                                                                   
CTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCGG  >  NC_002947/1819650‑1819942

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: