Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 210,699 T→A 9.0% L5402Q (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_002947210,6990TA 9.0% 437.6 / 31.8 166L5402Q (CTG→CAG) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (89/61);  new base A (8/7);  total (98/68)
Fisher's exact test for biased strand distribution p-value = 7.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.33e-01

Only 100 of 173 total aligned reads displayed.

GTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTG  >  NC_002947/210551‑210848
                                                                                                                                                    |                                                                                                                                                     
gTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTg                                                                                                                                                      >  2:2407945/1‑150 (MQ=34)
gTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTg                                                                                                                                                      <  2:3067099/150‑1 (MQ=34)
  tggctcGTCGGTGACCGACGTGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGt                                                                                                                                                    >  2:4703559/6‑150 (MQ=255)
  tggctcGTCGGTGACCGACGTGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGt                                                                                                                                                    >  2:1927709/6‑150 (MQ=255)
   aGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTc                                                                                                                                                   <  2:4466066/150‑1 (MQ=34)
     cACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGt                                                                                                                                                 <  2:3563823/150‑1 (MQ=34)
     cACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGt                                                                                                                                                 >  1:4452336/1‑150 (MQ=34)
     cACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGt                                                                                                                                                 >  2:3405903/1‑150 (MQ=34)
       cGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCAGGTCGTGa                                                                                                                                               <  1:2250824/150‑1 (MQ=21)
       cGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCAGGTCGTGa                                                                                                                                               <  1:346201/150‑1 (MQ=21)
        gTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGGCGTGAc                                                                                                                                              >  1:1942254/1‑150 (MQ=21)
        gTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCAGGTCGTGag                                                                                                                                              <  1:4174474/150‑2 (MQ=21)
         tCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGGTCACGACCAGGTCGTGagc                                                                                                                                             <  1:2632628/150‑3 (MQ=19)
         tCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCAGGTCGTGagc                                                                                                                                             <  2:584138/150‑3 (MQ=21)
             tCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGa                                                                                                                                         >  1:3463151/1‑150 (MQ=34)
             tCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGa                                                                                                                                         <  1:2876704/150‑1 (MQ=34)
               aCCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGc                                                                                                                                       >  2:1044057/1‑150 (MQ=34)
               aCCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGc                                                                                                                                       <  1:4709595/150‑1 (MQ=34)
               aCCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTATTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGc                                                                                                                                       <  1:4709591/150‑1 (MQ=21)
                  gacgaGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCTGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGcaac                                                                                                                                    >  2:545180/1‑150 (MQ=21)
                     gaGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAAt                                                                                                                                 >  1:3178721/1‑150 (MQ=35)
                     gaGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAAt                                                                                                                                 >  2:1758940/1‑150 (MQ=35)
                     gaGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAAt                                                                                                                                 >  2:5056073/1‑150 (MQ=35)
                      aGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATg                                                                                                                                >  1:686726/1‑150 (MQ=35)
                       gCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGc                                                                                                                               >  2:4850091/1‑150 (MQ=37)
                        ccAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGcc                                                                                                                              >  1:2413954/1‑150 (MQ=37)
                          aGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCa                                                                                                                            <  2:1747111/150‑1 (MQ=37)
                             tACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGt                                                                                                                         <  1:5232883/150‑1 (MQ=37)
                               ccccAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAAt                                                                                                                                 >  2:4779759/1‑140 (MQ=33)
                               ccccAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAAt                                                                                                                                 <  1:4779786/140‑1 (MQ=33)
                                cccAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCAc                                                                                                                      >  1:1655051/1‑150 (MQ=35)
                                cccAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGAAACAATGCCCAGGTCAc                                                                                                                      >  2:518542/1‑150 (MQ=25)
                                 ccAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCAcc                                                                                                                     <  2:5076420/150‑1 (MQ=34)
                                    ggCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATc                                                                                                                  <  2:2972567/150‑1 (MQ=32)
                                       aaCGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACAAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAg                                                                                                               >  2:5437501/1‑150 (MQ=17)
                                        aCGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGcc                                                                                                                              >  1:3398827/1‑134 (MQ=32)
                                        aCGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGcc                                                                                                                              <  2:3398827/134‑1 (MQ=32)
                                          gAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGcc                                                                                                            >  1:1057460/1‑150 (MQ=32)
                                          gAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGcc                                                                                                            >  1:3498486/1‑150 (MQ=32)
                                             ggCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCAACAAGGCCGGc                                                                                                         >  1:2285547/1‑150 (MQ=17)
                                               cGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCga                                                                                                       <  1:3656546/150‑1 (MQ=32)
                                               cGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCga                                                                                                       <  1:3209905/150‑1 (MQ=32)
                                                gACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACAGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCgag                                                                                                      >  2:1440653/1‑150 (MQ=17)
                                                      gTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGATCACGCCCTGGTCGTGACCCTGAGCAACAAGGCCCAGGTCACCATCAAGGCCGGCGAGACCAgc                                                                                                <  2:5111494/150‑1 (MQ=12)
                                                           ggTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGcc                                                                                           >  2:4067123/1‑150 (MQ=32)
                                                               aCCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTac                                                                                       <  1:4109538/150‑1 (MQ=34)
                                                                 cATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTacac                                                                                     >  1:2662157/1‑150 (MQ=34)
                                                                 cATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCACTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTacac                                                                                     >  1:3767411/1‑150 (MQ=21)
                                                                       cgccgACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGc                                                                               >  2:5151002/1‑150 (MQ=32)
                                                                        gccgACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCt                                                                              <  1:88441/150‑1 (MQ=32)
                                                                        gccgACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCGGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCt                                                                              >  1:244739/1‑150 (MQ=17)
                                                                         ccgACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCAcccgg                                                                             >  1:563827/1‑146 (MQ=17)
                                                                          cgACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTgc                                                                            <  1:2207384/150‑1 (MQ=32)
                                                                            aCCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTgcgc                                                                          >  2:3060163/1‑150 (MQ=32)
                                                                              cAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAg                                                                        >  2:646470/1‑150 (MQ=32)
                                                                              cAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAg                                                                        >  1:43045/1‑150 (MQ=32)
                                                                                   ttCGGTGGCCGAGAACTTCAATCCGTCCTTCACTGTTCTCGTCAACGAGCCGCTGGCTCACGACCTGGTCGTGAGCCTGAGCAACAATGCCCAGGTCAGCATCAAGGCCTGCGTGTCCAGCGCGCCGTTCACCCACGCTGCGCAGGGCGa                                                                   <  2:3792314/150‑1 (MQ=11)
                                                                                   ttCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGa                                                                   >  1:3124190/1‑150 (MQ=32)
                                                                                   ttCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCACTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGa                                                                   <  2:3767411/150‑1 (MQ=17)
                                                                                    tCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGAt                                                                  <  1:1820663/150‑1 (MQ=32)
                                                                                     cGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATg                                                                 <  1:3589314/150‑1 (MQ=32)
                                                                                            cGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTa                                                          >  1:651418/1‑150 (MQ=32)
                                                                                            cGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGGACACCCACGCTGCGCAGGGCGATGACGTTTa                                                          >  1:651543/1‑150 (MQ=17)
                                                                                             gagaACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTAc                                                         <  1:4224590/150‑1 (MQ=32)
                                                                                             gagaACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTAc                                                         >  2:4169190/1‑150 (MQ=32)
                                                                                               gaACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACaa                                                       <  2:1124429/150‑1 (MQ=32)
                                                                                                   gTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGAt                                                   >  2:318565/1‑150 (MQ=32)
                                                                                                          cGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGcc                                            >  2:2204019/1‑150 (MQ=32)
                                                                                                          cGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGcc                                            >  2:3891878/1‑150 (MQ=32)
                                                                                                          cGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGcc                                            <  2:4952362/150‑1 (MQ=32)
                                                                                                           gACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCa                                           >  1:3362595/1‑150 (MQ=32)
                                                                                                               ttCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATc                                       >  2:2549707/1‑150 (MQ=32)
                                                                                                                tCACTGTTCACGTCAACCAGCCGCTGGCTCACGCCCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCa                                      <  1:1418019/150‑1 (MQ=17)
                                                                                                                    tGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCt                                  >  1:3852967/1‑150 (MQ=32)
                                                                                                                     gtgcACGTCTACCCGCCGCTGGCCCTGGACCTGGTCGTGACCCTGAGCAACACTGCGCAGGTCACCATCACTGCCGGCGAGTCCAGCGCGCCGTTCACCCACGCTGGGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTg                                 <  2:782193/147‑1 (MQ=11)
                                                                                                                      ttCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTgg                                >  2:5361904/1‑150 (MQ=32)
                                                                                                                           gTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATc                           >  2:1725769/1‑150 (MQ=34)
                                                                                                                            tCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCCGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCCTCa                          >  1:3406134/1‑150 (MQ=18)
                                                                                                                              aaCCAGCCGCTGGCTCACGACCTGGTCGTGACCCAGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCAGCGCAGGGCGCTGACGTTTACAACGATGCTGGCCAGATCAGCCAGGGCATCAAc                        >  1:5158800/1‑150 (MQ=255)
                                                                                                                                ccAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTc                      >  2:5304613/1‑150 (MQ=34)
                                                                                                                                 cAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTcg                     >  1:3648179/1‑150 (MQ=34)
                                                                                                                                 cAGCCGCTGGCTCACGACCTGGTCGTGACCCTAAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTcg                     >  1:1112301/1‑150 (MQ=21)
                                                                                                                                     cGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCGTGGGCATCAACTcggcgg                 >  2:3541947/1‑150 (MQ=21)
                                                                                                                                      tcagggtCACGACCAGGTCGTGAGCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGt                >  2:34678/7‑150 (MQ=255)
                                                                                                                                       cTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTgg                                >  2:330696/1‑133 (MQ=32)
                                                                                                                                       cTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTa               <  1:1436322/150‑1 (MQ=34)
                                                                                                                                        agggtCACGACCAGGTCGTGAGCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGc                                    >  2:2805528/5‑128 (MQ=255)
                                                                                                                                        agggtCACGACCAGGTCGTGAGCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGc                                    <  1:2805528/124‑1 (MQ=255)
                                                                                                                                        aggCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAg              >  2:1217228/2‑150 (MQ=21)
                                                                                                                                         gggtCACGACCAGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGa             >  1:3395695/4‑150 (MQ=21)
                                                                                                                                         ggCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTcggcg                  <  1:3931356/145‑1 (MQ=34)
                                                                                                                                          ggtCACGACCAGGTCGTGAGCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGAc            >  1:1936907/3‑150 (MQ=19)
                                                                                                                                           gtCACGACCAGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACg           >  1:4583012/2‑150 (MQ=21)
                                                                                                                                            tCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGc          <  2:2508816/150‑1 (MQ=34)
                                                                                                                                                gcccTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGAGGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACt      <  2:4452336/148‑1 (MQ=21)
                                                                                                                                                   cTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGt   >  2:1555234/1‑150 (MQ=35)
                                                                                                                                                    tGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCCCGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTg  <  2:2662237/150‑1 (MQ=25)
                                                                                                                                                    tGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTg  >  1:1271553/1‑150 (MQ=35)
                                                                                                                                                    tGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTg  >  2:6990/1‑150 (MQ=35)
                                                                                                                                                    tGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTg  >  2:3653027/1‑150 (MQ=35)
                                                                                                                                                    |                                                                                                                                                     
GTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTG  >  NC_002947/210551‑210848

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