Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 133,296 C→G 5.8% P25R (CCG→CGG)  PP_0128 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947133,2960CG 5.8% 280.9 / 6.0 86P25R (CCG→CGG) PP_0128hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (39/42);  new base G (2/3);  total (41/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.25e-01

CCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGC‑CAAGTAGGCGGCGTCACCC‑ATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGG  >  NC_002947/133151‑133443
                                                                                                                                                      |                                                                                                                                                   
ccGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGC‑CAAGTAGGCGGCGTCACCC‑ATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCggg                                                                                                                                                  >  2:319953/1‑149 (MQ=255)
 cGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGC‑CAAGTAGGCGGCGTCACCC‑ATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCGGGGc                                                                                                                                                 <  1:158453/149‑1 (MQ=255)
 cGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGC‑CAAGTAGGCGGCGTCACCC‑ATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGc                                                                                                                                                 <  2:184527/149‑1 (MQ=255)
 cGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGC‑CAAGTAGGCGGCGTCACCC‑ATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGc                                                                                                                                                 >  4:373791/1‑149 (MQ=255)
   gAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGC‑CAAGTAGGCGGCGTCACCC‑ATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCt                                                                                                                                               >  6:269247/1‑149 (MQ=255)
      ggCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGC‑CAAGTAGGCGGCGTCACCC‑ATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGcc                                                                                                                                            <  2:58516/149‑1 (MQ=255)
       gCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGC‑CAAGTAGGCGGCGTCACCC‑ATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCGGGGCCTGcag                                                                                                                                           <  1:145993/149‑3 (MQ=255)
         gTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGC‑CAAGTAGGCGGCGTCACCC‑ATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGccc                                                                                                                                           >  8:229533/1‑147 (MQ=255)
              gAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGC‑CAAGTAGGCGGCGTCACCC‑ATGCCCCGCGGATTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAg                                                                                                                                        >  1:330328/1‑148 (MQ=255)
               aaTGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGC‑CAAGTAGGCGGCGTCACCC‑ATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGAt                                                                                                                                     <  5:274320/147‑1 (MQ=255)
               aaTGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGC‑CAAGTAGGCGGCGTCACCC‑ATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATg                                                                                                                                    <  5:18008/148‑1 (MQ=255)
                     gCCGACCAGCTGATCGACAAGGAGCGCGCTGC‑CAAGTAGGCGGCGTCACCC‑ATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTct                                                                                                                             <  1:357362‑M2/149‑3 (MQ=255)
                      ccGACCAGCTGATCGACAAGGAGCGCGCTGCGCAAGTAGGCGGCGTCACCC‑ATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTct                                                                                                                             <  1:42349‑M2/149‑3 (MQ=255)
                          ccAGCTGATCGACAAGTAGCGCGCTGC‑CAAGTAGGCGGCGTCACCC‑ATGCCCCGC‑‑‑TTCAGTGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTctggggc                                                                                                                        >  4:115126‑M2/1‑142 (MQ=255)
                          ccAGCTGATCGACAAGGAGCGCGCTGC‑CAAGTAGGCGGCGTCACCC‑ATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCGGGGCCTGcaggtga                                                                                                                                       <  7:226453/134‑7 (MQ=255)
                          ccAGCTGATCGACAAGGAGCGCGCTGC‑CAAGTAGGCGGCGTCACCC‑ATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCGGGGCCTGcaggtga                                                                                                                                       >  8:226453/1‑128 (MQ=255)
                           cAGCTGATCGACAAGGAGCGCGCTGC‑CAAGTAGGCGGCGTCACCC‑ATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTctgcggct                                                                                                                       <  5:38788‑M2/149‑9 (MQ=255)
                            aGCTGATCGACAAGGAGCGCGCTGC‑CAAGTAGGCGGCGTCACCC‑ATGCCCAGC‑‑‑TTGAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACGTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTctgcgtctg                                                                                                                      >  4:353011‑M2/1‑140 (MQ=255)
                              cTGATCGACAAGGAGCGCGCTGC‑CAAGTAGGCGGCGTCACCC‑ATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGgctgct                                                                                                                    >  5:326612/1‑149 (MQ=255)
                                gATCGACAAGGAGCGCGCTGC‑CAAGTAGGCGGCGTCACCC‑ATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGgctgct                                                                                                                    >  6:76878/1‑147 (MQ=255)
                                 aTCGACAAGGAGCGCGCTGC‑CAAGTAGGCGGCGTCACCC‑ATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAg                                                                                                                 <  3:109816/149‑1 (MQ=255)
                                 aTCGACAAGGAGCGCGCTGC‑CAAGTAGGCGGCGTAACCC‑ATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAg                                                                                                                 >  4:258505/1‑149 (MQ=255)
                                    gACAAGGAGCGCGCTGC‑CAAGTAGGCGGCGTCACCC‑ATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGttt                                                                                                              <  3:373791/149‑1 (MQ=255)
                                    gACAAGGAGCGCGCTGC‑CAAGTAGGCGGCGTCACCC‑ATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGttt                                                                                                              <  6:121168/149‑1 (MQ=255)
                                     aCAAGGAGCGCGCTGC‑CAAGTAGGCGGCGTCACCC‑ATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGttt                                                                                                              >  3:143189/1‑148 (MQ=255)
                                        atgAGCGCGCTTC‑CAAGTAGTCGGCGTCAACC‑ATTCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTTGCCTTCACGAGCCGCATGTCAATGAGCATCACCTTTAGGCCCCGGTCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAAc                                                                                                          <  7:255332/147‑1 (MQ=255)
                                        aGGAGCGCGCTGC‑CAAGTAGGCGGCGTCACCCAATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCac                                                                                                           >  1:48368/1‑148 (MQ=255)
                                        aGGAGCGCGCTGC‑CAAGTAGGCGGCGTCACCC‑ATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAAc                                                                                                          <  7:267953/149‑1 (MQ=255)
                                            gcgcgcTGC‑CAAGTAGGCGGCGTCACCC‑ATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCCGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATcg                                                                                                      <  6:280883/149‑2 (MQ=255)
                                                 cTGC‑CAAGTAGGCGGCGTCACCC‑ATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTg                                                                                                 <  2:355635/149‑1 (MQ=255)
                                                  tGC‑CAAGTAGGCGGCGTCACCC‑ATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTgg                                                                                                >  2:379175/1‑149 (MQ=255)
                                                    c‑cAAGTAGGCGGCGTCACCC‑ATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAAc                                                                                                          <  5:7154/137‑1 (MQ=255)
                                                    c‑cAAGTAGGCGGCGTCACCC‑ATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAAc                                                                                                          >  6:7154/1‑137 (MQ=255)
                                                    c‑cAAGTAGGCGGCGTCACCC‑ATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTgg                                                                                                <  5:82397/147‑1 (MQ=255)
                                                       aaGTAGGCGGCGTCACCC‑ATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCAt                                                                                            >  5:220862/1‑149 (MQ=255)
                                                       aaGTAGGCGGCGTCACCC‑ATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGAAGGCCGTCTGCGGCTGCTCAGTTTCAACATGCAGGTGGGCAt                                                                                            >  5:349654/1‑149 (MQ=255)
                                                         gTAGGCGGCGTCACCC‑ATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGTGCATCa                                                                                          <  8:83377/149‑1 (MQ=255)
                                                         gTAGGCGGCGTCACCC‑ATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCa                                                                                          <  2:357362/149‑1 (MQ=255)
                                                                  gTCACCC‑ATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAgc                                                                                 <  6:326612/149‑1 (MQ=255)
                                                                       cc‑ATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTAc                                                                            <  7:229533/149‑1 (MQ=255)
                                                                       cc‑ATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTAc                                                                            >  1:140571/1‑149 (MQ=255)
                                                                       cc‑ATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTAc                                                                            >  4:119273/1‑149 (MQ=255)
                                                                             ccccGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGcc                                                                        <  2:42349/148‑1 (MQ=255)
                                                                              cccGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCAt                                                                      >  2:291505/1‑149 (MQ=255)
                                                                                      ttCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCt                                                                 <  2:330328/149‑1 (MQ=255)
                                                                                         agagCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGAc                                                              <  8:241016/149‑1 (MQ=255)
                                                                                            gCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGAcgcg                                                           >  4:258582/1‑149 (MQ=255)
                                                                                                  cggggCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTg                                                     <  4:227591/146‑1 (MQ=255)
                                                                                                   gTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTgg                                                    <  3:258582/149‑1 (MQ=255)
                                                                                                      gCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGgca                                                  <  6:38788/148‑1 (MQ=255)
                                                                                                           ggCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACct                                            >  7:174525/1‑149 (MQ=255)
                                                                                                           ggCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACct                                            >  2:73905/1‑149 (MQ=255)
                                                                                                            gCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAg                                                        <  2:257265/136‑1 (MQ=255)
                                                                                                            gCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAg                                                        >  1:257265/1‑136 (MQ=255)
                                                                                                              cTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACctgct                                         >  8:336651/1‑149 (MQ=255)
                                                                                                               tGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCGACCGCCATTACCTGACGCGTAGCTGGCAGCACctgctg                                        >  4:223606/1‑149 (MQ=255)
                                                                                                                 cACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGcc                                      >  2:298637/1‑149 (MQ=255)
                                                                                                                 cACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGcc                                      <  1:113144/149‑1 (MQ=255)
                                                                                                                  aCGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCg                                     >  8:76553/1‑149 (MQ=255)
                                                                                                                    gAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTcatcggctcgtgcaggccaatgacacgcgtacctctgaagcggggcaagggtgccgccgcctacttggcagcccgctccttgtcgatcagct                                   <  2:343949‑M2/149‑93 (MQ=255)
                                                                                                                    gAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCa                                   <  5:269247/149‑1 (MQ=255)
                                                                                                                      gCCGCAGGTCGACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTctgcggctcgtgcaggccaatgccacgcgtagctctgaagcggggcatgggtgacgccgcctacttggcagcgcgctccttgtcgatcagctgg                                 <  3:346396‑M2/149‑95 (MQ=255)
                                                                                                                      gCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCAt                                  >  6:87995/1‑148 (MQ=255)
                                                                                                                      gCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCAt                                  >  7:333036/1‑148 (MQ=255)
                                                                                                                       ccGCAGGTCAACGAGCATCACCTGCAGGCCGCGGGCCTGCCCGAGGATGGCCGTctgcggctcgtgcaggccaatgccacgcgtagctctgaagcggggcatgggtgacgccgcctacttggcagcgcgctccttgtcgatcagctggt                                <  1:33096‑M2/149‑96 (MQ=255)
                                                                                                                         gCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCAc                                              <  1:163392/133‑1 (MQ=255)
                                                                                                                         gCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCAc                                              >  2:163392/1‑133 (MQ=255)
                                                                                                                         gCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACg                              <  8:174525/149‑1 (MQ=255)
                                                                                                                          cAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACgg                             <  8:333036/149‑1 (MQ=255)
                                                                                                                           aGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACggg                            <  3:258505/149‑1 (MQ=255)
                                                                                                                               gaaCGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTg                        >  7:70592/2‑149 (MQ=255)
                                                                                                                               cAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTg                        >  6:183379/1‑149 (MQ=255)
                                                                                                                                aaCGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGc                       >  2:5400/1‑149 (MQ=255)
                                                                                                                                aaCGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGAGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGc                       >  6:128935/1‑149 (MQ=255)
                                                                                                                                  cGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCGCCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCg                     >  4:262655/1‑149 (MQ=255)
                                                                                                                                  cGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCg                     <  6:21036/149‑1 (MQ=255)
                                                                                                                                  cGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCg                     <  4:143189/149‑1 (MQ=255)
                                                                                                                                   gAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGGCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCgg                    <  1:73905/149‑1 (MQ=255)
                                                                                                                                      cATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTaa                 >  1:393126/1‑149 (MQ=255)
                                                                                                                                           ccTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTg             >  8:55054/1‑148 (MQ=255)
                                                                                                                                           ccTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGc            <  7:55054/149‑1 (MQ=255)
                                                                                                                                              gCAGGCCCGGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCag          >  3:230098/1‑148 (MQ=255)
                                                                                                                                                 ggCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGt                         <  1:33295/130‑1 (MQ=255)
                                                                                                                                                 ggCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGt                         >  2:33295/1‑130 (MQ=255)
                                                                                                                                                 ggCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCagaaga      <  1:160431/149‑1 (MQ=255)
                                                                                                                                                  gCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCag          >  5:14571/1‑144 (MQ=255)
                                                                                                                                                  cccccGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCag          >  7:275679/2‑144 (MQ=255)
                                                                                                                                                   ccccGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATc    <  1:291505/149‑1 (MQ=255)
                                                                                                                                                    cccGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCag          >  2:23793/1‑142 (MQ=255)
                                                                                                                                                    cccGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCg   <  5:76878/149‑1 (MQ=255)
                                                                                                                                                    cccGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAACTGGCAGTACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCag          >  5:238853/1‑142 (MQ=255)
                                                                                                                                                      cGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCgg  >  5:58577/1‑148 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
CCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGC‑CAAGTAGGCGGCGTCACCC‑ATGCCCCGC‑‑‑TTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGG  >  NC_002947/133151‑133443

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