Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 197,089 T→C 17.1% N865N (AAT→AAC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947197,0890TC17.1% 124.6 / 21.1 71N865N (AAT→AACPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (22/34);  new base C (4/7);  total (28/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.46e-01

CGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTA  >  NC_002947/196957‑197224
                                                                                                                                    |                                                                                                                                       
cGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCg                                                                                                                                      <  2:301733/136‑1 (MQ=34)
cGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCg                                                                                                                                      <  1:224345/136‑1 (MQ=34)
cGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCg                                                                                                                                      <  1:193999/136‑1 (MQ=34)
      aGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCg                                                                                                                                <  1:345270/136‑1 (MQ=37)
      aGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCg                                                                                                                                <  1:385928/136‑1 (MQ=37)
      aGCCGGCAAGACCCAGGGCAGCGTCGATTTCAAGACCACGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCg                                                                                                                                <  1:70547/136‑1 (MQ=255)
       gCCGGCAAGACCCGGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGGAATTTCGAGCAACTGACCCCGAATCCGACCCCgg                                                                                                                               >  1:166541/1‑136 (MQ=255)
       gCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCgg                                                                                                                               >  2:103626/1‑136 (MQ=37)
       gCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCgg                                                                                                                               >  1:156912/1‑136 (MQ=37)
       gCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGAAAATTCAAACAACTGCCACCGAAACCGACCCCgg                                                                                                                               >  1:64834/1‑136 (MQ=11)
             aaaaCCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGa                                                                                                                                     <  1:136457/121‑1 (MQ=17)
             aaGACCCAGGGCAGCGTCGATTTCCAGACCGCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGa                                                                                                                         <  1:162455/136‑1 (MQ=37)
                       gCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCaa                                                                                                               >  2:196458/1‑136 (MQ=34)
                         aGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTTACCCCGAATCCGACCCCGGCTCAGACCACGATCAACg                                                                                                             <  1:90726/136‑1 (MQ=21)
                         aGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACg                                                                                                             >  2:18056/1‑136 (MQ=34)
                         aGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGACAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACg                                                                                                             <  1:358082/136‑1 (MQ=21)
                              cGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCTTCACCATCGAGAACGCCACTGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCg                                                                                                        <  1:92167/136‑1 (MQ=17)
                                aTTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACCGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGt                                                                                                      <  1:140979/136‑1 (MQ=21)
                                   tCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGa                                                                                                   >  2:410616/1‑136 (MQ=34)
                                   tCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGa                                                                                                   >  1:70961/1‑136 (MQ=34)
                                    ccAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGAt                                                                                                  <  1:18056/136‑1 (MQ=35)
                                    ccAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGAt                                                                                                  >  2:337101/1‑136 (MQ=35)
                                        aCCCCGGCGAATGACGTCTACAACAACGGATCGCCCGTCACCCAAACCATCGAGAACCCAACAGGCGGTAAATTAGAGCAACGGACCCCGAATCCGACCCACGCTCAGACCACGATCAACGACTCGGTCGATGccc                                                                                              >  1:254278/1‑135 (MQ=11)
                                          cccGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGccacc                                                                                            >  1:118705/1‑136 (MQ=21)
                                            cGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCcaccac                                                                                          <  1:13829/136‑1 (MQ=40)
                                             ggCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCaccccc                                                                                         >  2:113960/1‑136 (MQ=255)
                                                          tACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACgg                                                                            >  1:306146/1‑136 (MQ=21)
                                                              acaacGGTTCGACCGTCAGCGTCACCAGCGACAACCCAACCGGCGTTAATTCCCACCAACGCCCCCCGAACCCCCCCCCGGCTCACACCACGACCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAg                                                                        <  2:83790/136‑1 (MQ=11)
                                                                aacGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGcc                                                                      <  2:217068/136‑1 (MQ=37)
                                                                  cGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCg                                                                    >  2:527079/1‑136 (MQ=37)
                                                                  cGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCg                                                                    >  2:11183/1‑136 (MQ=37)
                                                                   ggTTCGACCGTCAGCTTCACCATCGACAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGt                                                                   <  2:356023/136‑1 (MQ=12)
                                                                   ggTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGt                                                                   <  2:274376/136‑1 (MQ=37)
                                                                    gtgcgcccGTCAGCGTCACAATCCAGAACGCCACAGGCGGTAATTTCGAGCAACTCACCCCGAATCCGACCCCGGCTCAGCCCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTc                                                                  <  2:53007/130‑1 (MQ=11)
                                                                    gTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTc                                                                  >  2:137720/1‑136 (MQ=37)
                                                                    gTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTc                                                                  >  1:38214/1‑136 (MQ=37)
                                                                    gTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAAATTCGAGCACCTGACCCCGAAACCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGGGCTCGta                                                                  >  1:73473/1‑135 (MQ=11)
                                                                       cGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGt                                                               >  2:59770/1‑136 (MQ=37)
                                                                             tCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGa                                                         >  2:651/1‑136 (MQ=37)
                                                                                gCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAgg                                                      >  2:300983/1‑136 (MQ=37)
                                                                                  gccACCATCGAGAACGCCACAGGCGGTAATTCCGAACAACTGACCCCGAATCCGCCCCCGGCTCACACCACGACCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcg                                                    <  1:112137/134‑1 (MQ=11)
                                                                                  gTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGCCCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcg                                                    <  2:429531/136‑1 (MQ=21)
                                                                                  gTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCCACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcg                                                    <  1:300055/136‑1 (MQ=21)
                                                                                  gTCACCATCGACAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcg                                                    <  1:146697/136‑1 (MQ=21)
                                                                                    cACCTTCGACACCCCCACAGGCGTTAATTCCGCGCAACTCCCCCCTAAGCCCACCCCGGCTCAGACCACTAACAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                                                  <  1:410616/136‑1 (MQ=11)
                                                                                    cACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                                                  <  1:103626/136‑1 (MQ=34)
                                                                                    cACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                                                  <  1:159754/136‑1 (MQ=34)
                                                                                    cACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                                                  <  2:143944/136‑1 (MQ=34)
                                                                                    cACCATCGACAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                                                  <  1:239837/136‑1 (MQ=21)
                                                                                      ccATCGAGACCGCCACAGGCGGTAATTTCGACCAACTCACCCCGAATCCCACCCCGGCTCACACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGt                                                <  2:38214/136‑1 (MQ=11)
                                                                                      ccATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGCCCGCGACTCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGt                                                <  2:320636/136‑1 (MQ=255)
                                                                                      ccATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGt                                                <  2:306146/136‑1 (MQ=21)
                                                                                          cGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATc                                            >  1:464475/1‑136 (MQ=17)
                                                                                            ataaCGCCACAGGGGTTAATTTCGAGCAACTGCCCCCGAATCCGCCCCGGGCTCTGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCAc                                          <  1:501862/134‑1 (MQ=11)
                                                                                            agaACGCCACCGGGGTTAATTTCGAGCAACTCACCCCGAAACCCACCCCGGCACAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCAc                                          <  1:363315/136‑1 (MQ=11)
                                                                                            agaACGCCACAGGGGTTAATTGCCAGCAACTGACCCCTAATCCGCCCCCGGCTCACACCACTACCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCAc                                          <  1:440880/136‑1 (MQ=11)
                                                                                            agaACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCAc                                          <  1:157869/136‑1 (MQ=37)
                                                                                            agaACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCAc                                          <  2:118705/136‑1 (MQ=37)
                                                                                              aaCGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCt                                        >  2:342032/1‑136 (MQ=37)
                                                                                              aaCGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCt                                        >  2:509778/1‑136 (MQ=37)
                                                                                               aCGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTa                                       >  1:407251/1‑136 (MQ=37)
                                                                                                 gCCACTGGCGGAAATTTCGACCACCTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATTCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTaca                                     <  2:380072/136‑1 (MQ=11)
                                                                                                 cccACTGGCGTTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCATACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTaca                                     <  1:170291/135‑1 (MQ=255)
                                                                                                      aGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTg                                <  1:651/136‑1 (MQ=35)
                                                                                                      aGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTg                                <  2:407251/136‑1 (MQ=35)
                                                                                                           ttAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCt                           <  1:137720/135‑1 (MQ=21)
                                                                                                           taaaTTTCGACAAACTCCCCCCCAACCCGCCCCCGGCTCACCCCACTATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCt                           <  1:196458/134‑1 (MQ=11)
                                                                                                                     gCAACTGACCCCGCACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGAAATCACCTACACCGTGACCCTGAGCAATCCt                 >  1:524406/1‑136 (MQ=255)
                                                                                                                       aaCTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGc               >  2:544370/1‑136 (MQ=34)
                                                                                                                       aaCTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGc               >  1:149502/1‑136 (MQ=34)
                                                                                                                              cccTAATCCGACCCCGGCTCAGACCACTATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACg        <  1:267300/136‑1 (MQ=21)
                                                                                                                                    tCCGACCCCGGCTCAGACCACGATAAACGACTAGGTCGATGACACCACCGCGACCCTGACGGCGAGCCCGTCGGTCAACGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGtg  >  2:292301/1‑135 (MQ=255)
                                                                                                                                    |                                                                                                                                       
CGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTA  >  NC_002947/196957‑197224

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