Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 195,629 A→G 12.4% T379A (ACC→GCC)  PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947195,6290AG12.4% 190.5 / 18.0 92T379A (ACC→GCC) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base A (36/42);  new base G (6/5);  total (43/49)
Fisher's exact test for biased strand distribution p-value = 7.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

GCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCA  >  NC_002947/195483‑195777
                                                                                                                                                  |                                                                                                                                                    
gCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGTAATTTCGAGCAACTGCCCCCGAATACCACCCCGGCTCAGACCACGAACAACGACTCGGTCGacac                                                                                                                                                     >  3:123273/1‑148 (MQ=11)
     gATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGACTCCGACCCCGGCTCAGCCCACGATCAACCGCTCGGTCGAcaccaccac                                                                                                                                               >  8:246372/1‑149 (MQ=255)
     gATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGCCCCCGGCTCAGACCACGATCAACGACTCGGTCGAcaccaccac                                                                                                                                               >  2:240368/1‑149 (MQ=255)
            agtccgcGGACAATGACGTCTTCATCATCGGCTCGACCGTCAGCGTCACGATGGTGAGCGGCACGGGCGGTATTTTCGAGCAACTGTGCCCGAATCCGACCCCGGCTCAGACCACGTTGAACGAGTCGGTCGACACCACGACCGCGAcc                                                                                                                                        <  4:123273/143‑1 (MQ=2)
                  cGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACg                                                                                                                                  <  5:61030/149‑1 (MQ=39)
                    gcgaaTGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGAAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGc                                                                                                                                >  3:269202/4‑149 (MQ=0)
                    gCCAATGAAGTATACAACAACGGTTCGCCCGTCAGCGTCACCATCGATAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGc                                                                                                                                <  7:47438/149‑1 (MQ=255)
                      gaaTGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGa                                                                                                                              >  4:10758/2‑149 (MQ=14)
                        aTGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGTTCGACACCACTACCGCGACCCTGACGGCGAGc                                                                                                                            <  8:20752/149‑1 (MQ=12)
                            cGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGt                                                                                                                        <  2:6743/149‑1 (MQ=17)
                             gTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTc                                                                                                                       >  6:309730/1‑149 (MQ=17)
                                tACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGt                                                                                                                    >  1:233133/1‑149 (MQ=37)
                                 acaacaACGGTTCGACCGTCAGCGTCACCATCGAGAAACCCACAGGCGGTAAAATCGAGCAACAGACCCCGAATCCGACCCCGGCTCAGCCCACGATCAACGACTCCGTCGACACCACAACCGCGCCCCTGACGGCGAGCCCCTCggg                                                                                                                    >  8:92425/1‑147 (MQ=11)
                                  caacaaCGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCa                                                                                                                  <  7:252072/149‑1 (MQ=17)
                                   aacaacGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGCCCCCGGCGCAGACCACGATCACCGACTCGGTCGACCCCACCCCCGCGACCCTGACGGCGCGCCCGTCGGTCAc                                                                                                                 >  2:132352/1‑149 (MQ=11)
                                    acaacGGTTCGACCGTCAGCGTCACTATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCAcc                                                                                                                >  5:169730/1‑149 (MQ=12)
                                    acaacGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAACAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCAcc                                                                                                                >  6:133315/1‑149 (MQ=12)
                                     caacGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACGCCGAATCCGACCCCGGCTCAGACCACGCTCAACGACTCGGTCGACACCCCACCCGCGACCCTGACGGCGAGCCCGTCGGTccac                                                                                                                >  7:156686/1‑145 (MQ=11)
                                          gTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggc                                                                                                          <  4:159422/149‑1 (MQ=37)
                                              gACCGTCAGCTTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCGCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGgcggcg                                                                                                      <  8:267633/149‑1 (MQ=255)
                                              gACCGTCAGCGTCACCATCGAGACCGCCACTGGCGGTAATTTCGACCATCTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                                                                                                       <  6:41640/148‑1 (MQ=11)
                                              gACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGgcggcg                                                                                                      <  2:140771/149‑1 (MQ=37)
                                              gACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                                                                                                       <  2:149919/148‑1 (MQ=17)
                                                ccGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTg                                                                                                    <  2:178754/149‑1 (MQ=17)
                                                 cGTCAGCGTCACCATCGAGAACGCCACAGGCGATAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGa                                                                                                   >  8:173214/1‑149 (MQ=12)
                                                     aGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCAc                                                                                               <  2:233133/149‑1 (MQ=39)
                                                     aGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCa                                                                                                <  7:56489/148‑1 (MQ=17)
                                                     aGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCAc                                                                                               <  2:239539/149‑1 (MQ=17)
                                                        gTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGt                                                                                                     <  2:226178/140‑1 (MQ=17)
                                                        gTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACAACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGt                                                                                                     >  1:226178/1‑140 (MQ=17)
                                                             cATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCg                                                                                       >  7:27260/1‑149 (MQ=17)
                                                               tCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTg                                                                                     >  7:259751/1‑149 (MQ=39)
                                                                 gagaACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGa                                                                                    <  3:46148/148‑1 (MQ=17)
                                                                 gagaACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGAc                                                                                   <  3:10758/149‑1 (MQ=17)
                                                                      cGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTTACCGAAGGCGGCGTGATCACTTACACCGTGACCCTgt                                                                              <  2:292852/149‑2 (MQ=12)
                                                                      cGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGa                                                                              <  6:325477/149‑1 (MQ=17)
                                                                        ccACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCAc                                                                                                                 <  6:298699/112‑1 (MQ=12)
                                                                         cacaGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCTTGAgga                                                                           >  6:292412/1‑147 (MQ=12)
                                                                         cacaGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTAGAAACAACTAACGCGACACTGAAGGCGAGCACGTCGGTCACCGAAGAAGGCGTGATCACTTACACAGAGACCCTCAGc                                                                            >  3:266571/1‑148 (MQ=11)
                                                                          acaggcggTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCaa                                                                          <  5:185470/146‑1 (MQ=17)
                                                                          aCCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCaa                                                                          <  8:259751/149‑1 (MQ=39)
                                                                            aggcggTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAAcc                                                                        <  7:16373/148‑1 (MQ=18)
                                                                            aggcggTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAAcc                                                                        >  8:113843/2‑149 (MQ=18)
                                                                                ggTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGAGACCCTGAGCAACCCTGc                                                                    >  2:36072/1‑149 (MQ=21)
                                                                                 gTATTTTCGAGCTAGTGAGCGCGAATCGGACCCCGGCTCAGACCTCGATCAACGACTCGGTCGACTCCACTTCTGAGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGcc                                                                   <  4:266571/149‑1 (MQ=11)
                                                                                 gTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGcc                                                                   <  2:12394/149‑1 (MQ=17)
                                                                                 gTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGcc                                                                   <  8:111307/149‑1 (MQ=21)
                                                                                  tAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGcc                                                                   >  4:38604/1‑148 (MQ=17)
                                                                                   aaattCGAGCAACTGCCCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCCCCACCGCGACCCTGACCGCGAGCCCCACCGTCCCCCCAGGCGGCGTGATCACTTACACCCAGACCCTGAGCAACCCAGcc                                                                   >  6:156963/4‑147 (MQ=11)
                                                                                    aTTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCa                                                                 <  6:169730/148‑1 (MQ=21)
                                                                                    aTTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCa                                                                 <  1:179281/148‑1 (MQ=21)
                                                                                     tttCGAGCAACTGACCGCGTATCCGACCCCGGCTCAGACCACGATCAACGAGTCGGTCGATGCCAGCACCGCGACCCTGACGGCGAGCCCGTCGGTCTCCGTAGGGGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAg                                                                <  5:156963/148‑1 (MQ=11)
                                                                                         gAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGc                                                            <  5:309730/148‑1 (MQ=21)
                                                                                          aGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCg                                                          >  5:127180/1‑149 (MQ=17)
                                                                                          aGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCg                                                          <  5:152153/149‑1 (MQ=39)
                                                                                           gCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCgg                                                         >  1:62794/1‑149 (MQ=39)
                                                                                            cAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGt                                                        <  2:62794/149‑1 (MQ=39)
                                                                                              aCTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGtga                                                      <  3:245960/149‑3 (MQ=2)
                                                                                              aCTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTaa                                                      <  7:173214/149‑1 (MQ=25)
                                                                                                tGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAAcc                                                    >  1:323304/1‑149 (MQ=38)
                                                                                                  aCCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGAcc                                                    <  3:38604/147‑1 (MQ=2)
                                                                                                  aCCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACAACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAATGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCAGTGAcc                                                    <  6:99165/147‑1 (MQ=11)
                                                                                                  aCCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGt                                                  >  3:318490/1‑149 (MQ=255)
                                                                                                   ccccGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGt                                                  >  2:108272/1‑148 (MQ=255)
                                                                                                    cccGAATCCGACCCCGGCTCAGACCACGATCAACTACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTg                                                 >  8:41826/1‑148 (MQ=255)
                                                                                                    cccGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGAcaccactac                                                                                                                                               <  8:48362/54‑1 (MQ=255)
                                                                                                    cccGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGa                                                >  8:199823/1‑149 (MQ=12)
                                                                                                    cccGAATCCGACCCCGGCTCAGACAACGATCAACGACTCGGTCGAcaccactac                                                                                                                                               >  7:48362/1‑54 (MQ=11)
                                                                                                     ccGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGaa                                               >  6:162889/1‑149 (MQ=255)
                                                                                                     ccGAATCCGACCCCGGCTCAGACCACGATCAACGACTAGGTCGACAACACTACCGCGACCCTGACGGCGATCCCGTCGCTCACCGAAGGCGGCGTGATCACTTACACCATGACCCTGAGCAACACTTCCCAGCCGCCGGTAACCGtca                                                >  4:130119/1‑146 (MQ=14)
                                                                                                      cGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGCCCGTGacc                                              >  2:76758/1‑147 (MQ=255)
                                                                                                        aaTCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCt                                            <  4:269202/149‑1 (MQ=255)
                                                                                                          tCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGCGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGt                                          <  1:36072/149‑1 (MQ=255)
                                                                                                                cccGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGAACCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACg                                    >  2:298653/1‑149 (MQ=255)
                                                                                                                 ccGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACgg                                   <  8:27260/149‑1 (MQ=37)
                                                                                                                 ccGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCCTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACgg                                   >  7:150617/1‑149 (MQ=255)
                                                                                                                  cGGATCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTTATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGc                                  <  6:85616/149‑1 (MQ=255)
                                                                                                                          aCCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGAt                          >  1:196995/1‑149 (MQ=255)
                                                                                                                            cACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAAc                                     >  6:307845/1‑136 (MQ=37)
                                                                                                                            cACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAAc                                     <  5:307845/136‑1 (MQ=37)
                                                                                                                            cACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTa                        >  7:157032/1‑149 (MQ=255)
                                                                                                                                atgaaCGACTCGGTCGATTCCACCCGCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGt                    <  8:248295/146‑1 (MQ=11)
                                                                                                                                 tCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCCGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGtt                   >  3:275133/1‑149 (MQ=255)
                                                                                                                                   aaagACTCGGTCGAAACCACCACCGCGACCCTGACGTCGAGCCCGTCGGTCACCGAAGTCGGCGTGATCAATTACACCGTTACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGTCCAAACGATTACCGTTGa                 <  8:156686/146‑1 (MQ=14)
                                                                                                                                   aaCGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTg                  >  7:218036/1‑148 (MQ=255)
                                                                                                                                    aCGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGAGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGATGaa                >  5:222548/1‑149 (MQ=255)
                                                                                                                                     cGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAg               >  6:281244/1‑149 (MQ=39)
                                                                                                                                       aCTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCg             >  8:171838/1‑149 (MQ=2)
                                                                                                                                       aCTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCg             <  1:108272/149‑1 (MQ=255)
                                                                                                                                        cTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg            <  4:318490/149‑1 (MQ=255)
                                                                                                                                           ggTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGt          <  7:113843/148‑1 (MQ=255)
                                                                                                                                           ggTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGAACCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCggg           <  1:298653/147‑1 (MQ=255)
                                                                                                                                            gTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTa         >  4:273720/1‑148 (MQ=39)
                                                                                                                                              cGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTaaaa      >  7:235112/1‑149 (MQ=255)
                                                                                                                                               gACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACt                                                                                              >  6:171608/1‑60 (MQ=31)
                                                                                                                                               gACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACt                                                                                              <  5:171608/60‑1 (MQ=31)
                                                                                                                                               gACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAAc     >  5:283373/1‑149 (MQ=39)
                                                                                                                                                  accactaccGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCa  <  8:218036/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCA  >  NC_002947/195483‑195777

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: