Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,403,814 A→G 5.6% S78G (AGC→GGC)  PP_3014 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,403,8140AG 5.6% 205.4 / 2.4 71S78G (AGC→GGC) PP_3014hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (35/32);  new base G (2/2);  total (37/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 = 8.16e-01

GCAACCCCCGGCAAACCCACTGCCAGCCAACTGACCAGGAAATCGCCGCAGCAGTACGCCGCCTGTGTCCTGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTT  >  NC_002947/3403670‑3403962
                                                                                                                                                |                                                                                                                                                    
gcAACCCCCGGCAAACCCACTGCCAGCCAACTGACCAGGAAATCGCCGCAGCAGTACGCCGCCTGTGTCCTGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAgcgc                                                                                                                                                  >  8:214757/1‑149 (MQ=255)
gcAACCCCCGGCAAACCCACTGCCAGCCAACTGACCAGGAAATCGCCGCAGCAGTACGCCGCCTGTGTCCTGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAgcgc                                                                                                                                                  >  4:31559/1‑149 (MQ=255)
gcAACCCCCGGCAAACCCACTGCCAGCCAACTGACCAGGAAATCGCCGCAGCAGTACGCCGCCTGTGTCCTGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAgcgc                                                                                                                                                  >  8:151666/1‑149 (MQ=255)
 cAACCCCCGGCAAACCCACTGCCAGCCAACTGACCAGGAAATCGCCGCAGCAGTACGCCGCCTGTGTCCTGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGcc                                                                                                                                                 >  4:41716/1‑149 (MQ=255)
 cAACCCCCGGCAAACCCACTGCCAGCCAACTGACCAGGAAATCGCCGCAGCAGTACGCCGCCTGTGTCCTGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGcc                                                                                                                                                 >  5:152940/1‑149 (MQ=255)
  aaCCCCCGGCAAACCCACTGCCAGCCAACTGACCAGGAAATCGCCGCAGCAGTACGCCGCCTGTGTCCTGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCg                                                                                                                                                >  2:202083/1‑149 (MQ=255)
     ccccGGCAAACCCACTGCCAGCCAACTGACCAGGAAATCGCCGCAGCAGTACGCCGCCTGTGTCCTGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCTTGGCGCCGt                                                                                                                                               >  2:59140/1‑147 (MQ=255)
     ccccGGCAAACCCACTGCCAGCCAACTGACCAGGAAATCGCCGCAGCAGTACGCCGCCTGTGTCCTGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCTTGGCGCCGTga                                                                                                                                             <  1:59140/149‑2 (MQ=255)
     ccccGGCAAACCCACTGCCAGCCAACTGACCAGGAAATCGCCGCAGCAGTACGCCGCCTGTGTCCTGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTgg                                                                                                                                             >  7:277641/1‑149 (MQ=255)
     ccccGGCAAACCCACTGCCAGCCAACTGACCAGGAAATCGCCGCAGCAGTACGCCGCCTGTGTCCTGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCCCCGTgg                                                                                                                                             >  8:170855/1‑149 (MQ=255)
       ccGGCAAACCCACTGCCAGCCAACTGACCAGGAAATCGCCGCAGCAGTACGCCGCCTGTGTCCTGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCTTGGCGCCGTgagc                                                                                                                                           <  8:272247/149‑4 (MQ=255)
       ccGGCAAACCCACTGCCAGCCAACTGACCAGGAAATCGCCGCAGCAGTACGCCGCCTGTGTCCTGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCa                                                                                                                                           <  1:283770/149‑1 (MQ=255)
        cGGCAAACCCACTGCCAGCCAACTGACCAGGAAATCGCCGCAGCAGTACGCCGCCTGTGTCCTGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCaa                                                                                                                                          <  7:172688/149‑1 (MQ=255)
        cGGCAAACCCACTGCCAGCCAACTGACCAGGAAATCGCCGCAGCAGTACGCCGCCTGTGTCCTGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCaa                                                                                                                                          <  3:47554/149‑1 (MQ=255)
                      ccagccaACTGACCAGGAAATCGCCGCAGCAGTACGCCGCCTGTGTCCTGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGc                                                                                                                              <  1:67681/147‑1 (MQ=255)
                        agccaACTGACCAGGAAATCGCCGCAGCAGTACGCCGCCTGTGTCCTGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGa                                                                                                                          <  2:174299/149‑1 (MQ=255)
                          ccaACTGACCAGGAAATCGCCGCAGCAGTACGCCGCCTGTGTCCTGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCt                                                                                                                             <  4:204575/144‑1 (MQ=255)
                          ccaACTGACCAGGAAATCGCCGCAGCAGTACGCCGCCTGTGTCCTGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCt                                                                                                                             >  3:204575/1‑144 (MQ=255)
                            aaCTGACCAGGAAATCGCCGCAGCAGTACGCCGACTGTGTCCTGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGAtgtg                                                                                                                      >  1:44913/1‑149 (MQ=255)
                             aCTGACCAGGAAATCGCCGCAGCAGTACACCGCCTGTGTCCTGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCCGCTCACGGCCCCAAGCGCCGTGGAAAGCGACGATGTGCTTGATGTgg                                                                                                                     >  4:33549/1‑149 (MQ=255)
                               tGACCAGGAAATCGCCGCAGCAGTACGCCGCCTGTGTCCTGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGtt                                                                                                                   >  5:263412/1‑149 (MQ=255)
                                 aCCAGGAAATCGCCGCAGCAGTACGCCGCCTGTGTCCTGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGa                                                                                                                 <  7:214757/149‑1 (MQ=255)
                                   cAGGAAATCGCCGCAGCAGTACGCCGCCTGTGTCCTGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACg                                                                                                               >  4:19328/1‑149 (MQ=255)
                                   cAGGAAATCGCCGCAGCAGTACGCCGCCTGTGTCCTGACGTAATGGCAGGCTCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACg                                                                                                               <  3:252693/149‑1 (MQ=255)
                                        aaTCGCCGCAGCAGTACGCCGCCTGTGTCCTGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCAc                                                                                                          <  5:212939/149‑1 (MQ=255)
                                          tCGCCGCAGCAGTACGCCGCCTGTGTCCTGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGa                                                                                                        <  1:202083/149‑1 (MQ=255)
                                                    agTACGCCGCCTGTGTCCTGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCcgc                                                                                              <  6:263412/149‑1 (MQ=255)
                                                    agTACGCCGCCTGTGTCCTGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCcgc                                                                                              <  7:151666/149‑1 (MQ=255)
                                                      tACGCCGCCTGTGTCCTGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCcgcgc                                                                                            <  3:19328/149‑1 (MQ=255)
                                                      tACGCCGCCTGTGTCCTGCCGAAATGGCAGGCGCTGCCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCcgcgc                                                                                            >  6:117174/1‑149 (MQ=255)
                                                           cgccTGTGTCCTGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGcc                                                                                                >  5:215934/1‑140 (MQ=255)
                                                           cgccTGTGTCCTGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGcc                                                                                                <  6:215934/140‑1 (MQ=255)
                                                            gccTGTGTCCTGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCtt                                                                                      <  6:110601/149‑1 (MQ=255)
                                                              cTGTGTCCTGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCt                                                                                    <  7:104875/149‑1 (MQ=255)
                                                                gtgtCCTGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAggg                                                                                                     <  1:40191/130‑1 (MQ=255)
                                                                gtgtCCTGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAggg                                                                                                     >  2:40191/1‑130 (MQ=255)
                                                                gtgtCCTGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTa                                                                                   >  5:51550/1‑148 (MQ=255)
                                                                    ccTGCCGAAATGGCAGGCGCTGGCGCCACTGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAgg                                                                              >  5:331089/1‑149 (MQ=255)
                                                                     cTGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAgg                                                                              <  6:152940/148‑1 (MQ=255)
                                                                      tGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCcgc                                                                                              >  3:127467/1‑131 (MQ=255)
                                                                      tGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCcgc                                                                                              <  4:127467/131‑1 (MQ=255)
                                                                        ccGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCa                                                                           <  7:231120/148‑1 (MQ=255)
                                                                           aaaTGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCtt                                                                       >  2:187949/1‑149 (MQ=255)
                                                                              tGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTtct                                                                    <  2:206069/149‑1 (MQ=255)
                                                                                gCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTtctgt                                                                  <  3:41716/149‑1 (MQ=255)
                                                                                  aGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTtctgtct                                                                >  3:223703/1‑149 (MQ=255)
                                                                                   ggCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTCCACGGGCAGCTTtctgtctg                                                               >  7:240154/1‑149 (MQ=255)
                                                                                    gcgcTGGCGCCACAGGCCACGCAGTAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTgg                                                              <  3:33549/149‑1 (MQ=255)
                                                                                    gcgcTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTgg                                                              >  4:101634/1‑149 (MQ=255)
                                                                                          gcgcCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACaa                                                        <  1:187949/149‑1 (MQ=255)
                                                                                                     cACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAgg                                             >  2:220418/1‑149 (MQ=255)
                                                                                                          aGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCa                                               <  4:108445/142‑1 (MQ=255)
                                                                                                          aGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCa                                               >  3:108445/1‑142 (MQ=255)
                                                                                                           gAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAG‑TACGCCAGGCTGTCAg                                      >  8:165777/1‑149 (MQ=255)
                                                                                                                      cccACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTgg                                                              >  1:77431/1‑115 (MQ=255)
                                                                                                                      cccACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTgg                                                              <  2:77431/115‑1 (MQ=255)
                                                                                                                      cccACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTg                            <  3:101634/149‑1 (MQ=255)
                                                                                                                       ccACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTg                            >  5:1477/1‑148 (MQ=255)
                                                                                                                        cACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGa                          <  5:117174/149‑1 (MQ=255)
                                                                                                                        cACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGATTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGa                          >  1:251779/1‑149 (MQ=255)
                                                                                                                         aCGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTgg                                                                                                                     <  6:251569/57‑1 (MQ=255)
                                                                                                                         aCGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTgg                                                                                                                     >  5:251569/1‑57 (MQ=255)
                                                                                                                            gCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTg                       >  7:283120/1‑148 (MQ=255)
                                                                                                                            gCTACAGGCTAACGGCTCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGAAGCTTTCTCTCTGCAGACAAGTTACGACAGGCTGTCAGGCAATGCCTGGATTg                       >  5:11858/1‑148 (MQ=255)
                                                                                                                                ccGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAAt                    >  6:244766/1‑147 (MQ=255)
                                                                                                                                  ggCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGt                <  8:277641/149‑1 (MQ=255)
                                                                                                                                    gccACGGCGCTTGGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGa              >  2:257110/3‑149 (MQ=255)
                                                                                                                                    cTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGa              <  8:170275/149‑1 (MQ=255)
                                                                                                                                       aCGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGAccc           >  7:197199/1‑149 (MQ=255)
                                                                                                                                       aCGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGAc             <  2:44913/147‑1 (MQ=255)
                                                                                                                                        cGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCg          >  7:269022/1‑149 (MQ=255)
                                                                                                                                                aGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGtt  <  4:223703/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GCAACCCCCGGCAAACCCACTGCCAGCCAACTGACCAGGAAATCGCCGCAGCAGTACGCCGCCTGTGTCCTGCCGAAATGGCAGGCGCTGGCGCCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTT  >  NC_002947/3403670‑3403962

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: ‑