Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,622,913 C→T 6.7% L148L (TTG→TTA) ‡ PP_1422 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,622,9130CT 6.7% 220.5 / 3.1 74L148L (TTG→TTA) ‡PP_1422hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (33/36);  new base T (2/3);  total (35/39)
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.03e-01

GCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGATCACCAGGC  >  NC_002947/1622773‑1623061
                                                                                                                                             |                                                                                                                                                    
gCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGc                                                                                                                                              <  1:16613/149‑1 (MQ=255)
   ttGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTg                                                                                                                                           >  7:11875/1‑149 (MQ=255)
      aGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCTTGGGTGGGCTGGt                                                                                                                                         <  6:45559/148‑1 (MQ=255)
       ggATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCTTGGGTGGGCTGGTca                                                                                                                                       <  5:55608/149‑2 (MQ=255)
        gATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGGCCACCTTGGGTGGGCTGGTca                                                                                                                                       <  4:361815/148‑2 (MQ=255)
         aTCTTGTCCCCCAGCACCGGCTGGAGTTGATCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACAGG‑CAGCGCTCCAATCTGGTTGACCAGCCCACCCAAGCTGGCCTGGACggg                                                                                                                                     >  7:182831/1‑149 (MQ=255)
          gcTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGt                                                                                                                                    <  6:396404/148‑1 (MQ=255)
              gTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGttt                                                                                                                                  >  3:45596/1‑147 (MQ=255)
               tCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTga                                                                                                                                <  7:76718/148‑1 (MQ=255)
                 ccccAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGagag                                                                                                                             >  4:399297/1‑149 (MQ=255)
                 ccccAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGagag                                                                                                                             <  4:45596/149‑1 (MQ=255)
                   ccAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGgt                                                                                                                           <  6:311604/149‑1 (MQ=255)
                    gaGCACGGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTTCTCTTCGTCATAGATGAGTACCGT‑CAGCGCGCTAATGTGGTTGACCAGCCCACCCATGTTGGCCTGGTCGGGTTTGAGAGgtg                                                                                                                          <  8:182831/148‑1 (MQ=255)
                     aGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGgtgt                                                                                                                         <  2:72879/149‑1 (MQ=255)
                             cTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGAc                                                                                                                 >  2:223326/1‑149 (MQ=255)
                                gAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACgg                                                                                                               <  8:11875/148‑1 (MQ=255)
                                gAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACgg                                                                                                               <  4:113831/148‑1 (MQ=255)
                                gAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACg                                                                                                                <  6:192350/147‑1 (MQ=255)
                                       tCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTc                                                                                                       >  1:257868/1‑149 (MQ=255)
                                         gATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGGTTGAGAGGTGTGAAGCGCGGGCAGATTTCGt                                                                                                     >  4:235868/1‑149 (MQ=255)
                                          aTCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTa                                                                                                    <  2:257868/149‑1 (MQ=255)
                                           tCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAg                                                                                                   >  4:369878/1‑149 (MQ=255)
                                           tCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAg                                                                                                   >  1:220633/1‑149 (MQ=255)
                                              caTAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGt                                                                                                     >  4:251545/1‑144 (MQ=255)
                                              caTAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGt                                                                                                     <  3:251545/144‑1 (MQ=255)
                                                 aaTGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATg                                                                                             >  8:276180/1‑149 (MQ=255)
                                                    gCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGccc                                                                                          >  3:207642/1‑149 (MQ=255)
                                                        gCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCa                                                                                      >  4:101745/1‑149 (MQ=255)
                                                         cAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCAc                                                                                     >  1:164525/1‑149 (MQ=255)
                                                              cGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCa                                                                                <  3:101745/149‑1 (MQ=255)
                                                              cGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCa                                                                                <  8:346240/149‑1 (MQ=255)
                                                                    ggACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTga                                                                                                                                >  7:311315/1‑95 (MQ=255)
                                                                    ggACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTga                                                                                                                                <  8:311315/95‑1 (MQ=255)
                                                                      acaacaGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGcc                                                                                           <  2:153158/130‑1 (MQ=255)
                                                                      acaacaGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGcc                                                                                           >  1:153158/1‑130 (MQ=255)
                                                                      acaacaGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTg                                                                        <  3:235868/149‑1 (MQ=255)
                                                                          caGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGcc                                                                     >  6:174329/1‑148 (MQ=255)
                                                                          caGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGcc                                                                     >  5:359836/1‑148 (MQ=255)
                                                                          caGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGcc                                                                     >  2:324373/1‑148 (MQ=255)
                                                                           aGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGAGCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCa                                                                   <  6:134911/149‑1 (MQ=255)
                                                                            gTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGTCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGAt                                                                                              >  1:54073/1‑122 (MQ=255)
                                                                            gTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGAt                                                                                              <  2:54073/121‑1 (MQ=255)
                                                                                    tgCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGTTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATcgcg                                                          <  2:164525/149‑1 (MQ=255)
                                                                                      ctctTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTATTGGCCCAGATATcgcg                                                          >  5:33281/1‑147 (MQ=255)
                                                                                           cGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGttt                                                    <  7:276180/148‑1 (MQ=255)
                                                                                            ggCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTc                                                   <  6:189266/148‑1 (MQ=255)
                                                                                                tAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTc                                              <  1:324373/149‑1 (MQ=255)
                                                                                                    tGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACt                                          >  6:107288/1‑149 (MQ=255)
                                                                                                    tGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACt                                          >  7:232116/1‑149 (MQ=255)
                                                                                                       gTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGcc                                       <  3:399297/149‑1 (MQ=255)
                                                                                                        tACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGat                                                                >  6:134995/1‑123 (MQ=255)
                                                                                                        tACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGat                                                                <  5:134995/123‑1 (MQ=255)
                                                                                                         aCCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTATTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCg                                      <  6:33281/148‑1 (MQ=255)
                                                                                                          ccGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCt                                                  >  1:211283/1‑135 (MQ=255)
                                                                                                          ccGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCt                                                  <  2:211283/135‑1 (MQ=255)
                                                                                                          ccGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGc                                    >  7:128509/1‑149 (MQ=255)
                                                                                                                aGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCg                               >  3:367404/1‑149 (MQ=255)
                                                                                                                 gcgcgcCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCgg                              <  2:220633/149‑1 (MQ=255)
                                                                                                                   gcgcCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGATGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCGGTTTCTGCTCCACTGCCGGCGAGCGGGc                            <  8:74366/149‑1 (MQ=255)
                                                                                                                   gcgcCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCggg                             >  5:331919/1‑148 (MQ=255)
                                                                                                                    cgcCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCa                           <  3:369878/149‑1 (MQ=255)
                                                                                                                     gcCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCAt                          >  8:381755/1‑149 (MQ=255)
                                                                                                                     gcCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCAt                          <  5:107288/149‑1 (MQ=255)
                                                                                                                         aTCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTTAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGa                       >  7:83309/1‑148 (MQ=255)
                                                                                                                          tCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCTCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGaa                      >  6:235046/1‑148 (MQ=255)
                                                                                                                              gTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTgg                                                                       <  8:382549/95‑1 (MQ=255)
                                                                                                                              gTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTgg                                                                       >  7:382549/1‑95 (MQ=255)
                                                                                                                              gTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATc                  <  6:331919/148‑1 (MQ=255)
                                                                                                                              gTTGACCAGCCAACACAAGGTGGCCTGGTCTGGTTTGAGAGGTGTGAAGCGAAGGCAGATTTCGTAGCTGATGCCCTCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGATCCACTGCAGGCGAGAGGGCATGGAAActct                 >  8:287940/1‑149 (MQ=255)
                                                                                                                               ttGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTg                <  8:128509/149‑1 (MQ=255)
                                                                                                                               ccgACCAGGCCACCTTGGGTGGGCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTg                >  1:25032/3‑149 (MQ=255)
                                                                                                                               ccgACCAGGCCACCTTGGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTg                >  8:98499/3‑149 (MQ=255)
                                                                                                                                      gCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGATc         <  6:359836/149‑1 (MQ=255)
                                                                                                                                      gCCCACACAAGGTGGCCTGGTCGGTTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGATc         <  7:287940/149‑1 (MQ=255)
                                                                                                                                          acccaAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGATCAc       >  8:286549/1‑147 (MQ=255)
                                                                                                                                          acccaAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGATCAc       >  8:159401/1‑147 (MQ=255)
                                                                                                                                          acccaAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGAATGGAAAATCTGCTCGATCAc       >  7:179634/1‑147 (MQ=255)
                                                                                                                                           cccaAGGTGGCGTGGTGGGGTTTGAGAGGTGTGAAGCGACGTCAGATTTCGTAGATGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTATCTGCTACACTGCAGGCGAGATGGCATGGAAAATCTGCTCGATCACCAg    >  6:348612/1‑149 (MQ=255)
                                                                                                                                            ccaAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTCCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGATCACCAgg   >  4:163656/1‑149 (MQ=255)
                                                                                                                                             caAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGACATGGAAAATCTGCTCGATCACCAGGc  <  1:223326/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
GCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGATCACCAGGC  >  NC_002947/1622773‑1623061

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