Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 195,628 C→T 13.2% D378D (GAC→GAT PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947195,6280CT13.2% 193.8 / 20.6 92D378D (GAC→GATPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (36/42);  new base T (6/6);  total (43/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAAC  >  NC_002947/195480‑195774
                                                                                                                                                    |                                                                                                                                                  
gCAGGGTCTATTTCAAGACCCCGGCCAATGAAGTCTACAACAACGTTTCGACCGTCAGCGTCACCATCGAGAAGGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGac                                                                                                                                                    <  8:299566/149‑1 (MQ=11)
   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                                  <  5:307845/136‑1 (MQ=37)
                                                                                                                               cACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAAc                                  >  6:307845/1‑136 (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)
                                                                                                                                                    |                                                                                                                                                  
GCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAAC  >  NC_002947/195480‑195774

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: