Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,446,203 T→G 6.0% A116A (GCT→GCG thiE → thiamine‑phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,446,2030TG 6.0% 229.8 / 4.0 80A116A (GCT→GCGthiEthiamine‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base T (32/43);  new base G (3/2);  total (35/45)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

TGACGCCGAACTGGCCGCCCGCCTGGGTGTCGGCGTGCACCTGGGCCAGACCGACGGCCCGCTGACCCCGGCCCGCGCACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCAC  >  NC_002947/5446059‑5446346
                                                                                                                                                |                                                                                                                                               
tGACGCCGAACTGGCCGCCCGCCTGGGTGTCGGCGTGCACCTGGGCCAGACCGACGGCCCGCTGACCCCGGCCCGCGCACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCa                                                                                                                                             <  8:67021/149‑1 (MQ=255)
 gACGCCGAACTGGCCGCCCGCCTGGGTGTCGGCGTGCACCTGGGCCAGACCGACGGCCCGCTGACCCCGGCCCGCGCACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAg                                                                                                                                            >  2:171100/1‑149 (MQ=255)
  aCGCCGAACTGGCCGCCCGCCTGGGTGTCGGCGTGCACCTGGGCCAGACCGACGGCCCGCTGACCCCGGCCCGCGCACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAg                                                                                                                                            <  1:171100/148‑1 (MQ=255)
  aCGCCGAACTGGCCGCCCGCCTGGGTGTCGGCGTGCACCTGGGCCAGACCGACGGCCCGCTGACCCCGGCCCGCGCACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGc                                                                                                                                           >  6:289390/1‑149 (MQ=255)
   cGCCGAACTGGCCGCCCGCCTGGGTGTCGGCGTGCACCTGGGCCAGACCGACGGCCCGCTGACCCCGGCCCGCGCACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCg                                                                                                                                          <  7:269182/149‑1 (MQ=255)
    gCCGAACTGGCCGCCCGCCTGGGTGTCGGCGTGCACCTGGGCCAGACCGACGGCCCGCTGACCCCGGCCCGCGCCCTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCg                                                                                                                                          >  7:250660/1‑148 (MQ=255)
     ccGAACTGGCCGCCCGCCTGGGTGTCGGCGTGCACCTGGGCCAGACCGACGGCCCGCTGACCCCGGCCCGCGCACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGaa                                                                                                                                        <  5:159301/149‑1 (MQ=255)
       gAACTGGCCGCCCGCCTGGGTGTCGGCGTGCACCTGGGCCAGACCGACGGCCCGCTGACCCCGGCCCGCGCACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAgg                                                                                                                                      <  8:316575/149‑1 (MQ=255)
       gAACTGGCCGCCCGCCTGGGTGTCGGCGTGCACCTGGGCCAGACCGACGGCCCGCTGACCCCGGCCCGCGCACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCAGCCTGGGCGGCCAgttcgag                                                                                                                                      <  4:350408/149‑7 (MQ=255)
          cTGGCCGCCCGCCTGGGTGTCGGCGTGCACCTGGGCCAGACCGACGGCCCGCTGACCCCGGCCCGCGCACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGgcgc                                                                                                                                   >  8:208053/1‑149 (MQ=255)
             gccgccCGCCTGGGTGTCGGCGTGCACCTGGGCCAGACCGACGGCCCGCTGACCCCGGCCCGCGCACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCa                                                                                                                                 <  5:228595/148‑1 (MQ=255)
                 cccgccTGGGTGTCGGCGTGCACCTGGGCCAGACCGACGGCCCGCTGACCCCGGCCCGCGCACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTAc                                                                                                                            <  5:324719/149‑1 (MQ=255)
                  ccgccTGGGTGTCGGCGTGCACCTGGGCCAGACCGACGGCCCGCTGACCCCGGCCCGCGCACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCCAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACg                                                                                                                           >  6:17472/1‑149 (MQ=255)
                    gccTCGCTGTCGCCGTGCACCTCGCCCAGACCGACGGCCCGCTGACCCCGGCCCGCGCACTCCTCGGCCCTAAGGCAATCATCGGCTCCACCTGCCACGCCCGCCTCGAAATCGCCGCCCAGGCTGCCACCGAACGCGCCAGCTACgag                                                                                                                         >  8:114059/1‑147 (MQ=255)
                          gtgtCGGCGTGCACCTGGGCCAGACCGACGGCCCGCTGACCCCGGCCCGCGCACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCttt                                                                                                                   <  6:293287/149‑1 (MQ=255)
                          gtgtCGGCGTGCACCTGGGCCAGACCGACGGCCCGCTGACCCCGGCCCGCGCACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCttt                                                                                                                   <  5:258588/149‑1 (MQ=255)
                               ggCGTGCACCTGGGCCAGACCGACGGCCCGCTGACCCCGGCCCGCGCACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCg                                                                                                              >  2:225854/1‑149 (MQ=255)
                                    gCACCTGGGCCAGACCGACGGCCCGCTGACCCCGGCCCGCGCACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGcttct                                                                                                         >  2:47434/1‑149 (MQ=255)
                                    gCACCTGGGCCAGACCGACGGCCCGCTGACCCCCGCCCGCGCACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGcttct                                                                                                         <  8:129929/149‑1 (MQ=255)
                                        cTGGGCCAGACCGACGGCCCGCTGACCCCGGCCCGCGCACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCaa                                                                                                     >  2:197013/1‑149 (MQ=255)
                                        cTGGGCCAGACCGACGGCCCGCTGACCCCGGCCCGCGCACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCaa                                                                                                     <  8:290856/149‑1 (MQ=255)
                                          gggCCAGACCGACGGCCCGCTGACCCCGGCCCGCGCACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAAtt                                                                                                   <  8:256722/149‑1 (MQ=255)
                                             ccAGACCGACGGCCCGCTGACCCCGGCCCGCGCACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTg                                                                                                <  4:203942/149‑1 (MQ=255)
                                              cAGACCGACGGCCCGCTGACCCCGGCCCGCGCACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCttt                                                                                                                   >  3:24928/1‑129 (MQ=255)
                                              cAGACCGACGGCCCGCTGACCCCGGCCCGCGCACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCttt                                                                                                                   <  4:24928/129‑1 (MQ=255)
                                               aGACCGACGGCCCGCTGACCCCGGCCCGCGCACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTc                                                                                              <  7:114059/149‑1 (MQ=255)
                                                gACCGACGGCCCGCTGACCCCGGCCCGCGCACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCa                                                                                             >  6:87281/1‑149 (MQ=255)
                                                   cGACGGCCCGCTGACCCCGGCCCGCGCACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGcttct                                                                                                         >  1:140663/1‑134 (MQ=255)
                                                   cGACGGCCCGCTGACCCCGGCCCGCGCACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGcttct                                                                                                         <  2:140663/134‑1 (MQ=255)
                                                      cGGCCCGCTGACCCCGGCCCGCGCACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGc                                                                                       <  5:289390/149‑1 (MQ=255)
                                                       ggCCCGCTGACCCCGGCCCGCGCACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGcc                                                                                      >  6:33124/1‑149 (MQ=255)
                                                        gCCCGCTGACCCCGGCCCGCGCACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCg                                                                                     >  5:327310/1‑149 (MQ=255)
                                                         cccGCTGACCCCGGCCCGCGCACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCgg                                                                                    <  1:225854/149‑1 (MQ=255)
                                                           cGCTGACCCCGGCCCGCGCACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGgc                                                                                  <  7:208053/149‑1 (MQ=255)
                                                           cGCTGACCCCGGCCCGCGCACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGgc                                                                                  <  3:270667/149‑1 (MQ=255)
                                                           cGCTGACCCCGGCCCGCGCACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGgc                                                                                  >  4:89509/1‑149 (MQ=255)
                                                                 ccccGGCCCGCACACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGccc                                                                            <  7:190214/149‑1 (MQ=255)
                                                                     ggCCCGCGCACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTg                                                                        <  5:87281/149‑1 (MQ=255)
                                                                      gCCCGCGCACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGc                                                                       <  3:89509/149‑1 (MQ=255)
                                                                        ccGCGCACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCACGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATACTGTCACCAACCCGGGCGCGCCCGCTGcc                                                                      >  8:1483/1‑148 (MQ=255)
                                                                            gcACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCACCGt                                                                 >  2:329363/1‑149 (MQ=255)
                                                                              aCTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCg                                                               <  5:71791/149‑1 (MQ=255)
                                                                                tGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGc                                                             >  3:52555/1‑149 (MQ=255)
                                                                                  cTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGccg                                                           >  5:211257/1‑148 (MQ=255)
                                                                                  cTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGc                                                             <  6:327310/147‑1 (MQ=255)
                                                                                  cTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTAACCAAGCCGGGCGCGCACGCTGCCAACGTCGGCCt                                                           >  4:294333/1‑149 (MQ=255)
                                                                                     ggCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCttt                                                                                                                   <  5:77642/90‑1 (MQ=255)
                                                                                     ggCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCttt                                                                                                                   >  6:77642/1‑90 (MQ=255)
                                                                                         gTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAg                                                    <  4:52555/149‑1 (MQ=255)
                                                                                                 atcatcGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCcg                                            <  3:208306/149‑1 (MQ=255)
                                                                                                    atcGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGcttctt                                                                                                        <  4:7814/86‑1 (MQ=255)
                                                                                                    atcGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGcttctt                                                                                                        >  3:7814/1‑86 (MQ=255)
                                                                                                    atcGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCaa                                                                    <  2:67961/122‑1 (MQ=255)
                                                                                                    atcGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCaa                                                                    >  1:67961/1‑122 (MQ=255)
                                                                                                    atcGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGccgacg                                          <  5:17472/148‑5 (MQ=255)
                                                                                                    atcGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCAGCGCGAGCACGCCCgcgc                                         >  3:327487/1‑149 (MQ=255)
                                                                                                     tcGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGcc                                        >  7:272181/1‑149 (MQ=255)
                                                                                                       ggCTCCACCTGCCCCGCCAGCCTCGAACTGGCCGCCCCGGCTGCAAGCGAAGGCGCGAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCCAGCGGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGccc                                       <  4:93426/148‑1 (MQ=255)
                                                                                                        gCTCCACCTGCCACGCCAGCGTCGAACTGGCCGCCCTGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTGTGTCACCTAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAg                                     <  7:1483/149‑1 (MQ=255)
                                                                                                        ccgccACGTGCCACGCCCGCCTCCACCTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAg                                     <  2:310139/146‑1 (MQ=255)
                                                                                                            cACCTGCCACGCCAGCCTCGAACTGGCTGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGGCGGCCTGCTCGAGCAGGCCCGCGCCCAGGt                                   <  5:272866/147‑1 (MQ=255)
                                                                                                            cACCTGCCACGCCAGCCTCGAACGGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGt                                   <  6:138918/147‑1 (MQ=255)
                                                                                                             aCCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGaa                                <  8:272181/149‑1 (MQ=255)
                                                                                                             aCCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCGGACAACGTCCGCCTGCTCGAGCAGGCCCGCGCCCAGGGGa                                 >  8:311590/1‑148 (MQ=255)
                                                                                                                 gCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGGTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCGGCCAACGTCGGCCGGGGGGCGCGGGCCCCCGCCCCGGGGaagcgg                            >  1:200723/1‑145 (MQ=255)
                                                                                                                  ccACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGc                           <  1:47434/149‑1 (MQ=255)
                                                                                                                            cTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGAGCCGGCCCGCGCCCCGGTGAAACAGCCGATCGCCGt                 >  5:323924/1‑149 (MQ=255)
                                                                                                                             cttcgcTGGCAGCCTGGGCGGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCGCGCGCCGGCCCGCGCCCAGGTGGAACTGCCGATCGCCgg                 >  7:289926/6‑147 (MQ=255)
                                                                                                                              ttcgcTGGCAGCCTGGGCGGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGCGCATGGCCGCGCCCAGGTGAAACTTCCGATCGCCGTGa               >  7:258114/5‑149 (MQ=255)
                                                                                                                              cGAACTGGCCGCCCAGGCTGCCTGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGGCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGa               <  2:249594/149‑1 (MQ=255)
                                                                                                                                aaCTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTg                                  <  3:206275/128‑1 (MQ=255)
                                                                                                                                aaCTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTg                                  >  4:206275/1‑128 (MQ=255)
                                                                                                                                 aCTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCACGCCGGGCGCGCCCGCTGCCAACGGCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCg            <  5:33124/149‑1 (MQ=255)
                                                                                                                                  cTGGCAGCCTGGGCGGCCAGCGAAGGAGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCgg           >  7:12575/1‑149 (MQ=255)
                                                                                                                                    ggCCGCCCAGGCGGCCCCCGAAGGCGCCCGGTTCGGGGGCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTg         <  2:200723/149‑1 (MQ=255)
                                                                                                                                     gcggccCTGGGTGCCAGCGATGGCGCAAGCTAGGTGGCCTTTGGCCGCTTTTTAAATTTTGTCACCATGCGGGGCGCGCCGGCTGCCAACGTCGGCGTGCTCTAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTg         <  4:327487/145‑1 (MQ=255)
                                                                                                                                     gccgccCAGGCTGCCAGGGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGCGCAGGCCCCCGCCCAGGTGATACTGCCGATCGCCGTGATCCGTgg        >  3:229849/1‑149 (MQ=255)
                                                                                                                                     gccgccCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGCGCAGGCCCGCGCCCCGGGGAAACAGCCGATCGCCGTGATCGGTgg        >  1:188780/1‑149 (MQ=255)
                                                                                                                                         cccAGGCTCCCCGCTACGGCCCCTGCTAGGGGGCTTTGGCCCGCTTATCCAATGCTCCCACCCCCCGGGGCGCCCCCGCGCCCACCGGGGCCGGGCCCCCGCCGCCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATc    <  2:188780/149‑1 (MQ=255)
                                                                                                                                           cAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGCGCAGGCCCCCCCCCAGGGGCAACTGCCCATCGCCCGGATCGGTGGCAtacc  >  5:337532/1‑146 (MQ=255)
                                                                                                                                                |                                                                                                                                               
TGACGCCGAACTGGCCGCCCGCCTGGGTGTCGGCGTGCACCTGGGCCAGACCGACGGCCCGCTGACCCCGGCCCGCGCACTGCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCAC  >  NC_002947/5446059‑5446346

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