Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 195,634 C→T 46.2% T380T (ACC→ACT PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947195,6340CT46.2% 61.5 / 182.3 145T380T (ACC→ACTPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (41/35);  new base T (30/36);  total (73/72)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.45e-01

Only 100 of 157 total aligned reads displayed.

CGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGG  >  NC_002947/195487‑195780
                                                                                                                                                        |                                                                                                                                                  
     tCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACCACcgc                                                                                                                                               >  4:162798/1‑148 (MQ=39)
        aGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACCACCGCGAcc                                                                                                                                           <  3:162798/149‑1 (MQ=39)
               ggCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACgg                                                                                                                                    >  1:108561/1‑149 (MQ=39)
               ggCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACgg                                                                                                                                    >  7:375808/1‑149 (MQ=39)
                   aaTGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAg                                                                                                                                <  6:362503/149‑1 (MQ=17)
                    aTGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGc                                                                                                                               <  8:217109/149‑1 (MQ=37)
                     tGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCTATCATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAg                                                                                                                                <  1:126413/149‑1 (MQ=11)
                      gACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGccc                                                                                                                             >  8:355104/1‑149 (MQ=17)
                         gTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTc                                                                                                                          >  5:241792/1‑149 (MQ=37)
                              caacaaCGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCa                                                                                                                     >  4:315687/1‑149 (MQ=17)
                                   acGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCg                                                                                                                  >  8:258619/1‑147 (MQ=37)
                                   acGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGaa                                                                                                                <  4:220032/149‑1 (MQ=37)
                                   acGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGaa                                                                                                                <  1:203824/149‑1 (MQ=37)
                                    cGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAg                                                                                                               <  1:211270/149‑1 (MQ=37)
                                       ttCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcc                                                                                                            >  1:282707/1‑148 (MQ=17)
                                        tCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCTCCGAAggcgg                                                                                                           >  8:183571/1‑149 (MQ=37)
                                          gACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCACGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                                                                                                          <  7:258619/148‑1 (MQ=37)
                                          gACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCACGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACAACCGCGACCCTGACGGCGAGCCCGTCGGACACCCAAGgcgccg                                                                                                         >  1:314163/1‑149 (MQ=11)
                                              gTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACTACCGCGAcc                                                                                                                                           <  6:202814/111‑1 (MQ=12)
                                              gTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGAt                                                                                                     <  3:97253/149‑1 (MQ=17)
                                               tCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGAt                                                                                                     >  8:125930/1‑148 (MQ=39)
                                               tCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGAGGTAACGATCAACGACTAGGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGt                                                                                                        >  1:329981/1‑148 (MQ=11)
                                                cAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCg                                                                                                                         >  3:320868/1‑127 (MQ=17)
                                                cAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCg                                                                                                                         <  4:320868/127‑1 (MQ=17)
                                                 aGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCa                                                                                                   <  7:355104/148‑1 (MQ=17)
                                                 aGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCAc                                                                                                  <  8:35626/149‑1 (MQ=17)
                                                 aGCGTCACCATCCAGAACGCCACCGGCGGTAAATTCGAGCAACTGACCCCGA‑‑CTCCGACCCCGGCTGAGA‑‑CCACGATCAACGACTCGGTCGACACCACCTCCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATctc                                                                                                  <  8:116046/149‑3 (MQ=14)
                                                  gCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACt                                                                                                 <  8:375808/149‑1 (MQ=39)
                                                  gCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGACCAACTGACCCCGA‑‑ATCCGACCCCCGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACt                                                                                                 <  5:373231/149‑1 (MQ=255)
                                                     tCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTa                                                                                               >  1:169692/1‑148 (MQ=17)
                                                       aCCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTacac                                                                                            >  5:139473/1‑149 (MQ=39)
                                                       aCCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTacac                                                                                            >  6:220925/1‑149 (MQ=17)
                                                         cATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCCCGATCAACGACTCGGTCGACACCACCCCCGCGACCCTGCCGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCg                                                                                          >  5:119897/1‑149 (MQ=255)
                                                         cATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACAcc                                                                                           <  6:398069/148‑1 (MQ=17)
                                                          aTCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACAcc                                                                                           <  2:312634/147‑1 (MQ=17)
                                                          aTCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGt                                                                                         >  3:344313/1‑149 (MQ=17)
                                                           tCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGt                                                                                         >  6:124361/1‑148 (MQ=39)
                                                           tCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                                                                                                          >  3:225599/1‑131 (MQ=17)
                                                           tCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGt                                                                                         >  2:32490/1‑148 (MQ=17)
                                                             gagaACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGa                                                                                       <  3:300868/148‑1 (MQ=17)
                                                                aaCGCCACCGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTTGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGAt                                                                                                     >  3:138771/1‑131 (MQ=37)
                                                                     cacaGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGC‑‑CCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGc                                                                               >  6:193668/1‑148 (MQ=12)
                                                                        cggcggTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGAt                                                                                                     <  4:25712/123‑1 (MQ=38)
                                                                        cggcggTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGAt                                                                                                     >  3:25712/1‑123 (MQ=38)
                                                                        aggcggTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAAcc                                                                           <  3:6210/148‑1 (MQ=11)
                                                                         ggcggTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAAccc                                                                          <  5:53376/149‑1 (MQ=17)
                                                                            ggTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGc                                                                       <  2:294665/149‑1 (MQ=17)
                                                                              tAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGccc                                                                     <  6:369798/149‑1 (MQ=17)
                                                                                aTTTCGTGCAAGTGTCCCCGA‑‑ATGCGAGCCGGGCTCAGT‑‑CCTCTATCAACGTCTCGGTCGATGCCACCAGCGCGACGCTGTAGTCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCTCTTACACCGTGACCCTGAGCAACCGTGCCCa                                                                    <  4:6259/148‑1 (MQ=2)
                                                                                aTTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCa                                                                    <  3:270346/148‑1 (MQ=17)
                                                                                 tttCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGa                                                                  >  6:98977/1‑149 (MQ=17)
                                                                                 tttCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGa                                                                  >  2:187246/1‑149 (MQ=39)
                                                                                   tCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGAc                                                                 <  1:187246/148‑1 (MQ=39)
                                                                                    cGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGc                                                               >  3:167698/1‑149 (MQ=39)
                                                                                     gAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGcc                                                              >  7:11585/1‑149 (MQ=17)
                                                                                        cAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGt                                                           <  5:372975/149‑1 (MQ=21)
                                                                                         aaCTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGt                                                           >  5:43236/1‑148 (MQ=21)
                                                                                                cccGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACGCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGa                                                   <  4:212821/149‑1 (MQ=255)
                                                                                                cccGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTg                                                    >  3:212821/1‑148 (MQ=255)
                                                                                                cccGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTg                                                    >  2:321053/1‑148 (MQ=255)
                                                                                                 ccGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAAcc                                                       <  7:27670/144‑1 (MQ=40)
                                                                                                 ccGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAAcc                                                       >  8:27670/1‑144 (MQ=40)
                                                                                                  cGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTAACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGacc                                                 >  4:7380/1‑147 (MQ=255)
                                                                                                  cGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAAc                                                 <  1:321053/149‑1 (MQ=255)
                                                                                                       aTCCGACCCCGGCTCAGA‑‑CCAAGATCAACGACTCGGACGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTg                                              >  8:81749/1‑149 (MQ=255)
                                                                                                           gACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCg                                                             >  8:244434/1‑130 (MQ=17)
                                                                                                           gACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCg                                                             <  7:244434/130‑1 (MQ=17)
                                                                                                           gACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCAGGTAACCGTGAACCTGTcc                                           >  4:357121/1‑148 (MQ=255)
                                                                                                             ccccGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAAc                                        <  5:220925/149‑1 (MQ=38)
                                                                                                              cccGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACg                                       >  1:184035/1‑149 (MQ=40)
                                                                                                                     cAGA‑‑CCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCTTTAACCTGTCCAACGGCCAAAc                                <  1:227809/149‑1 (MQ=255)
                                                                                                                      aGA‑‑CCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAAc                                <  2:23821/148‑1 (MQ=37)
                                                                                                                      aGA‑‑CCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAAc                                <  6:43236/148‑1 (MQ=37)
                                                                                                                      aGA‑‑CCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCcaaa                                >  2:227809/1‑147 (MQ=255)
                                                                                                                       gA‑‑CCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGa                              >  8:69656/1‑149 (MQ=38)
                                                                                                                       gA‑‑CCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGa                              <  4:356460/149‑1 (MQ=39)
                                                                                                                        a‑‑CCACTATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGAt                             <  1:106725/149‑1 (MQ=255)
                                                                                                                        a‑‑CCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGAt                             >  1:124093/1‑149 (MQ=39)
                                                                                                                           ccACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGAtt                            <  5:193668/149‑1 (MQ=255)
                                                                                                                           ccACGATCAACGACTCGGTCGACACCACCACCGCGAGCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGAtt                            <  7:89536/149‑1 (MQ=255)
                                                                                                                            cACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTa                           <  2:108561/149‑1 (MQ=39)
                                                                                                                            cACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTa                           >  4:135289/1‑149 (MQ=39)
                                                                                                                               gATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCg                        >  1:94987/1‑149 (MQ=39)
                                                                                                                               gATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTAACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGAc                                                                 >  1:394988/1‑108 (MQ=25)
                                                                                                                                 tgaaCGAGTCGGTCGAGACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCTTGACCCTGAGCTACCCTGCCCAGACGCCTGTAACCGTGAACGTGTCCAACGGCCAAACGa                              <  4:150533/139‑1 (MQ=14)
                                                                                                                                    aCGACTCGGTCGACACCACCACCGCGACCCTGACGCCGAGCCCGTCGGTCACCGAACGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTCCCCAGACGCCCGAAACCGTGAACCTGTCCAACGCCCAAACGATTACCGTTGaa                   >  7:129209/1‑149 (MQ=14)
                                                                                                                                     cGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAg                  >  6:127597/1‑149 (MQ=255)
                                                                                                                                        cTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCTGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg               <  7:69656/149‑1 (MQ=255)
                                                                                                                                            ggTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCggg              <  2:329981/147‑1 (MQ=255)
                                                                                                                                             gTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTa            >  2:374046/1‑148 (MQ=39)
                                                                                                                                               cGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTaaaa         >  7:281138/1‑149 (MQ=255)
                                                                                                                                               cGACACCACAACAGCGACCCTGTCGGCGAGACCGTCGGGCACCGAACGCGGCGAGATCACATTCACACTCTCCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTaaaa         <  3:138487/149‑1 (MQ=11)
                                                                                                                                                   accactaccGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCAAAACGATTACCGTTGAAGCGGGTAAAAcc       >  1:173182/1‑147 (MQ=255)
                                                                                                                                                     cactaccGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTaaaa         >  2:177522/1‑143 (MQ=255)
                                                                                                                                                     cactaccGCGACCCTGACGGCGACCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAgg   <  5:382816/149‑1 (MQ=255)
                                                                                                                                                     caccaccGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCcgg              >  1:247634/1‑138 (MQ=12)
                                                                                                                                                     caccaccGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCcgg              <  2:247634/138‑1 (MQ=12)
                                                                                                                                                      accaccGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCcggc             >  1:389201/1‑137 (MQ=9)
                                                                                                                                                      accaccGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCcggc             <  2:389201/138‑2 (MQ=9)
                                                                                                                                                      accaccGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAggg  >  8:275799/1‑149 (MQ=39)
                                                                                                                                                        |                                                                                                                                                  
CGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGA‑‑ATCCGACCCCGGCTCAGA‑‑CCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGG  >  NC_002947/195487‑195780

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: