Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 214,014 T→A 5.6% L6507Q (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_002947214,0140TA 5.6% 529.3 / 14.7 162L6507Q (CTG→CAG) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (88/64);  new base A (4/5);  total (93/69)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.31e-01

Only 100 of 163 total aligned reads displayed.

AGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAATTCGGCGGTAGACGCCACCGGCG  >  NC_002947/213869‑214163
                                                                                                                                                 |                                                                                                                                                     
aGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTc                                                                                                                                                   <  2:3067566/150‑1 (MQ=255)
aGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTc                                                                                                                                                   <  2:2678991/150‑1 (MQ=255)
   aCTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTg                                                                                                                                                >  2:2514159/1‑150 (MQ=255)
    cTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCAGGTCGTgg                                                                                                                                               <  2:975226/150‑2 (MQ=255)
     ttCGGTCACCGACGAGCCGGGTACCCCAGGCATCGAAGGCGACCCGGTCAAGGTCACGATCAGCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGAc                                                                                                                                              <  2:4252216/150‑1 (MQ=255)
     ttCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGAc                                                                                                                                              >  2:4907328/1‑150 (MQ=255)
     ttCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGAc                                                                                                                                              >  1:2129788/1‑150 (MQ=255)
     ttCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGGCGttgt                                                                                                                                              <  2:3341824/150‑4 (MQ=255)
     ttCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCAGGTCGTggg                                                                                                                                              <  2:1160724/150‑3 (MQ=255)
       cGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGAccc                                                                                                                                            <  1:3068243/150‑1 (MQ=255)
        ggTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCt                                                                                                                                           <  2:811642/150‑1 (MQ=255)
          tCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGa                                                                                                                                         >  1:1996998/1‑150 (MQ=255)
          tCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCAGGTCGTgggccagc                                                                                                                                         <  1:803125/150‑8 (MQ=255)
              cgacgaGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGcaa                                                                                                                                     >  1:4162777/1‑150 (MQ=255)
                     ccGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGcc                                                                                                                              >  2:2320553/1‑150 (MQ=255)
                     ccGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGcc                                                                                                                              >  2:1929493/1‑150 (MQ=255)
                     ccGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGcc                                                                                                                              >  2:1929480/1‑150 (MQ=255)
                           aCCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTc                                                                                                                        <  2:2076911/150‑1 (MQ=255)
                             cccAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCGCTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCAc                                                                                                                      >  1:3420567/1‑150 (MQ=255)
                             cccAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCAc                                                                                                                      >  2:899399/1‑150 (MQ=40)
                             cccAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCAc                                                                                                                      >  1:1272167/1‑150 (MQ=40)
                              ccAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCAcc                                                                                                                     <  2:4049732/150‑1 (MQ=39)
                               cAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCa                                                                                                                    <  1:2419290/150‑1 (MQ=38)
                               cAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCa                                                                                                                    <  1:3917472/150‑1 (MQ=38)
                                  gCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCTAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCa                                                                                                                 >  2:2008890/1‑150 (MQ=37)
                                  gCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCa                                                                                                                 >  1:3197102/1‑150 (MQ=37)
                                   cAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCaa                                                                                                                <  2:1585846/150‑1 (MQ=37)
                                       gAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGcc                                                                                                            >  1:4921284/1‑150 (MQ=38)
                                              aCCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGaaa                                                                                                     <  1:2957810/150‑1 (MQ=38)
                                              aCCCGGTCAACGTCACCAACACCGCCGACCAGACCTCGCTGGCCGACAACGACAAACCGACCTTCACTGTGCACGACAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACCCCCAGGTCACCAACAAGGCCGGCGaaa                                                                                                     <  2:1237574/150‑1 (MQ=11)
                                               cccGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAAc                                                                                                    <  2:341090/150‑1 (MQ=37)
                                                ccGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAAcc                                                                                                   >  1:2930409/1‑150 (MQ=35)
                                                 cGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCa                                                                                                  >  2:1544820/1‑150 (MQ=35)
                                                   gTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAgc                                                                                                >  2:2041550/1‑150 (MQ=34)
                                                     cAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAgcgc                                                                                              >  2:4037826/1‑150 (MQ=34)
                                                       aGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAgc                                                                                                >  2:544593/1‑146 (MQ=34)
                                                       aGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAgc                                                                                                <  1:544432/146‑1 (MQ=34)
                                                        ggTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCg                                                                                                           >  1:4526832/1‑134 (MQ=34)
                                                         gTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCg                                                                                          >  2:1643140/1‑150 (MQ=35)
                                                          tCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGt                                                                                         >  2:4471276/1‑150 (MQ=35)
                                                           cACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTa                                                                                        <  2:1500865/150‑1 (MQ=37)
                                                                tCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACAccc                                                                                   >  1:3638409/1‑150 (MQ=37)
                                                                 cACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCa                                                                                  <  1:4745694/150‑1 (MQ=37)
                                                                  aCCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCac                                                                                 <  2:2148862/150‑1 (MQ=37)
                                                                    cgccgACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAg                                                                                                                           >  1:2636940/1‑106 (MQ=35)
                                                                    cgccgACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCacga                                                                               <  1:3906303/150‑1 (MQ=38)
                                                                          ccAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCa                                                                         >  1:4574186/1‑150 (MQ=255)
                                                                           cAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCaa                                                                        <  1:672579/150‑1 (MQ=255)
                                                                            aGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAg                                                                       <  1:2734505/150‑1 (MQ=40)
                                                                              aCCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCCCGACCTGGTCGTGACCCTGAGCACCAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGc                                                                     <  2:4128109/150‑1 (MQ=255)
                                                                               ccTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCg                                                                    >  1:683961/1‑150 (MQ=38)
                                                                               ccTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCg                                                                    >  1:1371770/1‑150 (MQ=38)
                                                                                cTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGGCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGa                                                                   <  1:1556643/150‑1 (MQ=38)
                                                                                cTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGa                                                                   <  2:1348398/150‑1 (MQ=38)
                                                                                 tCGGTGGCCGAGATCGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACATCGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTTCACCCACGACTCGCAAGGCgtt                                                                  <  2:459759/150‑3 (MQ=255)
                                                                                 tCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGAt                                                                  <  1:2852757/150‑1 (MQ=37)
                                                                                   ggTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGa                                                                <  1:549316/150‑1 (MQ=37)
                                                                                      ggCCGAGAACGACAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGt                                                             <  1:4872096/150‑1 (MQ=37)
                                                                                               cGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGa                                                    >  1:4642831/1‑150 (MQ=255)
                                                                                                gTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGAt                                                   >  1:1033375/1‑150 (MQ=255)
                                                                                                gTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGAt                                                   <  2:4810066/150‑1 (MQ=255)
                                                                                                   aaaCCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCt                                                <  1:585337/150‑1 (MQ=255)
                                                                                                       cGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGcc                                            >  2:4964330/1‑150 (MQ=255)
                                                                                                       cGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGcc                                            >  2:2239373/1‑150 (MQ=255)
                                                                                                        gACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTg                                               <  1:716046/146‑1 (MQ=255)
                                                                                                          ccTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGa                                         >  2:2346618/1‑150 (MQ=255)
                                                                                                          ccTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGa                                         >  2:336642/1‑150 (MQ=255)
                                                                                                           cTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCa                                           <  2:2913923/147‑1 (MQ=255)
                                                                                                           cTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGAt                                        >  1:4598220/1‑150 (MQ=255)
                                                                                                            ttCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATc                                       >  1:3065923/1‑150 (MQ=255)
                                                                                                             tCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCa                                      <  2:3887573/150‑1 (MQ=255)
                                                                                                               aCTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGc                                    >  1:2729978/1‑150 (MQ=255)
                                                                                                                 tgtgCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGccg                                  >  1:3568757/1‑149 (MQ=40)
                                                                                                                 tgtgCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCt                                  >  2:1011504/1‑150 (MQ=40)
                                                                                                                   tgCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTgg                                <  2:5065271/150‑1 (MQ=39)
                                                                                                                       cGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCAt                            >  2:1083903/1‑150 (MQ=39)
                                                                                                                        gTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCCTc                           >  2:3848097/1‑150 (MQ=255)
                                                                                                                        gTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATc                           >  2:3347307/1‑150 (MQ=39)
                                                                                                                        gTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATc                           >  1:355088/1‑150 (MQ=39)
                                                                                                                           aaCCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAAt                        >  1:838168/1‑150 (MQ=40)
                                                                                                                              cAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAATTcg                     >  2:667896/1‑150 (MQ=255)
                                                                                                                               aGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAATTcgg                    <  1:1432961/150‑1 (MQ=255)
                                                                                                                               aGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAATTcgg                    <  2:4121369/150‑1 (MQ=255)
                                                                                                                                gCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAATTcggt                   >  2:5092179/1‑149 (MQ=255)
                                                                                                                                 ccGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAATTcggcg                  <  1:1041505/150‑1 (MQ=255)
                                                                                                                                 ccGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAATTcggcg                  <  2:1394541/150‑1 (MQ=255)
                                                                                                                                  cGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAATTcggcgg                 <  2:1648647/150‑1 (MQ=255)
                                                                                                                                   gCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAATTCGGCGGt                <  2:723781/150‑1 (MQ=255)
                                                                                                                                    cTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAATTCGGCGGTa               <  1:4938714/150‑1 (MQ=255)
                                                                                                                                     agggtcACGACCAGGTCGTGGGCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAATTCGGCGGTAg              >  1:2584398/6‑150 (MQ=255)
                                                                                                                                     agggtcACGACCAGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAATTCGGCGGTAg              >  2:1798759/6‑150 (MQ=255)
                                                                                                                                       gCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAATTCGGCGGTAGAc            <  1:1544820/150‑1 (MQ=255)
                                                                                                                                         ccACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAATTCGGCGGTAGACGc          >  2:373727/1‑150 (MQ=255)
                                                                                                                                         ccACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAATTCGGCGGTAGACGc          >  2:2098798/1‑150 (MQ=255)
                                                                                                                                         ccACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAATTCGGCGGTAGACGc          >  1:3233135/1‑150 (MQ=255)
                                                                                                                                          cACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAATTCGGCGGTAGACGcc         >  1:1842815/1‑150 (MQ=255)
                                                                                                                                           aCGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAATTCGGCGGTAGACGCCa        >  1:4343542/1‑150 (MQ=255)
                                                                                                                                               ccTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAATTCGGCGGTAGACGCCACCgg    <  2:3006666/150‑1 (MQ=255)
                                                                                                                                                cTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAATTCGGCGGTAGACGCCACCGgc   >  2:4456058/1‑150 (MQ=255)
                                                                                                                                                 tGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAATTCGGCGGTAGACGCCACCGgcg  >  1:3213942/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
AGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAATTCGGCGGTAGACGCCACCGGCG  >  NC_002947/213869‑214163

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