Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 223,870 G→C 7.5% G233G (GGG→GGC) ‡ PP_0171 → ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947223,8700GC 7.5% 255.5 / 9.2 80G233G (GGG→GGC) ‡PP_0171ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (36/38);  new base C (4/2);  total (40/40)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.09e-01

CCCTGGATGCGCTGACCCGGCGCAAGATGCAGGAAGAGCTGTTGCTGCTGTGGGAGGAAGTGCGGTTCACCCTGTTGTTCGTCACCCACTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTCCAGGCCAGTGCCCGGCGTATCCACCGTTTGCTGTTCGATGAAGCAGATATCCCGGAACAGCCTGACAACCTTGGGTTCAAC  >  NC_002947/223722‑224014
                                                                                                                                                    |                                                                                                                                                           
cccTGGATGCGCTGACCCGGCGCAAGATGCAGGAAGAGCTGTTGCTGCTGTGGGAGGAAGTGCGGTTCACCCTGTTGTTCGTCACCCACTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCgggg                                                                                                                                                             >  6:158708/1‑149 (MQ=255)
cccTGGATGCGCTGACCCGGCGCAAGATGCAGGAAGAGCTGTTGCTGCTGTGGGAGGAAGTGCGGTTCACCCTGTTGTTCGTCACCCACTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCgggg                                                                                                                                                             >  1:334286/1‑149 (MQ=255)
     gATGCGCTGACCCGGCGCAAGATGCAGGAAGAGCTGTTGCTGCTGTGGGAGGAAGTGCGGTTCACCCTGTTGTTCGTCACCCACTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCgg                                                                                                                                                          <  3:44947/147‑1 (MQ=255)
     gATGCGCTGACCCGGCGCAAGATGCAGGAAGAGCTGTTGCTGCTGTGGGAGGAAGTGCGGTTCACCCTGTTGTTCGTCACCCACTCCATCGAAGAAGCACTGGTGGGTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGt                                                                                                                                                        >  8:81913/1‑149 (MQ=255)
     gATGCGCGGACCCGGCGCAAGATGCAGGAAGAGCTGTTGCTGCTGTGGGAGGAAGTGCGGTTCACCCTGTTGTTCGTCACCCACTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCgg                                                                                                                                                          <  2:34147/147‑1 (MQ=255)
       tGCGCTGACCCGGCGCAAGATGCAGGAAGAGCTGTTGCTGCTGTGGGAGGAAGTGCGGTTCACCCTGTTGTTCGTCACCCACTCCATCGAAGAAGCACAGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTgc                                                                                                                                                      <  8:178618/149‑1 (MQ=255)
         cgcTGACCCGGCGCAAGATGCAGGAAGAGCTGTTGCTGCTGTGGGAGGAAGTGCGGTTCACCCTGTTGTTCGTCACCCACTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTgcgc                                                                                                                                                    <  2:59776/149‑1 (MQ=255)
         cgcTGACCCGGCGCAAGATGCAGGAAGAGCTGTTGCTGCTGTGGGAGGAAGTGCGGTTCACCCTGTTGTTCGTCACCCACTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGCCCCGGGTgcgg                                                                                                                                                    <  5:60114/149‑2 (MQ=255)
         cgcTGACCCGGCGCAAGATGCAGGAAGAGCTGTTGCTGCTGTGGGAGGAAGTGCGGTTCACCCTGTTGTTCGTCACCCACTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGCCCCGGGTgcgg                                                                                                                                                    <  4:121845/149‑2 (MQ=255)
          gcTGACCCGGCGCAAGATGCAGGAAGAGCTGTTGCTGCTGTGGGAGGAAGTGCGGTTCACCCTGTTGTTCGTCACCCACTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCg                                                                                                                                                           <  1:282755/141‑1 (MQ=255)
          gcTGACCCGGCGCAAGATGCAGGAAGAGCTGTTGCTGCTGTGGGAGGAAGTGCGGTTCACCCTGTTGTTCGTCACCCACTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCg                                                                                                                                                           >  2:282755/1‑141 (MQ=255)
            tGACCCGGCGCAAGATGCAGGAAGAGCTGTTGCTGGTGTGGGAGGAAGTGGGGTTCACCCTGTTGTTCGTCACCCACTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGcc                                                                                                                                                 >  4:298619/1‑149 (MQ=255)
            tGACCCGGCGCAAGATGCAGGAAGAGCTGTTGCTGCTGTGGGAGGAAGTGCGGTTCACCCTGTTGTTCGTCACCCACTCCATCAAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGcc                                                                                                                                                 >  3:42638/1‑149 (MQ=255)
             gACCCGGCGCAAGATGCAGGAAGAGCTGTTGCTGCTGTGGGAGGAAGTGCGGTTCACCCTGTTGTTCGTCACCCACTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGcc                                                                                                                                                 >  6:348030/1‑148 (MQ=255)
               cccGGCGCAAGATGCAGGAAGAGCTGTTGCTGCTGTGGGCGGAAGTGCGGTGCACCCTGTTGTTCGTCACCCACTCCATCGAAGAAGCACTGGTGGTGGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAg                                                                                                                                              <  4:311463/149‑1 (MQ=255)
                ccGGCGCAAGATGCAGGAAGAGCTGTTGCTGCTGTGGGAGGAAGTGCGGTTCACCCTGTTGTTCGTCACCCACTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAgg                                                                                                                                             >  3:205823/1‑149 (MQ=255)
                 cGGCGCAAGATGCAGGAAGAGCTGTTGCTGCTGTGGGAGGAAGTGCGGTTCACCCTGTTGTTCGTCACCCACTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGt                                                                                                                                            <  2:269937/149‑1 (MQ=255)
                 cGGCGCAAGATGCAGGAAGAGCTGTTGCTGCTGTGGGAGGAAGTGCGGTTCACCCTGTTGTTCGTCACCCACTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGt                                                                                                                                            >  7:18833/1‑149 (MQ=255)
                 cGGCGCAAGATGCAGGAAGAGCTGTTGCTGCTGTGGGAGGAAGTGCGGTTCACCCTGTTGTTCGTCACCCACTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGt                                                                                                                                            >  4:166487/1‑149 (MQ=255)
                   gcgcAAGATGCAGGAAGAGCTGTTGCTGCTGTGGGAGGAAGTGCGGTTCACCCTGTTGTTCGTCACCCACTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGc                                                                                                                                          <  6:349270/149‑1 (MQ=255)
                       aaGATGCAGGAAGAGCTGTTGCTGCTGTGGGAGGAAGTGCGGTTCACCCTGTTGTTCGTCACCCACTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAg                                                                                                                                      >  7:121024/1‑149 (MQ=255)
                         gATGCAGGAAGAGCTGTTGCTGCTGTGGGAGGAAGTGCGGTTCACCCTGTTGTTCGTCACCCACTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGcc                                                                                                                                    <  3:298619/149‑1 (MQ=255)
                         gATGCAGGAAGAGCTGTTGCGGCTGTGGGAGGAAGTGCGGTTCACCCTGTTGTTCGTCACCCACTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGcc                                                                                                                                    <  4:109840/149‑1 (MQ=255)
                             cAGGAAGAGCTGTTGCTGCTGTGGGAGGAAGTGCGGTTCACCCTGTTGTTCGTCACCCACTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATc                                                                                                                                 <  2:121347/148‑1 (MQ=255)
                               ggAAGAGCTGTTGCTGCTGTGGGAGGAAGTGCGGTTCACCCTGTTGTTCGTCACCCACTCCATCAAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAAt                                                                                                                              <  4:42638/149‑1 (MQ=255)
                                 aaGAGCTGTTGCTGCTGTGGGAGGAAGTGCGGTTCACCCTGTTGTTCGTCACCCACTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAc                                                                                                                            <  3:67841/149‑1 (MQ=255)
                                     gCTGTTGCTGCTGTGGGAGGAAGTGCGGTTCACCCTGTTGTTCGTCACCCACTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCC                                                                                                              <  1:39369‑M1/149‑1 (MQ=255)
                                         ttGCTGCTGTGGGAGGAAGTGCGGTTCACCCTGTTGTTCGTCACCCACTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTggg                                                                                                          <  5:158708/149‑1 (MQ=255)
                                         ttGCTGCTGTGGGAGGAAGTGCGGTTCACCCTGTTGTTCGTCACCCACTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTG                                                                                                            <  4:119470‑M1/147‑1 (MQ=255)
                                             tgctgTGGGAGGAAGTGCGGTTCACCCTGTTGTTCGTCACCCACTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGc                                                                                                      >  8:32818/1‑149 (MQ=255)
                                                    ggaggaAGTGCGGTTCACCCTGTTGTTCGTCACCCACTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCaa                                                                                                                               >  2:178752/1‑127 (MQ=255)
                                                    ggaggaAGTGCGGTTCACCCTGTTGTTCGTCACCCACTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCaa                                                                                                                               <  1:178752/127‑1 (MQ=255)
                                                     gaggaAGTGCGGTTCACCCTGTTGTTCGTCACCCACTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTgg                                                                                              <  7:32818/149‑1 (MQ=255)
                                                           gTGCGGTTCACCCTGTTGTTCGTCACCCACTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGa                                                                                        <  5:348030/149‑1 (MQ=255)
                                                                 ttCACCCTGTTGTTCGTCACCCACTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTcca                                                                                  >  1:77881‑M2/1‑146 (MQ=255)
                                                                   cACCCTGTTGTTCGTCACCCACTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCAATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTCGGTTGAAGCGACTTccagg                                                                                <  2:285909‑M2/149‑6 (MQ=255)
                                                                      ccTGTTGTTCGTCACCCACTCCATCGAAGAAGCACTGGTGGTTGGCAACCTTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGaagtcgctgccacccaggctgcccaggccg                                                                             >  3:294853‑M1/1‑119 (MQ=255)
                                                                      ccTGTTGTTCGTCACCCACTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTCCAGGCCa                                                                             <  8:121024/149‑1 (MQ=255)
                                                                      ccTGTTGTTCGTCACCCACTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAAAAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGACTGGGTGGCAGCGACTTccaggc                                                                               <  4:256098‑M2/147‑7 (MQ=255)
                                                                      ccTGTTGTTCGTCACCCACTCCATCGAAGAAGCACTGGTGGGTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGGGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTCCAGGCCa                                                                             <  1:372977/149‑1 (MQ=255)
                                                                          ttgttCGTCACCCACTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGaagtcgctgccacccaggctgcccaggccgta                                                                           <  8:18833‑M1/147‑33 (MQ=255)
                                                                             ttCGTCACCCACTCCATCGAAGAAGCACTGGTGGTTGGCAACCTTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGaagtcgctgccacccaggctgcccaggccgtattgat                                                                      <  4:294853‑M1/149‑38 (MQ=255)
                                                                                   aCCCACTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTCCAGGCCAGTGCCCGGCGTAt                                                               >  7:81368/1‑149 (MQ=255)
                                                                                    cccACTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTCCAGGCCAGTGCCCGGCGTATc                                                              >  1:60034/1‑149 (MQ=255)
                                                                                       aCTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATACGCCGGGCACTGGCCTGGaagtcgctgccacccaggctgcccaggccgtattgat                                                                      <  2:84298‑M1/149‑38 (MQ=255)
                                                                                       aCTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTccaggc                                                                               <  7:261897‑M2/130‑7 (MQ=255)
                                                                                       aCTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTccaggc                                                                               >  8:261897‑M2/1‑124 (MQ=255)
                                                                                           cATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTCCAGGCCAGTGCCCGGCGTATCCACCGtt                                                       <  3:166487/149‑1 (MQ=255)
                                                                                             tCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTCCAGGCCAGTGCCCGGCGTATCCACCGTTTg                                                     >  7:307703/1‑149 (MQ=255)
                                                                                               gaagaaGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTCCAGGCCAGTGCCCGGCGTATCCACCGTTTGCt                                                   <  2:334286/149‑1 (MQ=255)
                                                                                               gaagaaGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTCCAGGCCAGTGCCCGGCGTATCCACCGTTTGCt                                                   >  8:254807/1‑149 (MQ=255)
                                                                                                aagaagCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTCCAGGCCAGTGCCCGGCGTATCCACCGTTTGCTg                                                  >  2:273662/1‑149 (MQ=255)
                                                                                                  gaagCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTCCAGGCCAGTGCCCGGCGTATCCACCGTTTGCTGtt                                                <  8:307703/149‑1 (MQ=255)
                                                                                                      cACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATTGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTCCAGGCCAGTGCCCGGCGTATCCACCGTTTGCTGTTCGAt                                            >  7:253989/1‑149 (MQ=255)
                                                                                                      cACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTCCAGGCCAGTGCCCGGCGTATCCACCGTTTGCTGTTCGAt                                            <  8:104588/149‑1 (MQ=255)
                                                                                                             ggtTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTCCAGGCCAGCTGCCCGGCGTATCCACCGTTTGCt                                                   <  7:237990/136‑1 (MQ=255)
                                                                                                             ggtTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTCCAGGCCAGTGCCCGGCGTATCCACCGTTTGCt                                                   >  8:237990/1‑135 (MQ=255)
                                                                                                             ggtTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTCCAGGCCAGTGCCCGGCGTATCCACCGTTTGCTGTTCGATGAAg                                        >  4:305194/1‑146 (MQ=255)
                                                                                                              gtTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTCCAGGCCAGTGCCCGGCGTATCCACCGTTTGCTGTTCGATGAAGCAg                                     <  3:305194/148‑1 (MQ=255)
                                                                                                                     aCCGTATCCTGCTGTTGTCGCCGCACCCGCGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTCCAGGCCAGTGCCCGGCGTATCCAGCGTTTGCTGTTCGATGAAGCAGATATCCCg                             >  4:235002/1‑149 (MQ=255)
                                                                                                                      ccGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTCCAGGCCAGTGCCCGGCGTATCCACCGTTTGCTGTTCGATGAAGCAGATATCCCgg                            <  8:81368/149‑1 (MQ=255)
                                                                                                                            ccTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTCCAGGCCAGTGCCCGGCGTATCCACCGTTTGCTGTTCGATGAAGCAGATATCCCGGAACAGc                      >  6:154229/1‑149 (MQ=255)
                                                                                                                             ctgctgTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTCCAGGCCAGTGCCCGGCGTATCCACCGTTTGCTGTTCGATGAAg                                        >  1:156563/1‑130 (MQ=255)
                                                                                                                             ctgctgTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTCCAGGCCAGTGCCCGGCGTATCCACCGTTTGCTGTTCGATGAAGCAGATATCCCGGAACAGcc                     <  4:205823/149‑1 (MQ=255)
                                                                                                                             ctgctgTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTCCAGGCCAGTGCCCGGCGTATCCACCGTTTGCTGTTCGATGAAGCAGATATCCCGGAACAGc                      <  7:81913/148‑1 (MQ=255)
                                                                                                                              tgctgTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTCCAGGCCAGTGCCCGGCGTATCCACCGTTTGCTGTTCGATGAAGCAGATATCCCGGAACAGCCt                    >  7:183559/1‑149 (MQ=255)
                                                                                                                              tgctgTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTCCAGGCCAGTGCCCGGCGTATCCACCGTTTGCTGTTCGATGAAGCAGATATCCCGGAACAGCCt                    >  2:104754/1‑149 (MQ=255)
                                                                                                                                ctgTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTCCAGGCCAGTGCCCGGCGTATCCACCGTTTGCTGTTCGATGAAGCAGATATCCCGGAACAGCCTGa                  >  2:35833/1‑149 (MQ=255)
                                                                                                                                   ttgtCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTCCAGGCCAGTGCCCGGCGTATCCACCGTTTGCTGTTCGATGAAGCAGATATCCCGGAACAGCCTGACaa               >  4:168109/1‑149 (MQ=255)
                                                                                                                                      tCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTCCAGGCCAGTGCCCGGCGTATCCACCGTTTGCTGTTCGATGAAGCAGATATCCCGGAACAGCCTGACAACCt            >  2:326888/1‑149 (MQ=255)
                                                                                                                                        gccgcACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTCCAGGCCAGTGCCCGGCGTATCCACCGTTTGCTGTTCGATGAAGCAGATATCCCGGAACAGCCTGACAACCTTg          >  4:15543/1‑149 (MQ=255)
                                                                                                                                        gccgcACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTG                                                                                                            >  4:140671‑M1/1‑52 (MQ=255)
                                                                                                                                        gccgcACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTG                                                                                                            <  3:140671‑M1/52‑1 (MQ=255)
                                                                                                                                        cgcgcACCCGCCCCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTCCAGGCCAGTGCCCGGCGTATCCACCGTTTGCTGTTCGATGAAGCAGATATCCCGGAACAGCCTGACAACCTTg          >  3:87519/3‑149 (MQ=255)
                                                                                                                                         gcgcACCCGCCCCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTCCAGGCCAGTGCCCGGCGTATCCACCGTTTGCTGTTCGATGAAGCAGATATCCCGGAACAGCCTGACAACCTTgg         >  5:326391/2‑149 (MQ=255)
                                                                                                                                         gcgcACCCGCCCCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTCCAGGCCAGTGCCCGGCGTATCCACCGTTTGCTGTTCGATGAAGCAGATATCCCGGAACAGCCTGACAACCTTgg         >  2:218893/2‑149 (MQ=255)
                                                                                                                                         ccgcACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTCCAGGCCAGTGCCCGGCGTATCCACCGTTTGCTGTTCGATGAAGCAGATATCCCGGAACAGCCTGACAACCTTgg         <  2:60034/149‑1 (MQ=255)
                                                                                                                                         ccgcACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTCCAGGCCAGTGCCCGGCGTATCCACCGTTTGCTGTTCGATGAAGCAGATATCCCGGAACAGCCTGACAACCTTgg         >  4:285721/1‑149 (MQ=255)
                                                                                                                                           gcACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTCCAGGCCAGTGCCCGGCGTATCCACCGTTTGCTGTTCGATGAAGCAGATATCCCGGAACAGCCTGACAACCTTGGGt       >  5:185420/1‑149 (MQ=255)
                                                                                                                                           gcACCCGCCCCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTCCAGGCCAGTGCCCGGCGTATCCACCGTTTGCTGTTCGATGAAGCAGATATCCCGGAACAGCCTGACAACCTTGGGt       >  1:137768/1‑149 (MQ=255)
                                                                                                                                              cccGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTCCAGGCCAGTGCCCGGCGTATCCACCGTTTGCTGTTCGATGAAGCAGATATCCCGGAACAGCCTGACAACCTTGGGTTc     >  1:308328/1‑148 (MQ=255)
                                                                                                                                                cGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTCCAGGCCAGTGCCCGGCGTATCCACCGTTTGCTGTTCGATGAAGCAGATATCCCGGAACAGCCTGACAACCTTGGGTTCAAc  >  4:214610/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                           
CCCTGGATGCGCTGACCCGGCGCAAGATGCAGGAAGAGCTGTTGCTGCTGTGGGAGGAAGTGCGGTTCACCCTGTTGTTCGTCACCCACTCCATCGAAGAAGCACTGGTGGTTGGCAACCGTATCCTGCTGTTGTCGCCGCACCCGGGGCGGGTGCGCGCCGAGGTGCATAGCCATCAATAC‑‑‑‑‑‑‑‑‑‑GGCCTGGGCAGCCTGGGTGGCAGCGACTTCCAGGCCAGTGCCCGGCGTATCCACCGTTTGCTGTTCGATGAAGCAGATATCCCGGAACAGCCTGACAACCTTGGGTTCAAC  >  NC_002947/223722‑224014

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: