Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,639,601 C→A 6.5% A231A (GCC→GCA PP_4950 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,639,6010CA 6.5% 248.4 / 5.8 77A231A (GCC→GCAPP_4950hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (38/34);  new base A (2/3);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.30e-01

GCCAGGCCTGGAGTGCCGACGACTTCACCGCACGCCTGCGCCAAGCACGAGAGCAGTTGGGTTTCCTCGGCCGCCCGGCAGTTACCCTCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGCCCGCTGGTCAGTTTTACCGAGCAGGTCAGCGGCTCGCTTAGCCCGGCATTTTCCG  >  NC_002947/5639459‑5639749
                                                                                                                                              |                                                                                                                                                    
gCCAGGCCTGGAGTGCCGACGACTTCACCGCACGCCTGCGCCAAGCACGAGAGCAGTTGGGTTTCCTCGGCCGCCCGGCAGTTACCCTCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCAtg                                                                                                                                                <  4:253980/149‑1 (MQ=255)
    ggCCTGGAGTGCCGACGACTTCACCGCACGCCTGCGCCAAGCACGAGAGCAGTTGGGTTTCCTCGGCCGCCCGGCAGTTACCCTCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCAtgttgt                                                                                                                                            <  1:175619/149‑1 (MQ=255)
    ggCCTGGAGTGCCGACGACTTCACCGCACGCCTGCGCCAAGCACGAGAGCAGTTGGGTTTCCTCGGCCGCCCGGCAGTTACCCTCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCAtgttgt                                                                                                                                            <  4:23660/149‑1 (MQ=255)
               ccGACGACTTAACCGCACGCCTGCGCCAAGCACGAGAGCAGTTGGGTTTCCTCGGCCGCCCGGCAGTTACCCTCAAGCCCGGCACCAACAGCGCCTACCTAGCACCTTCGGCCATGGACGAAGACGCCGGCATGTTGTGCTGGggcgg                                                                                                                                  >  8:131820/1‑148 (MQ=255)
                     acTTCACCGCACGCCTGCTCCAAGCACGAGAGCAGTTGGGTTTCCTCGGCCGCCCGGCAGTTACCCTCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGACGGCATGTTGTGCTGGGGCGGCTTTTcc                                                                                                                           <  7:131820/149‑1 (MQ=255)
                       ttCACCGCACGCCTGCGCCAAGCACGAGAGCAGTTGGGTTTCCTCGGCCGCCCGGCAGTTACCCTCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCcgc                                                                                                                         >  7:302419/1‑149 (MQ=255)
                       ttCACCGCACGCCTGCGCCAAGCACGAGAGCAGTTGGGTTTCCTCGGCCGCCCGGCAGTTACCCTCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCcgc                                                                                                                         >  5:293696/1‑149 (MQ=255)
                       ttCACCGCACGCCTGCGCCAAGCACGAGAGCAGTTGGGTTTCCTCGGCCGCCCGGCAGTTACCCTCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCcgc                                                                                                                         >  8:185038/1‑149 (MQ=255)
                                 gCCTGCGCCAAGCACGAGAGCAGTTGGGTTTCCTCGGCCGCCCGGCAGTTACCCTCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGATGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTg                                                                                                               >  2:158739/1‑149 (MQ=255)
                                              aCGAGAGCAGTTGGGTTTCCTCGGCCGCCCGGCAGTTACCCTCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTcc                                                                                                                           <  8:224796/124‑1 (MQ=255)
                                              aCGAGAGCAGTTGGGTTTCCTCGGCCGCCCGGCAGTTACCCTCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTcc                                                                                                                           >  7:224796/1‑124 (MQ=255)
                                              aCGAGAGCAGTTGGGTTTCCTCGGACGCCCGGCAGTTACCCTCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTgg                                                                                                                                       <  1:344049/112‑1 (MQ=255)
                                              aCGAGAGCAGTTGGGTTTCCTCGGACGCCCGGCAGTTACCCTCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTgg                                                                                                                                       >  2:344049/1‑112 (MQ=255)
                                                 agagCAGTTGGGTTTCCTCGGCCGCCCGGCAGTTACCCTCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGggc                                                                                                                                    >  5:407958/1‑112 (MQ=255)
                                                 agagCAGTTGGGTTTCCTCGGCCGCCCGGCAGTTACCCTCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGggc                                                                                                                                    <  6:407958/112‑1 (MQ=255)
                                                   agCAGTTGGGTTTCCTCGGCCGCCCGGCAGTTACCCTCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGct                                                                                             >  8:203793/1‑149 (MQ=255)
                                                     cAGTTGGGTTTCCTCGGCCGCCCGGCAGTTACCCTCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGctct                                                                                           <  4:149418/149‑1 (MQ=255)
                                                      aGTTGGGTTTCCTCGGCCGCCCGGCAGTTACCCTCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTg                                                                                          <  2:365361/149‑1 (MQ=255)
                                                          gggTTTCCTCGGCCGCCCGGCAGTTACCCTCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAgc                                                                                      >  6:391342/1‑149 (MQ=255)
                                                          gggTTTCCTCGGCCGCCCGGCAGTTACCCTCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAgc                                                                                      >  2:212962/1‑149 (MQ=255)
                                                             tttCCTCGGCCGCCCGGCAGTTACCCTCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCTTGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGcc                                                                                   <  2:310444/149‑1 (MQ=255)
                                                               tCCTCGGCCGCCCGGCAGTTACCCTCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGGCACTGGCAACAGCGCTCTGCAGCGCCTg                                                                                 >  7:167623/1‑149 (MQ=255)
                                                                    ggCCGCCCGGCAGTTACCCTCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTAc                                                                              <  7:203793/147‑1 (MQ=255)
                                                                       cgccCGGCAGTTACCCTCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACgg                                                                         >  2:273222/1‑149 (MQ=255)
                                                                       cgccCGGCAGTTACCCTCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACgg                                                                         >  2:267343/1‑149 (MQ=255)
                                                                             gCAGTTACCCTCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGGGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAAGAGCGCTCTGCAGCGCCTGTACAACGGCGATGc                                                                   <  1:267343/149‑1 (MQ=255)
                                                                             gCAGTTACCCTCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGc                                                                   <  1:273222/149‑1 (MQ=255)
                                                                               agctaaccTGAAGCCCTGCACCTACCGCGCCTAGCTTGCAACTGCGGCAATGGACGTAGTCGCAGGCATGTTGTGCTGGGTCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAAGAGCGCGCTGCAGCTCCTGTACAACGGCGATTCAc                                                                 <  8:177764/143‑1 (MQ=255)
                                                                               aGTTACCCTCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCAc                                                                 >  7:393132/1‑149 (MQ=255)
                                                                                cttACCCTCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATg                                                                    >  6:272651/2‑145 (MQ=255)
                                                                                cttACCCTCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATg                                                                    <  5:272651/144‑1 (MQ=255)
                                                                                   aCCCTCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCt                                                             >  7:96808/1‑149 (MQ=255)
                                                                                     ccTCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTa                                                           <  7:185038/149‑1 (MQ=255)
                                                                                     ccTCAAGCCCGGCACATACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGAGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCt                                                             >  6:89025/1‑147 (MQ=255)
                                                                                       tCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTgg                                                                                                       <  2:252279/103‑1 (MQ=255)
                                                                                       tCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTgg                                                                                                       >  1:252279/1‑103 (MQ=255)
                                                                                       tCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGc                                                         >  5:278332/1‑149 (MQ=255)
                                                                                       tCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGc                                                         >  6:114379/1‑149 (MQ=255)
                                                                                       tCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGc                                                         >  6:203508/1‑149 (MQ=255)
                                                                                       tCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCAGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGc                                                         >  7:28348/1‑149 (MQ=255)
                                                                                       tCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCAGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGc                                                         >  7:34606/1‑149 (MQ=255)
                                                                                         aaGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGc                                                         <  8:302419/147‑1 (MQ=255)
                                                                                          aGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGCCCg                                                      >  8:340643/1‑149 (MQ=255)
                                                                                             ccGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGCCCGCTg                                                   >  5:409259/1‑149 (MQ=255)
                                                                                                 cACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACAAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGCCCGCTGGt                                                 >  6:51658/1‑147 (MQ=255)
                                                                                                     tACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGCCCGCTGGTCAGttt                                           >  3:367772/1‑149 (MQ=255)
                                                                                                       ccGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGCCCGCTGGTCAGTTTTa                                         <  5:391342/149‑1 (MQ=255)
                                                                                                       ccGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGCCCGCTGGTCAGTTTTa                                         <  6:293696/149‑1 (MQ=255)
                                                                                                        cgcgCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGCCCGCTGGTCAGTTTTAc                                        <  1:21072/149‑1 (MQ=255)
                                                                                                        cgcgCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGCCCGCTGGTCAGTTTTAc                                        <  1:212962/149‑1 (MQ=255)
                                                                                                          cgcCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGCCCGCTGGTCAGTTTTAcc                                       >  5:152412/1‑148 (MQ=255)
                                                                                                           gcCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCtttt                                                                                                                             >  3:199224/1‑61 (MQ=255)
                                                                                                           gcCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCtttt                                                                                                                             <  4:199224/61‑1 (MQ=255)
                                                                                                               aCCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACgg                                                               >  6:247624/1‑119 (MQ=255)
                                                                                                               aCCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACgg                                                               <  5:247624/119‑1 (MQ=255)
                                                                                                                ccTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGCCCGCTGGTCAGTTTTACCGAGCAg                                 <  8:96808/148‑1 (MQ=255)
                                                                                                                 cTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGCCCGCTGGTCAGTTTTACCGAGCAgg                                <  7:340643/148‑1 (MQ=255)
                                                                                                                 cTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGCCCGCTGGTCAGTTTTACCGAGCAGGt                               >  8:166647/1‑149 (MQ=255)
                                                                                                                    gCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGCCCGCTGGTCAGTTTTACCGAGCa                                  <  2:92844/143‑1 (MQ=255)
                                                                                                                    gCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGCCCGCTGGTCAGTTTTACCGAGCa                                  >  1:92844/1‑143 (MQ=255)
                                                                                                                     cACCTGCGGCCATGGACGAAGTCTCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGCCCGCTGGTCAGTTTTACCGAGCAGGTCAGc                           <  1:27939/149‑1 (MQ=255)
                                                                                                                      aCCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTTGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGCCCGCTGGTCAGTTTTACCGAGCAGGTCAGCg                          <  8:393132/149‑1 (MQ=255)
                                                                                                                       ccTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGCCCGCTGGTCAGTTTTACCGAGCAGGTCAGCgg                         >  5:289996/1‑149 (MQ=255)
                                                                                                                        cTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGCCCGCTGGTCAGTTTTACCGAGCAg                                 >  4:115466/1‑140 (MQ=255)
                                                                                                                        cTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGCCCGCTGGTCAGTTTTACCGAGCAg                                 <  3:115466/140‑1 (MQ=255)
                                                                                                                          gCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGCCCGCTGGTCAGTTTTACCGAGCAGGTCAGCGGCTc                      >  6:65948/1‑149 (MQ=255)
                                                                                                                            ggCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGCCCGCTGGTCAGTTTTACCGAGCAGGTCAGCGGCTCGc                    <  6:409259/149‑1 (MQ=255)
                                                                                                                              ccATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGAACGGCTTAGCCCGCTGGTCAGTTTTACCGAGCAGGTCAGCGGCTCGCt                   >  1:214453/1‑148 (MQ=255)
                                                                                                                                aTGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGCCCGCTGGTCAGTTTTACCGAGCAGGTCAGCGGCTCGCTTAg                >  6:27683/1‑149 (MQ=255)
                                                                                                                                aTGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGCCCGCTGGTCAGTTTTACCGAGCAGGTCAGCGGCTCGCTTAg                >  2:100754/1‑149 (MQ=255)
                                                                                                                                aTGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTAGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGCCCGCTGGTCAGTTTTACCGAGCAGGTCAGCGGCTCGCTTAg                >  3:323304/1‑149 (MQ=255)
                                                                                                                                   gACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGCCCGCTGGTCAGTTTTACCGAGCAGGTCAGCGGCTCGCTTAGcc              <  6:289996/148‑1 (MQ=255)
                                                                                                                                         gTCGCAGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGCCCGCTGGTCAGTTTTACCGAGCAGGTCAGCGGCTCGCTTAGCCCGGCAtt       <  8:28348/149‑1 (MQ=255)
                                                                                                                                         gTCGCAGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGCCCGCTGGTCAGTTTTACCGAGCAGGTCAGCGGCTCGCTTAGCCCGGCAtt       <  8:34606/149‑1 (MQ=255)
                                                                                                                                            gCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTTGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGCCCGCTGGTCAGTTTTACCGAGCAGGTCAGCGGCTCGCTTAGCCCGGCATTTTc    <  4:323304/149‑1 (MQ=255)
                                                                                                                                            gCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGCCCGCTGGTCAGTTTTACCGAGCAGGTCAGCGGCTCGCTTAGCCCGGCATTTTc    <  1:100754/149‑1 (MQ=255)
                                                                                                                                            gCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGCCCGCTGGTCAGTTTTACCAAGCAGGTCAGCGGCTCGCTTAGCCCGGCAtttt     <  4:367772/148‑1 (MQ=255)
                                                                                                                                             ccGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGCCCGCTGGTCAGTTTTACCGAGCAGGTCAGCGGCTCGCTTAGCCCGGCATTTTcc   >  1:408045/1‑149 (MQ=255)
                                                                                                                                             ccGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGCCCGCTGGTCAGTTTTACCGAGCAGGTCAGCGGCTCGCTTAGCCCGGCATTTTcc   >  1:401765/1‑149 (MQ=255)
                                                                                                                                             ccGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGCCCGCTGGTCAGTTTTACCGAGCAGGTCAGCGGCTCGCTTAGCCCGGCATTTTcc   <  5:89025/149‑1 (MQ=255)
                                                                                                                                              cGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTTGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGCCCGCTGGTCAGTTTTACCGAGCAGGTCAGCGGCTCGCTTAGCCCGGCATTTTCCg  <  5:65948/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
GCCAGGCCTGGAGTGCCGACGACTTCACCGCACGCCTGCGCCAAGCACGAGAGCAGTTGGGTTTCCTCGGCCGCCCGGCAGTTACCCTCAAGCCCGGCACCTACCGCGCCTACCTTGCACCTGCGGCCATGGACGAAGTCGCCGGCATGTTGTGCTGGGGCGGCTTTTCCGCGCAGGCCCTGGCCACTGGCAACAGCGCTCTGCAGCGCCTGTACAACGGCGATGCACGGCTTAGCCCGCTGGTCAGTTTTACCGAGCAGGTCAGCGGCTCGCTTAGCCCGGCATTTTCCG  >  NC_002947/5639459‑5639749

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: