Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,202,981 C→A 5.6% L430M (CTG→ATG) ‡ PP_4583 → putative Peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,202,9810CA 5.6% 234.3 / 3.4 72L430M (CTG→ATG) ‡PP_4583putative Peptidase
Reads supporting (aligned to +/- strand):  ref base C (37/31);  new base A (2/2);  total (39/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-01

AAGCGCTCAACCGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCCTGCTGGAGCGCGGCGTGGCCTTT  >  NC_002947/5202835‑5203127
                                                                                                                                                  |                                                                                                                                                   
aaGCGCTCAACCGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTg                                                                                                                                                   >  2:204024/1‑149 (MQ=255)
   cgcTCAACCGCCCGGCCCAGGTGCACCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTg                                                                                                                                                >  6:318334/1‑149 (MQ=255)
        aaCCGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGAt                                                                                                                                          >  8:282417/1‑149 (MQ=255)
        aaCCGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCAGGGTGCCgtc                                                                                                                                          <  4:99013/149‑3 (MQ=255)
        aaCCGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCAGGGTGCCgtc                                                                                                                                          >  3:99013/1‑147 (MQ=255)
         aCCGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGAt                                                                                                                                          >  3:148289/1‑148 (MQ=255)
          ccGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGAt                                                                                                                                          >  2:92324/1‑147 (MQ=255)
          ccGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCa                                                                                                                                        <  1:204024/149‑1 (MQ=255)
           cGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCAGGGTGCCgtc                                                                                                                                          >  3:132784/1‑144 (MQ=255)
           cGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCAGGGTGCCgtc                                                                                                                                          <  4:132784/146‑3 (MQ=255)
            gCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGc                                                                                                                                      >  7:197886/1‑149 (MQ=255)
            gCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGc                                                                                                                                      <  6:110366/149‑1 (MQ=255)
             cccGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGc                                                                                                                                              <  7:34583/140‑1 (MQ=255)
             cccGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGc                                                                                                                                              >  8:34583/1‑140 (MQ=255)
               cGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTg                                                                                                                                   <  4:148289/149‑1 (MQ=255)
               cGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTg                                                                                                                                   <  8:197886/149‑1 (MQ=255)
                     aGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTAcgc                                                                                                                             >  8:227611/1‑149 (MQ=255)
                       gTGCGCCAACTGACGTTGCCCACGGGCGCCAAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTAcgccg                                                                                                                           >  4:280729/1‑149 (MQ=255)
                            ccAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACgccgccgcc                                                                                                                      <  4:217474/149‑1 (MQ=255)
                             cAACTGACGTTGCCCACGGGCGCCCAGGTGTTGCTCAAGCAATCCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCa                                                                                                                     <  7:227611/149‑1 (MQ=255)
                              aaCTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAg                                                                                                                    >  2:162016/1‑149 (MQ=255)
                                 tGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCa                                                                                                                     <  6:288258/145‑1 (MQ=255)
                                 tGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCa                                                                                                                     >  5:288258/1‑145 (MQ=255)
                                      ttGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTggg                                                                                                            <  1:144349/149‑1 (MQ=255)
                                           cccGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGa                                                                                                       <  1:79412/147‑1 (MQ=255)
                                           cACGGGCGCCCAGGTGGTGATCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGa                                                                                                       >  5:4166/1‑149 (MQ=255)
                                              gggCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCg                                                                                                    >  2:137684/1‑149 (MQ=255)
                                                  gcCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGc                                                                                                 >  3:190893/1‑148 (MQ=255)
                                                   cccAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCa                                                                                               <  6:43928/149‑1 (MQ=255)
                                                         tggtgCTCAAGCAAACCCCGGTGCCCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCActgtac                                                                                         <  7:282417/149‑1 (MQ=255)
                                                              cTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCActgtacctgt                                                                                     <  3:280729/148‑1 (MQ=255)
                                                                 aaGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGTACGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCActgtacctgtac                                                                                   <  1:92324/148‑1 (MQ=255)
                                                                 aaGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACgg                                                                                 <  1:162016/149‑1 (MQ=255)
                                                                   gCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGc                                                                                                           <  8:191739/121‑1 (MQ=255)
                                                                   gCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGc                                                                                                           >  7:191739/1‑121 (MQ=255)
                                                                       aCCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATgg                                                                           >  6:318196/1‑149 (MQ=255)
                                                                           cGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGcc                                                                       >  6:150342/1‑149 (MQ=255)
                                                                              tGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTAt                                                                    <  4:190893/149‑1 (MQ=255)
                                                                                cTCGGTGCATTCGATGCCGATGACTATGTCAGAGAGCGCCTGTGGCCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATgg                                                                  >  7:99351/1‑149 (MQ=255)
                                                                                         ttCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCt                                                         >  6:192627/1‑149 (MQ=255)
                                                                                         ttCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCt                                                         <  1:137684/149‑1 (MQ=255)
                                                                                            gATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCg                                                       <  5:190100/148‑1 (MQ=255)
                                                                                             aTGCCGATGACTATGTCATTGAGCTCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGTCTATGGTGCCTATGGCGAAAGCCTCGAc                                                     <  1:13039/149‑1 (MQ=255)
                                                                                                 cGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCt                                                                                                               >  5:7196/1‑87 (MQ=255)
                                                                                                 cGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCt                                                                                                               <  6:7196/87‑1 (MQ=255)
                                                                                                 cGATGACTATGTCAGTGAGCGACTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCt                                                 >  2:72226/1‑149 (MQ=255)
                                                                                                  gATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTg                                                >  7:67927/1‑149 (MQ=255)
                                                                                                   aTGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTgg                                               >  6:8551/1‑149 (MQ=255)
                                                                                                       cTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTtct                                           <  8:67927/149‑1 (MQ=255)
                                                                                                       cTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACACTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTtct                                           >  5:312729/1‑149 (MQ=255)
                                                                                                         aTGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTca                                         >  4:69473/1‑149 (MQ=255)
                                                                                                           gTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGAc                                                     >  1:73037/1‑135 (MQ=255)
                                                                                                           gTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGAc                                                     <  2:73037/135‑1 (MQ=255)
                                                                                                            tCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTcaca                                       <  5:318196/148‑1 (MQ=255)
                                                                                                                    cgcCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCt                              <  6:78408/149‑1 (MQ=255)
                                                                                                                     gcCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCt                              >  3:238646/1‑148 (MQ=255)
                                                                                                                     gcCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTg                             <  4:238646/149‑1 (MQ=255)
                                                                                                                       cTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACAt                                      >  5:33521/1‑138 (MQ=255)
                                                                                                                       cTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACAt                                      <  6:33521/138‑1 (MQ=255)
                                                                                                                       cTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAg                           >  1:273283/1‑149 (MQ=255)
                                                                                                                        tgtgGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGc                          >  1:17221/1‑149 (MQ=255)
                                                                                                                          tgGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCct                        >  2:227938/1‑149 (MQ=255)
                                                                                                                              cAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGa                            >  1:12984/1‑141 (MQ=255)
                                                                                                                              cAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGa                            <  2:12984/141‑1 (MQ=255)
                                                                                                                              cAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCctgctg                    <  6:312729/149‑1 (MQ=255)
                                                                                                                               aaCCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCctgctg                    >  5:97219/1‑148 (MQ=255)
                                                                                                                               aaCCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCctgctg                    >  3:177902/1‑148 (MQ=255)
                                                                                                                               aaCCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCGGGTTCTCACCTGCGCGCCTGAGCCTGCTgg                   >  5:92677/1‑149 (MQ=255)
                                                                                                                                  cGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACgg                                                                                 <  3:258274/84‑1 (MQ=255)
                                                                                                                                  cGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACgg                                                                                 >  4:258274/1‑84 (MQ=255)
                                                                                                                                  cGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCCTGCTGGAgc                <  6:4166/149‑1 (MQ=255)
                                                                                                                                              cACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCCTGCTGGAGCGCGGCGTGGCCt    >  3:267369/1‑149 (MQ=255)
                                                                                                                                                cccTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCCTGCTGGAGCGCGGCGTGGCCttt  <  4:177902/149‑1 (MQ=255)
                                                                                                                                                 ccTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTtctc                                          >  6:207872/1‑108 (MQ=255)
                                                                                                                                                 ccTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTtctc                                          <  5:207872/108‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AAGCGCTCAACCGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCCTGCTGGAGCGCGGCGTGGCCTTT  >  NC_002947/5202835‑5203127

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: