Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,291,249 T→A 5.7% Y283N (TAT→AAT)  PP_4664 → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,291,2490TA 5.7% 255.0 / 4.7 87Y283N (TAT→AAT) PP_4664LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (43/39);  new base A (2/3);  total (45/42)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.67e-01

GCGGGCGGGGCCGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGCGCGGCCGTTGGGCCTATGGTGAC  >  NC_002947/5291102‑5291393
                                                                                                                                                      |                                                                                                                                                
gcgggcggGGCCGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCt                                                                                                                                                  <  1:81501/148‑1 (MQ=255)
gcgggcggGGCCGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTa                                                                                                                                                 <  2:279683/149‑1 (MQ=255)
 gtggctggGCCGGAATCTGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTAt                                                                                                                                                <  6:62129/147‑1 (MQ=255)
   ggcgggGCCGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGc                                                                                                                                              <  1:92157/149‑1 (MQ=255)
     cgggGCCGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATAGCCCGGCAATGCCg                                                                                                                                            <  4:210851/149‑1 (MQ=255)
       gggCCGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCa                                                                                                                                          <  3:185150/149‑1 (MQ=255)
        ggCCGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCg                                                                                                                                            >  6:354713/1‑146 (MQ=255)
            ggAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATAGCCCGGCAATGCCGc                                                                                                                                           >  4:138312/1‑143 (MQ=255)
            ggAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATAGCCCGGCAATGCCGc                                                                                                                                           <  3:138312/143‑1 (MQ=255)
               aTCGGGGTGCTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGt                                                                                                                                  <  2:102003/149‑1 (MQ=255)
                 cGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAg                                                                                                                                >  2:232743/1‑149 (MQ=255)
                   gggTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGc                                                                                                                               >  8:119895/1‑148 (MQ=255)
                      tggtgCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCAcc                                                                                                                           <  1:11679/149‑1 (MQ=255)
                          gCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTa                                                                                                                       >  7:256269/1‑149 (MQ=255)
                             aCTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCAATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTAcc                                                                                                                     >  1:310180/1‑149 (MQ=255)
                             aCTGGCAGCAGTACAGAGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCC‑‑ATGGCCCAACAGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATACCCGGGATATGCCGCAGAACTGGTAGCCACCTTTACCg                                                                                                                    >  8:116865/1‑149 (MQ=255)
                                 gcagcaGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCAGATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGc                                                                                                                  >  1:305628/1‑149 (MQ=255)
                                 gcagcaGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCt                                                                                                                 >  8:217918/1‑148 (MQ=255)
                                 gcagcaGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCt                                                                                                                 >  4:168983/1‑148 (MQ=255)
                                 gcagcaGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCt                                                                                                                 <  3:61453/148‑1 (MQ=255)
                                   agcagTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGt                                                                                                                                  >  8:302863/1‑129 (MQ=255)
                                   agcagTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGt                                                                                                                                  <  7:302863/129‑1 (MQ=255)
                                   agcagTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGAc                                                                                                              >  7:92693/1‑149 (MQ=255)
                                      agTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGAcgcg                                                                                                           >  5:167574/1‑149 (MQ=255)
                                         aCAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACgcgc                                                                                                          >  6:107941/1‑147 (MQ=255)
                                          cAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCaa                                                                                                        <  7:119895/148‑1 (MQ=255)
                                             cgcTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGgc                                                                                                    <  6:167574/149‑1 (MQ=255)
                                              gcTGGCAAGGCATGGTGCCGTTGGACTGGGTGGCACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCGGGGCTATTCCGCAGAACTTGTAGCCACCTTTACCGGCTGACGCGCAAGGgcg                                                                                                   <  7:116865/149‑1 (MQ=255)
                                                tGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGca                                                                                                 >  6:181356/1‑149 (MQ=255)
                                                  gCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGcaca                                                                                               >  5:113592/1‑149 (MQ=255)
                                                        cATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCC‑‑ATGGGCCAACCGCAGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGAGCACACCACTg                                                                                         >  8:61059/1‑149 (MQ=255)
                                                         ctgctgcCAGTGGCACTGGGTGGCACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACt                                                                                          <  3:352099/147‑1 (MQ=255)
                                                         aTGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGt                                                                                        >  3:357053/1‑149 (MQ=255)
                                                         aTGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGt                                                                                        >  3:153462/1‑149 (MQ=255)
                                                            ctgcCGTTGCACTGGGTGGCACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGa                                                                                     <  6:113592/149‑1 (MQ=255)
                                                               ccGTTGCACTGGGTGGAACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGAAGAACTGGTAGCCACCTTTACAGGCTGACGCGCAATGGAGCACACCACTGTCGAgg                                                                                   >  7:16918/1‑148 (MQ=255)
                                                                     gCACTGGGTGGCACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAAc                                                                              <  1:44443/148‑1 (MQ=255)
                                                                       aCTGGGTGGCACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGc                                                                           >  6:278858/1‑149 (MQ=255)
                                                                           ggTGGCACTGCAGGGGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGc                                                                       >  7:327726/1‑149 (MQ=255)
                                                                           ggTGGCACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTg                                                                        <  6:237116/148‑1 (MQ=255)
                                                                           ggTGGCACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGc                                                                       <  8:327726/149‑1 (MQ=255)
                                                                            gTGGCACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGc                                                                       >  3:356273/1‑148 (MQ=255)
                                                                              ggCACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCa                                                                     >  3:129463/1‑148 (MQ=255)
                                                                                cACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGAc                                                                                                              <  2:263505/105‑1 (MQ=255)
                                                                                cACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCa                                                                     >  8:366004/1‑146 (MQ=255)
                                                                                cACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATa                                                                   <  7:366004/148‑1 (MQ=255)
                                                                                cACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAg                                                                  <  3:42801/149‑1 (MQ=255)
                                                                                 aCTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAgc                                                                 >  5:197558/1‑149 (MQ=255)
                                                                                  cTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACATCACTGTCGAGGCAAACCGCCTGCCATAgcg                                                                <  1:232743/149‑1 (MQ=255)
                                                                                     cAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGcgcg                                                              <  7:108675/148‑1 (MQ=255)
                                                                                       ggAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCggg                                                            >  3:164590/1‑148 (MQ=255)
                                                                                               aTGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCa                                                          <  5:76079/144‑1 (MQ=255)
                                                                                               aTGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCa                                                          >  6:76079/1‑144 (MQ=255)
                                                                                                 gggCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCt                                                   >  3:229408/1‑149 (MQ=255)
                                                                                                   gCCAACCGCCTGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCAGACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTgg                                                 <  3:168983/149‑1 (MQ=255)
                                                                                                   gCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTgg                                                 <  7:217918/149‑1 (MQ=255)
                                                                                                   gCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTgg                                                 <  4:33168/149‑1 (MQ=255)
                                                                                                   gCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTgg                                                 <  4:153462/149‑1 (MQ=255)
                                                                                                    ccAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGa                                                >  4:338667/1‑149 (MQ=255)
                                                                                                    ccAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGa                                                >  2:15408/1‑149 (MQ=255)
                                                                                                     cAACCGCCTGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCTGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGcc                                                                          <  4:345430/122‑1 (MQ=255)
                                                                                                     cAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGAATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGcc                                                                          >  3:345430/1‑122 (MQ=255)
                                                                                                         cgccgGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCg                                           <  5:354713/149‑1 (MQ=255)
                                                                                                           ccgGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAg                                         <  2:305628/149‑1 (MQ=255)
                                                                                                           ccgGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAg                                         >  1:120553/1‑149 (MQ=255)
                                                                                                           ccgGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGACTGCCATAGCGCGGGCAGACACCTGCACGTCGAg                                         >  4:93283/1‑149 (MQ=255)
                                                                                                              gTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGGCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGcc                                       >  1:255025/1‑148 (MQ=255)
                                                                                                                  cTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGc                                  <  4:229408/149‑1 (MQ=255)
                                                                                                                  cTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGc                                  <  1:181933/149‑1 (MQ=255)
                                                                                                                   tGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGc                                                           <  6:55152/123‑1 (MQ=255)
                                                                                                                   tGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGc                                                           >  5:55152/1‑123 (MQ=255)
                                                                                                                   tGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGc                                  >  5:143461/1‑148 (MQ=255)
                                                                                                                       gtgCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTg                             >  6:219677/1‑149 (MQ=255)
                                                                                                                       gtgCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTg                             <  8:256269/149‑1 (MQ=255)
                                                                                                                             aCGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTgcggcgc                       >  5:72582/1‑149 (MQ=255)
                                                                                                                             aCGCGATTATGCGGCATTGCCGGGCTACGCCGCAGAACTGCTAGCCACCTTTACCGGCGGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTgcggcgc                       <  6:300469/149‑1 (MQ=255)
                                                                                                                                 gATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGCGCGGcc                   >  4:43458/1‑149 (MQ=255)
                                                                                                                                     ctGCGGCATAGCCCGGCAATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGCGCGGCCGTTg               >  4:151323/2‑149 (MQ=255)
                                                                                                                                      tGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGCGCGGCCGTTg               <  2:310180/148‑1 (MQ=255)
                                                                                                                                        cGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGCGCGGCCGTTGGGc            <  5:278858/149‑1 (MQ=255)
                                                                                                                                        cGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGCGCGGCCGTTGGGc            <  8:92693/149‑1 (MQ=255)
                                                                                                                                        cGGCACTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGCGCGGCCGTTGGGc            <  2:127002/149‑1 (MQ=255)
                                                                                                                                         ggCATTGCCGGGCTATGCCGCAGAACTGTTAGCCACCTTTACCGGCTGACGCGCAAGGGAGCACACCACTGTCGAGGAAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGAGCGGCCGTTGGGcc           >  7:349601/1‑149 (MQ=255)
                                                                                                                                           cATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGCGCGGCCGTTGGTCCTa         >  6:387529/1‑149 (MQ=255)
                                                                                                                                           cATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGCGCGGCCGTTGGGCCTa         >  3:338326/1‑149 (MQ=255)
                                                                                                                                             ttGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGCGCGGCCGTTGGGCCTATg       <  6:197558/149‑1 (MQ=255)
                                                                                                                                             tagCCCGGCAATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGCGCGGCCGTTGGGCCTATg       <  3:151323/147‑1 (MQ=255)
                                                                                                                                                  gggCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGCGCGCCGGCAACTGCGGCGCGGCCGTTGGGCCTATGGTgcc  >  4:251246/1‑147 (MQ=255)
                                                                                                                                                      |                                                                                                                                                
GCGGGCGGGGCCGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCC‑‑ATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGCGCGGCCGTTGGGCCTATGGTGAC  >  NC_002947/5291102‑5291393

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