Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 212,694 T→A 9.3% L6067Q (CTG→CAG)  PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947212,6940TA 9.3% 158.3 / 9.9 66L6067Q (CTG→CAG) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (25/34);  new base A (3/3);  total (29/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.81e-01

GGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTG  >  NC_002947/212545‑212843
                                                                                                                                                     |                                                                                                                                                     
ggTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCt                                                                                                                                                       >  1:655740/1‑150 (MQ=35)
  tCAGCACGTCGGTTACCGACGAGTCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTgg                                                                                                                                                     <  1:1447163/150‑1 (MQ=25)
  tCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCAgg                                                                                                                                                     <  2:925169/150‑1 (MQ=25)
  tCAGCACGGGGCTTACCGACGAGCCATGTACGGCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTgg                                                                                                                                                     <  2:109125/150‑1 (MQ=11)
    aGCACGTCGGGTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTc                                                                                                                                                   <  1:1825756/150‑1 (MQ=25)
      cACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGt                                                                                                                                                 <  2:1490393/150‑1 (MQ=35)
      cACGTCGGTTACCGACGAGCCAGGGACGCCGGGCCACGTCGGGGAGCGGGTGAAGGTCGCGATCACCGCCGACCATACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGt                                                                                                                                                 <  2:592805/150‑1 (MQ=11)
        cGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGa                                                                                                                                               <  2:1621606/150‑1 (MQ=35)
         gTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCCATCTGCTGCTGGCCCACGACCAGGTCGTggg                                                                                                                                              <  1:51040/150‑3 (MQ=11)
              ttACCGACGAGACAGTTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGCCCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGa                                                                                                                                         <  2:123576/150‑1 (MQ=16)
                  cgacgaGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGcaa                                                                                                                                     <  2:1555686/150‑1 (MQ=32)
                  cgacgaGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGAACTCGGTGGCCGAGCACGTCAAGTCGAGCGATATCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGcaa                                                                                                                                     <  1:233587/150‑1 (MQ=11)
                   gactaGCCCGGTAGGGCGGGCAACGAAGGCCGCCTGGTCAACGGCGCCGTCACCGCCGCCCAGACCTCGGTGGCCTAGAACGTCAAGCCGGCCTTCACCGTCCACCTCAACCAGCCGCGGGCCCACGACCTGGTCGTGACCCTGAGcaac                                                                                                                                    <  2:786634/150‑1 (MQ=11)
                   gacgaGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGcaac                                                                                                                                    >  2:1691285/1‑150 (MQ=32)
                   gacgaGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGcaac                                                                                                                                    <  2:655740/150‑1 (MQ=32)
                      gaGCCAGGTACGCCGGGCAACGAAGGCGACCTTGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCaacaac                                                                                                                                 >  2:1084175/1‑150 (MQ=17)
                      gaGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGATCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCaacaac                                                                                                                                 <  1:1381331/150‑1 (MQ=17)
                         ccAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGcc                                                                                                                              >  1:1577578/1‑150 (MQ=32)
                         ccAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGcc                                                                                                                              >  1:1212065/1‑150 (MQ=32)
                         ccAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGCCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGcc                                                                                                                              <  2:256544/150‑1 (MQ=17)
                         ccAGGTACGAAGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGcc                                                                                                                              >  2:78957/1‑150 (MQ=17)
                         ccAAGTACGCTGGGCAACGAAGGCGCCCTGGTCAGAGTCACGCAGTCCGCCGGCGCGCCCTCGGGGGCCGAGAACGTAAAGCCGACCTTCACCGCGCACGCCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTCAGCAACAACGcc                                                                                                                              <  2:256593/150‑1 (MQ=11)
                          cAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGccc                                                                                                                             <  1:372784/150‑1 (MQ=32)
                           aGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCa                                                                                                                            <  1:1691285/150‑1 (MQ=32)
                            ggTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACCACGCCCAg                                                                                                                           <  2:1918307/150‑1 (MQ=17)
                              tACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGt                                                                                                                         <  2:91442/150‑1 (MQ=32)
                               aCGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTc                                                                                                                        <  2:1438133/150‑1 (MQ=32)
                                cGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCa                                                                                                                       >  2:1714484/1‑150 (MQ=32)
                                cGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCa                                                                                                                       <  1:1714484/150‑1 (MQ=32)
                                        aaCGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAg                                                                                                               >  1:192692/1‑150 (MQ=32)
                                          cGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGc                                                                                                             <  2:1224036/150‑1 (MQ=32)
                                          cGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGc                                                                                                             >  1:1353275/1‑150 (MQ=32)
                                           gAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGt                                                                                                                                                 <  1:198940/113‑1 (MQ=32)
                                           gAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGt                                                                                                                                                 >  2:198940/1‑113 (MQ=33)
                                           gAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCAGGTCGTgg                                                                                                                                               >  2:1736203/1‑114 (MQ=18)
                                           gAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCAGGTCGTgg                                                                                                                                               <  1:1736203/115‑2 (MQ=18)
                                            aaGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTg                                                                                                           <  2:422168/150‑1 (MQ=32)
                                              ggCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGc                                                                                                         <  1:144736/150‑1 (MQ=32)
                                              ggCGACCTGCTCCAGGCCACCATCACCGCCGACCAGACCTCGGTGGCCGAGCACGTCAAGCCGCCCTCCACCGCGCACGCCAACCAGCCGCTGGCACACCACCCGGCCGTCACCCCGAGCAACAACGCCCACGTCCCCATCAAGGcccgc                                                                                                         >  2:33048/1‑146 (MQ=11)
                                                 gACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGaa                                                                                                      >  2:1398392/1‑150 (MQ=32)
                                                  aCCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCCACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTAAGCAACAACTCCCAGGTCACCATCAAGGCTGGCGaaa                                                                                                     >  1:124149/1‑150 (MQ=12)
                                                   ccTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAAc                                                                                                    <  1:78957/150‑1 (MQ=31)
                                                                                             cGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTa                                                          >  1:351186/1‑150 (MQ=32)
                                                                                             cGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAACCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTcca                                                          >  2:2934/1‑148 (MQ=11)
                                                                                                    gTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAacgac                                                   >  1:1340694/1‑150 (MQ=32)
                                                                                                        agacGACCTTTACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTg                                               <  2:1573093/147‑1 (MQ=17)
                                                                                                        aGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTg                                               <  2:170855/150‑1 (MQ=32)
                                                                                                               cTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGAt                                        >  2:129735/1‑150 (MQ=32)
                                                                                                               cTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGAt                                        >  1:1593960/1‑150 (MQ=32)
                                                                                                               cTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGAt                                        <  2:1212065/150‑1 (MQ=32)
                                                                                                                ttCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATc                                       >  2:1126158/1‑150 (MQ=32)
                                                                                                                     cGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCt                                  <  2:1593960/150‑1 (MQ=32)
                                                                                                                       tGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCCGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCCgg                                >  2:1131284/1‑150 (MQ=17)
                                                                                                                           ctgcAACCAGCCGCTGGCCCACGACCTGGTCGGGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTAt                            <  2:678949/147‑1 (MQ=25)
                                                                                                                            gTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATc                           >  1:1544701/1‑150 (MQ=34)
                                                                                                                              cAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCCGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCCGGCCAGATCAGCCTGGGTATCaa                         <  2:948034/150‑1 (MQ=255)
                                                                                                                                   aGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTcgg                    >  2:426960/1‑150 (MQ=37)
                                                                                                                                         tGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTgg              <  2:748095/150‑1 (MQ=37)
                                                                                                                                          gggtcACGACCAGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGa             >  1:374425/5‑150 (MQ=255)
                                                                                                                                             tcACGACCAGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGc          >  1:174514/2‑150 (MQ=255)
                                                                                                                                             ccACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGc          >  2:776179/1‑150 (MQ=38)
                                                                                                                                             ccACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGc          >  1:1228957/1‑150 (MQ=39)
                                                                                                                                             ccACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGc          >  2:1663820/1‑150 (MQ=39)
                                                                                                                                              cACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGcc         >  2:1242099/1‑150 (MQ=37)
                                                                                                                                              cACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGcc         >  2:1242049/1‑150 (MQ=38)
                                                                                                                                                cGACCTGGTCGTGACCCTGAGCAACACCGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCAc       <  2:1257877/150‑1 (MQ=38)
                                                                                                                                                   ccTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTgg    <  2:1228957/150‑1 (MQ=39)
                                                                                                                                                     tGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTg  >  2:1526385/1‑150 (MQ=38)
                                                                                                                                                     |                                                                                                                                                     
GGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTG  >  NC_002947/212545‑212843

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