Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 943,611 C→G 5.1% G671R (GGG→CGG)  PP_0806 ← surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947943,6110CG 5.1% 294.7 / 2.7 92G671R (GGG→CGG) PP_0806surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (35/52);  new base G (3/2);  total (38/54)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

GGCGACGGTGCCGGTCGCGAGGATGTTGCCTGCCGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCGGGTTGAGGTTGAGGCTGAACACGCCGGTCGGGCCGACCACGGCAGTGCCCAGT  >  NC_002947/943468‑943759
                                                                                                                                                 |                                                                                                                                                    
ggcgACGGTGCCGGTCGCGAGGATGTTGCCTGCCGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGGCCCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCg                                                                                                                                                  <  1:100149/149‑1 (MQ=255)
     cGGTGCCGGTCGCGAGGATGTTGCCTGCCGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATcc                                                                                                                                            <  2:16654/149‑1 (MQ=255)
     cGGTGCCGGTCGCGAGGATGTTGCCTGCCGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATcc                                                                                                                                            <  4:311815/149‑1 (MQ=255)
      ggTGCCGGTCGCGAGGATGTTGCCTGCCGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGAGCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCg                                                                                                                                           >  3:159512/1‑149 (MQ=255)
       gTGCCGGTCGCGAGGATGTTGCCTGCCGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCgg                                                                                                                                          <  2:195954/149‑1 (MQ=255)
        tGCCGGTCGCGAGGATGTTGCCTGCCGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGt                                                                                                                                         <  7:46125/149‑1 (MQ=255)
        tGCCGGTCGCGAGGATGTTGCCTGCCGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGt                                                                                                                                         <  5:220251/149‑1 (MQ=255)
          ccGGTCGCGAGGATGTTGCCTGCCGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTg                                                                                                                                        <  6:75653/148‑1 (MQ=255)
          ccGGTCGCGAGGATGTTGCCTGCCGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTg                                                                                                                                        <  6:328469/148‑1 (MQ=255)
           cGGTCGCGAGGATGTTGCCTGCCGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCg                                                                                                                                      <  3:173424/149‑1 (MQ=255)
            ggTCGCGAGGATGTTGCCTGCCGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCgg                                                                                                                                     <  2:368606/149‑1 (MQ=255)
              tCGCGAGGATGTTGCCTGCCGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGt                                                                                                                                    >  3:355444/1‑148 (MQ=255)
                  gAGGATGTTGCCTGCCGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTg                                                                                                                               >  3:74533/1‑149 (MQ=255)
                  gAGGATGTTGCCTGCCGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTg                                                                                                                               >  7:248127/1‑149 (MQ=255)
                     gATGTTGCCTGCCGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTgcgc                                                                                                                            >  8:39282/1‑149 (MQ=255)
                      aTGTTGCCTGCCGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCgct                                                                                                                           <  4:355444/149‑1 (MQ=255)
                         ttGCCTGCCGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCgctgct                                                                                                                        >  2:318632/1‑149 (MQ=255)
                           gcctgccGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTg                                                                                                                                        <  2:254570/131‑1 (MQ=255)
                           gcctgccGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTg                                                                                                                      >  4:271467/1‑149 (MQ=255)
                           gcctgccGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGAACACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTg                                                                                                                                        >  1:254570/1‑131 (MQ=255)
                            cctgccGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTgg                                                                                                                     <  4:113021/149‑1 (MQ=255)
                            cctgccGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTgg                                                                                                                     >  4:244696/1‑149 (MQ=255)
                               gccGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGACTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGa                                                                                                                  >  4:271453/1‑149 (MQ=255)
                                ccGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGa                                                                                                                  <  3:309450/148‑1 (MQ=255)
                                ccGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGAt                                                                                                                 <  7:384796/149‑1 (MQ=255)
                                   cgTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGAtgtg                                                                                                              <  5:73132/149‑1 (MQ=255)
                                    gTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATgtgt                                                                                                             <  5:267470/149‑1 (MQ=255)
                                     tCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGtt                                                                                                            <  1:318632/149‑1 (MQ=255)
                                       aCGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCAAGGAATGGCGAGGTGATAGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTg                                                                                                           >  1:131438/1‑149 (MQ=25)
                                         gTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGc                                                                                                          >  8:182538/1‑147 (MQ=255)
                                         gTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGc                                                                                                          >  5:312538/1‑147 (MQ=255)
                                          tACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAg                                                                                                       <  7:324186/149‑1 (MQ=255)
                                               gCACGGTAGAGCCGGCTACGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCa                                                                                                  >  2:372429/1‑149 (MQ=255)
                                                 aCGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATc                                                                                                >  6:355388/1‑149 (MQ=255)
                                                   ggTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCgg                                                                                              >  4:359781/1‑149 (MQ=255)
                                                   ggTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCgg                                                                                              >  3:98530/1‑149 (MQ=255)
                                                   ggTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCgg                                                                                              >  6:281365/1‑149 (MQ=255)
                                                       gagCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGc                                                                                          <  3:271453/149‑1 (MQ=255)
                                                       gagCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGc                                                                                          >  2:383778/1‑149 (MQ=255)
                                                         gCCGGCTTTGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCg                                                                                        >  8:225817/1‑149 (MQ=255)
                                                         gCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCg                                                                                        >  5:307472/1‑149 (MQ=255)
                                                              cTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCt                                                                                   >  2:106770/1‑149 (MQ=255)
                                                              cTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAATCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCt                                                                                   <  7:263228/149‑1 (MQ=255)
                                                                  gCCTCGACCCCTGACTTGGCTGCCGGCGAGG‑CCGATCACCAAGTTTGTCTTGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCTCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCTGCGTGTCTTgca                                                                               <  8:148301/149‑1 (MQ=255)
                                                                  gCCGCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCTGTCGGGGCTTCCGGTGGCGTGGT‑ATCGACCCGGTCGATATCGGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGAGGCATCGGTCGCCGTGGCTTgc                                                                                <  7:61772/148‑1 (MQ=38)
                                                                   ccTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTgcag                                                                              <  8:248127/149‑1 (MQ=255)
                                                                            cTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGcc                                                                     >  1:370631/1‑149 (MQ=255)
                                                                               aGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGCCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCg                                                                    <  4:159512/147‑1 (MQ=255)
                                                                                 ttGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCaa                                                                >  1:347661/1‑149 (MQ=255)
                                                                                   ggCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTgg                                                                                     <  2:357862/126‑1 (MQ=255)
                                                                                   ggCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTgg                                                                                     >  1:357862/1‑126 (MQ=255)
                                                                                   ggCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAAc                                                               <  7:39282/148‑1 (MQ=255)
                                                                                   ggCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAAc                                                               <  3:271467/148‑1 (MQ=255)
                                                                                         cGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTtgt                                                           >  3:162769/1‑146 (MQ=255)
                                                                                         cGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTtgt                                                           <  4:162769/146‑1 (MQ=255)
                                                                                          ggCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCg                                                                    <  2:60104/136‑1 (MQ=255)
                                                                                           gCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCgg                                                      <  7:182538/149‑1 (MQ=255)
                                                                                            cGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCggg                                                     >  2:316917/1‑149 (MQ=255)
                                                                                                g‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCGggtt                                                   >  8:231377/1‑147 (MQ=255)
                                                                                                  cctaTCACCAAGTCAGTCGGGGCTGCGGGTGGCGTGTT‑ATCTACCCCGTCTATATCCGGTGCGTTAACTGCGCTGCTTGGCGATGTGTTGCCTGCGGCATCGGTCGCCGTGGTTTGCAACACTTCGCCGTCAACTTGTGGCGggttgag                                                <  1:242641/146‑1 (MQ=25)
                                                                                                  ccGATCACGAAGTCAGTCGGGGCTGCGGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCGggttgag                                                <  7:225817/149‑1 (MQ=255)
                                                                                                  ccGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCGggttgag                                                <  1:106770/149‑1 (MQ=255)
                                                                                                  ccGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCGggttgag                                                <  4:98530/149‑1 (MQ=255)
                                                                                                  ccGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCGggttgag                                                <  2:131438/149‑1 (MQ=255)
                                                                                                   cGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCATGGCTTgcag                                                                              >  2:35289/1‑118 (MQ=255)
                                                                                                   cGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTAGCCAGCGGCATCGGTCGCCATGGCTTgcag                                                                              <  1:35289/118‑1 (MQ=255)
                                                                                                      tCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCGggttgaggttg                                            <  6:307472/149‑1 (MQ=255)
                                                                                                        aCCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCGggttgaggttga                                           <  1:316917/148‑1 (MQ=255)
                                                                                                        aCCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCGGgttgaggttgag                                          <  6:312538/149‑1 (MQ=255)
                                                                                                          cAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCGGGTTGAGGTTGAGGc                                        <  4:74533/149‑1 (MQ=255)
                                                                                                             gtcagtcGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCGggttga                                                 >  4:329615/1‑137 (MQ=255)
                                                                                                             gtcagtcGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCGggttga                                                 <  3:329615/137‑1 (MQ=255)
                                                                                                              tcagtcGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCGGGTTGAGGTTGAGGCTGaa                                    <  2:347661/149‑1 (MQ=255)
                                                                                                                agtcGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCGGGTTGAGGTTGAGGCTGAac                                   <  2:370631/148‑1 (MQ=255)
                                                                                                                  tcGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCGGGTTGAGGTTGAGGCTGAACACg                                >  3:173715/1‑149 (MQ=255)
                                                                                                                    ggggCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCGGGTTGAGGTTGAGGCTGAACACGc                               <  1:372429/148‑1 (MQ=255)
                                                                                                                    ggggCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCGGGTTGAGGTTGAGGCTGAACACGc                               <  3:359781/148‑1 (MQ=255)
                                                                                                                    ggggCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCGGGTTGAGGTTGAGGCTGAACACGc                               <  5:355388/148‑1 (MQ=255)
                                                                                                                      ggCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCGGGTTGAGGTTGAGGCTGAACACGCCgg                            <  3:244696/149‑1 (MQ=255)
                                                                                                                       gCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCGGGTTGAGGTTGAGGCTGAACACGCCGGt                           <  4:173715/149‑1 (MQ=255)
                                                                                                                         tGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCGGGTTGAGGTTGAGGCTGAACACGCCGGTCg                         >  2:118978/1‑149 (MQ=255)
                                                                                                                                 gCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCGGGTTGAGGTTGAGGCTGAACACGCCGGTCGGGCCGAc                  <  5:281365/148‑1 (MQ=255)
                                                                                                                                 accggat‑ATCGACGGGGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCGGGTTGAGGTTGAGGCTGAACACGCCGGTCGGGCCGAc                  >  5:248614/7‑148 (MQ=18)
                                                                                                                                  ccggat‑ATCGACGGGGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCGGGTTGAGGTTGAGGCTGAACACGCCGGTCGGGCCGACCa                >  6:301928/6‑149 (MQ=255)
                                                                                                                                    tGGT‑ATCGACCCCGTCGATTTCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTGAACTTGTGGAGGGTTGAGGTTGAGGCTGAACACGCCGGTCGGGCCGACCACg              >  7:219620/1‑149 (MQ=255)
                                                                                                                                    tGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCGGGTTGAGGTTGAGGCTGAACACGCCGGTCGGGCCGACCACg              >  6:78182/1‑149 (MQ=255)
                                                                                                                                     gat‑ATCGACGGGGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCGGGttgaggttgagg                                         <  7:346020/119‑1 (MQ=255)
                                                                                                                                     gat‑ATCGACGGGGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCGGGttgaggttgagg                                         >  8:346020/3‑121 (MQ=255)
                                                                                                                                       t‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCGGGTTGAGGTTGAGGCTGAACACGCCGGTCGGGCCGACCACGGCa           <  7:231377/149‑1 (MQ=255)
                                                                                                                                            gACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCGGGTTGAGGTTGAGGCTGAACACGCCGGTCGGGCCGACCACGGCAGTGc       <  2:153073/149‑1 (MQ=255)
                                                                                                                                            gACCCCGTCGATATCCGGGGCGGCAACGGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCGGGTTGAGGTTGAGGCTGAACACGCCGGTCGGGCCGACCACGGCAGTg        <  8:271517/148‑1 (MQ=255)
                                                                                                                                             aCCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCGGGTTGAGGTTGAGGCTGAACACGCCGGTCGGGCCGACCACGGCAGTGcc      >  6:250761/1‑149 (MQ=255)
                                                                                                                                                 cGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCGGGTTGAGGTTGAGGCTGAACACGCCGGTCGGGCCGACCACGGCAGTGCCCAGt  >  8:261102/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GGCGACGGTGCCGGTCGCGAGGATGTTGCCTGCCGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGG‑CCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGT‑ATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCGGGTTGAGGTTGAGGCTGAACACGCCGGTCGGGCCGACCACGGCAGTGCCCAGT  >  NC_002947/943468‑943759

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