Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 353,105 T→G 13.7% intergenic (+69/+49) hisF → / ← cbcV imidazole glycerol phosphate synthase subunit HisF/choline/betaine/carnitine ABC transporter ATP binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947353,1050TG13.7% 377.4 / 9.2 164intergenic (+69/+49)hisF/cbcVimidazole glycerol phosphate synthase subunit HisF/choline/betaine/carnitine ABC transporter ATP binding protein
Reads supporting (aligned to +/- strand):  ref base T (75/62);  new base G (12/15);  total (87/77)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 187 total aligned reads displayed.

CAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTT  >  NC_002947/352958‑353254
                                                                                                                                                     |                                                                                                                                                     
caacATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCtgt                                                                                                                                                     >  2:1827331/4‑150 (MQ=255)
cAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCtgt                                                                                                                                                     <  1:1789169/150‑1 (MQ=255)
cAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCtgt                                                                                                                                                     <  2:204196/150‑1 (MQ=255)
cAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCAACAGGTTTGCTGAGATCTCtgt                                                                                                                                                     >  2:2468446/1‑150 (MQ=255)
 aGCATCGTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCtgtg                                                                                                                                                    <  1:889928/150‑1 (MQ=255)
  gCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTgg                                                                                                                                                   <  1:3759239/150‑1 (MQ=255)
  gCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTgg                                                                                                                                                   <  2:1019313/150‑1 (MQ=255)
  gCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTgg                                                                                                                                                   <  1:1149655/150‑1 (MQ=255)
  gCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTgg                                                                                                                                                   <  1:494754/150‑1 (MQ=255)
  gCATCTTCCACGTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTgg                                                                                                                                                   <  2:5505328/150‑1 (MQ=255)
   cATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGGACCCCCCCGCCCCCAGGGTTGCGGTGTTCTCGGTggg                                                                                                                                                  >  1:1494652/1‑150 (MQ=255)
       ttCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCACCCACAGGTTTGCTGAGATATCGGTGGGAGCg                                                                                                                                              >  1:3128673/1‑150 (MQ=255)
        tCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCgg                                                                                                                                             <  1:5493579/150‑1 (MQ=255)
        tCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCgg                                                                                                                                             >  1:5050979/1‑150 (MQ=255)
          cACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCGACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGc                                                                                                                                           <  2:5284749/150‑1 (MQ=255)
           aCTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCg                                                                                                                                          <  1:4024799/150‑1 (MQ=255)
                ggCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCGGTGGGAGCGGGCGTGccc                                                                                                                                     <  1:2990297/150‑1 (MQ=255)
                     gTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGTTTTGCTGAGATCTCTGTGGGAGCAGGCGTGCCCGCGaa                                                                                                                                <  2:590821/150‑1 (MQ=255)
                      tACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCTGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCAGACAAAGGTTTTGTTTGATATTTGTGGGAGCGGGCGTTGCcgcgtcg                                                                                                                               >  1:4702563/1‑147 (MQ=37)
                      tACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCGCCCACAGGTTTGCTGAGATCTAGGTGGGAGCGGGCGTGCCCCCGCAg                                                                                                                               >  1:1192178/1‑150 (MQ=39)
                        cacaGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATACATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGc                                                                                                                                       >  2:2913457/1‑140 (MQ=255)
                        cacaGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATACATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGc                                                                                                                                       >  2:2913605/1‑140 (MQ=255)
                        cacaGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATACATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGc                                                                                                                                       <  1:2913457/140‑1 (MQ=255)
                        cacaGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATACATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGc                                                                                                                                       <  1:2913605/140‑1 (MQ=255)
                         acaGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGc                                                                                                                            >  2:4317620/1‑150 (MQ=255)
                         acaGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGGTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGc                                                                                                                            >  2:4700467/1‑150 (MQ=255)
                           cctGCCGGAGCCCAACGCCTCCAGCGCTTCCCGCGGGAACGCCGCTCCCGCACCCAGTGACCCGCGCCGCGTTGCAGCCCTACGCGGCCACGCCCGCACCCACAGGTTTGCTCAGAACTCGGGCGGACCGGGCGTGCCCGCGAAGGGCcg                                                                                                                          <  2:1680829/148‑1 (MQ=255)
                               ccGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACCCCCCCGCCCCCCCGTGTTGTGGGGTCTCGGTGGGTGCGGGGGTGCCCGCGAAGGGCcgcgca                                                                                                                      >  1:3829185/1‑148 (MQ=255)
                                cGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACCCCCCCGCCCACCCGGTTGGTGGGGTCGCGGTGGGAGCGGGCGGGCCCGCGAAGGGCCGCGAAg                                                                                                                     >  1:2541304/1‑150 (MQ=255)
                                           gCCGACATGGCGTCGCGCGGGATCGTCGTTCGCGGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAAtt                                                                                                          <  2:4601328/150‑1 (MQ=255)
                                              tACATGGCTTCGCGCGGGATCGTCGTTCGCTGTTCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGGAGGCCCGCTCCCACAGGTTTTGTGTGGTATGTGTGGGAGCGGGCGTGGCCGCGGAGGGGCGCGAAGCGGCCCCAATTggt                                                                                                       >  1:1722533/1‑148 (MQ=37)
                                               acttGGCTGCGCGCGGGATCGGCGTTCGCTGATCCTTTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACCGGTTGGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCtt                                                                                                      <  2:28836/147‑1 (MQ=38)
                                               aCATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGGCCGCCCCCCCCCACACGGTTGGGGAGATATAGGTGGGTGCGGGCGTGCCCCCCAAGGGCCGCCAAGCGGCCCCCAATGCtt                                                                                                      >  1:1449299/1‑150 (MQ=255)
                                                      ttCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGAGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGt                                                                                               >  2:1140676/1‑150 (MQ=255)
                                                      ttCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCGTAGATCTCGGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGt                                                                                               >  2:1997857/1‑150 (MQ=255)
                                                      ttCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCCCACCCACAGGGTTGCTGAGATATCTGGGGGAGCGGGCGTGCCCGCGCAGGGCCGCGAAGCGGCCCCAAATGCTTCaaacct                                                                                               >  1:3550962/1‑146 (MQ=38)
                                                        cgcgcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTgg                                                                                             <  2:3758040/150‑1 (MQ=255)
                                                         gcgcgGGATCGTCGTTCGCTGAAACATTGACCGGGGCTGCTTTGCAGACCTTAGCCGGCAACCACGCACCCACAGGTTTGCTGAGAACACATGTCGATCGCGCGTCACCGCGAACGTCCCCCAAGCGCCCCCAATTGCGTCAACAGTGGt                                                                                            >  2:2824418/1‑150 (MQ=255)
                                                           gcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGGTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTc                                                                                          >  2:3243294/1‑150 (MQ=255)
                                                             gggATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCtt                                                                                        <  2:4373033/150‑1 (MQ=255)
                                                               gATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGc                                                                                      <  2:4255703/150‑1 (MQ=255)
                                                               gATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGc                                                                                      <  2:1264911/150‑1 (MQ=255)
                                                                   ggggtTCGCTGATCCATGGCCCGGGGCTGGTTTGCAGCCCTGCGCGGGCCCGCCCGCCCCCACCGGTTGGGGCAGCTCTGTGGGGGAGGGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGag                                                                                  <  1:5213480/147‑1 (MQ=25)
                                                                       ttCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACCGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGagcagc                                                                              <  2:1754937/150‑1 (MQ=255)
                                                                         cGCTTACCCATGGCCCGGGGGTGCTTTGCCCCCCGTCGGGGGCCCGCCCGCCCCCACTGGTTGGCTGTGCGCGCGGGGGGCGCGGGCGGCCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAgcgc                                                                            <  1:496528/150‑1 (MQ=255)
                                                                         cGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGTTCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAAATGCTTCCACCGGGGTTCGTGCCGAGCAgcgc                                                                            >  2:4378571/1‑150 (MQ=38)
                                                                          gCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAAt                                                                                                           <  2:4129411/118‑1 (MQ=255)
                                                                           cTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGt                                                                          <  2:655323/150‑1 (MQ=255)
                                                                                ccATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACCCCCCCACCCACAGGTTTGCGGGGATCTCGGGGGGTGCGGGCGGGCCCGCGAAGGGCCGCGAAGCGGCCCCAAATGCTTCAACCGGGGGTCTTGCCGGGCAGCGCGTGGTaa                                                                     >  2:1419314/1‑150 (MQ=37)
                                                                                     gACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTGGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTc                                                                <  1:4194902/150‑1 (MQ=255)
                                                                                      aCCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACCGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg                                                               <  2:887912/150‑1 (MQ=255)
                                                                                      aCCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGCTCTCGGGGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg                                                               <  1:2294988/150‑1 (MQ=255)
                                                                                         ggggCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCCCAGGGTTGCTGAGATCTCGGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCCATTGCGTCAACCGTGGTTCTTGCCGAGCCGCGCGTGGTGAAGCGCGGtt                                                            >  2:3999384/1‑149 (MQ=37)
                                                                                          gggCTGCTTTGCAGCCCTTCGCGGGGACGCCCGCTCCCACAGGTTTGGTGTTTTGTGTGTGGGAGCGGGCGTGCCCGCGTAGGGGGGGGTAGCGGGCCCAATTGGTTGAAGCGTGGTTGTTGGCGGGGAGGGGGTGGTAAAGGTGGGttt                                                           >  2:1397441/1‑148 (MQ=16)
                                                                                           ggCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACCGGTTTGCTGGGATCTCGGTGGGGGCGGGCGGGGCCGCGAAGGGCCGCGAAGCGGCCCCCAATGCTTCAACCGGGGGTCTTGCCGAGCAGCGCGTGGGAAAGCTCGGGATc                                                          >  1:2190211/1‑150 (MQ=37)
                                                                                              tgctTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGGTGTGTTGTGTGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTGTTGGCGAGGAGGGGGTGGTAAAGGTCGGTATcccc                                                       >  2:3310605/1‑150 (MQ=37)
                                                                                                ggttGCCGCCCTTCGGGGGCACCCCCGCTCCCCCCGGTTGGCTGCGATCTGGGTGGGCGGGGGCGTCCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGa                                                     <  1:2478517/148‑1 (MQ=255)
                                                                                                ctTTGCAGCCCTGCGCGGGCACGCCCGCTCCCCCCGTTCTCTCCCCAAACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGa                                                     <  2:1943872/150‑1 (MQ=21)
                                                                                                 tttGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCCATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGGATCCCCGAg                                                    >  2:4911275/1‑150 (MQ=255)
                                                                                                 tttGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCGGAGATCTCTGGGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAg                                                    >  2:4878727/1‑150 (MQ=255)
                                                                                                     cAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATc                                                <  2:1494652/150‑1 (MQ=255)
                                                                                                       gCCCTCCGCGGGCACGCCCGCTCCCACCGTTTTGCTGATATCTCGGTGGTAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATccc                                              <  1:5297623/150‑1 (MQ=39)
                                                                                                       ccccTTCGCGGGCACGCCCGCTCCCACCGGTTTGCTGAGATCGCGGTGGGAGCGGGCGGGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATccc                                              <  1:2095946/149‑1 (MQ=255)
                                                                                                        cccTTCGCGGGCACGCCCGCTCCCACAGAGATCGCAGCAAAACGTGGGGAGCGCGCGTGACCGCGAAGGGCCGCGAAGCGGCCCCAATTGCATCAACCGGGGTTCTTGCCGAGCAGCGCGGGGTAAAGCTCGGGATCCCCGAGTATcccc                                             >  2:1346448/1‑150 (MQ=18)
                                                                                                           ttCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCccacc                                          >  2:2648130/1‑150 (MQ=255)
                                                                                                           ttCGCGGGCACGCCCGCTCCCACAGGGTTGCTGCGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGGGGTACAGCTCGGTATCCCCGAGTATCCccacc                                          >  2:4988946/1‑150 (MQ=255)
                                                                                                                   cACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTc                                  >  2:2242414/1‑150 (MQ=255)
                                                                                                                       cccGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTa                              >  1:1338966/1‑150 (MQ=255)
                                                                                                                        ccGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGGATCCCCACCACCCTGTCGTTAt                             >  1:1071069/1‑150 (MQ=255)
                                                                                                                         cGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATc                            >  2:1441723/1‑150 (MQ=255)
                                                                                                                         cGCACCCCCCGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATc                            <  1:1588072/150‑1 (MQ=255)
                                                                                                                          gCTCCCCCCGGGTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATcc                           >  2:3802086/1‑150 (MQ=255)
                                                                                                                            tCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGAAAAGCTCGGTATCCCCGAGTATCCCCACCACCCAGACGTTATCCTg                         >  2:1590646/1‑150 (MQ=255)
                                                                                                                              ccACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATcc                                               <  1:1274865/126‑1 (MQ=255)
                                                                                                                              ccACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTaca                       <  1:2742838/150‑1 (MQ=255)
                                                                                                                               cacaGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTgcag                      >  2:1262267/1‑150 (MQ=255)
                                                                                                                               cacaGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTgcag                      >  1:986610/1‑150 (MQ=255)
                                                                                                                                accggTTTGCTGATATCTCTGTGGGCGGGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTgcagc                     <  1:4317620/147‑1 (MQ=255)
                                                                                                                                acaGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGa                                                     <  2:2791082/118‑1 (MQ=255)
                                                                                                                                acaGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGa                                                     >  1:2791007/1‑118 (MQ=255)
                                                                                                                                acaGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTgcagc                     >  2:3450362/1‑150 (MQ=255)
                                                                                                                                acaGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTgcagc                     >  2:4019640/1‑150 (MQ=255)
                                                                                                                                       ttGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGc               <  2:4527088/150‑1 (MQ=255)
                                                                                                                                           tGAGATCTGTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCttgt           >  2:1095067/1‑150 (MQ=255)
                                                                                                                                           tGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTgca                       <  2:1288441/138‑1 (MQ=255)
                                                                                                                                           tGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCttgt           >  1:1464462/1‑150 (MQ=255)
                                                                                                                                           tGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCttgt           >  1:5529811/1‑150 (MQ=255)
                                                                                                                                           tGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCttgt           >  2:5140580/1‑150 (MQ=255)
                                                                                                                                           tGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCttgt           >  2:4414195/1‑150 (MQ=255)
                                                                                                                                            gagaTCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCttgtt          >  1:2639263/1‑150 (MQ=255)
                                                                                                                                            gagaTCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCttgtt          <  1:2792836/150‑1 (MQ=255)
                                                                                                                                                tctcTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCg      <  2:2322104/150‑1 (MQ=255)
                                                                                                                                                  tctGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGt    >  1:1047049/1‑150 (MQ=255)
                                                                                                                                                  tctGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGt    >  1:2802802/1‑150 (MQ=255)
                                                                                                                                                  tctGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGt    >  2:3247022/1‑150 (MQ=255)
                                                                                                                                                   ctGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGtt   >  2:2515117/1‑150 (MQ=255)
                                                                                                                                                   ctGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGtt   <  2:2639263/150‑1 (MQ=255)
                                                                                                                                                     tgtgGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGttt  <  2:3620895/150‑1 (MQ=255)
                                                                                                                                                     tgtgGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGttt  >  2:800215/1‑150 (MQ=255)
                                                                                                                                                     tgtgGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGttt  >  1:585533/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTT  >  NC_002947/352958‑353254

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