Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,539,045 G→C 11.4% intergenic (+52/+51) nadE → / ← PP_4870 NH(3)‑dependent NAD(+) synthetase/azurin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,539,0450GC11.4% 265.9 / 13.6 98intergenic (+52/+51)nadE/PP_4870NH(3)‑dependent NAD(+) synthetase/azurin
Reads supporting (aligned to +/- strand):  ref base G (38/46);  new base C (5/7);  total (44/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTCGAA  >  NC_002947/5538898‑5539194
                                                                                                                                                    |                                                                                                                                                     
cTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGAc                                                                                                                                                     <  2:4628260/150‑1 (MQ=255)
cTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGAc                                                                                                                                                     <  2:4507750/150‑1 (MQ=255)
  gCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGc                                                                                                                                                   >  2:365504/1‑150 (MQ=255)
  gCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAAAACGAGGGTGAACCCACCCACGc                                                                                                                                                   >  2:365424/1‑150 (MQ=255)
     tGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCAtt                                                                                                                                                <  2:1677502/150‑1 (MQ=255)
      gggCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTc                                                                                                                                               <  1:2457780/150‑1 (MQ=255)
        gCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTcgc                                                                                                                                             <  1:41377/150‑1 (MQ=255)
        gCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTcgc                                                                                                                                             <  1:3909669/150‑1 (MQ=255)
        gCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCGCCCGCAAGACCCAGCACCAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTcgc                                                                                                                                             <  2:655608/150‑1 (MQ=255)
          aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACCAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                           <  2:4264681/150‑1 (MQ=255)
          aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACCAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                           <  2:4264955/150‑1 (MQ=255)
          aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                           <  1:2853735/150‑1 (MQ=25)
          aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                           <  2:3169965/150‑1 (MQ=25)
          aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                           <  1:4493638/150‑1 (MQ=25)
          aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                           <  2:3158009/150‑1 (MQ=25)
          aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                           <  1:2853732/150‑1 (MQ=25)
          aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                           <  1:5106313/150‑1 (MQ=255)
          aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                           <  2:3168976/150‑1 (MQ=255)
          aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCGACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                           <  2:4896953/150‑1 (MQ=255)
            ccGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGt                                                                                                                                         <  2:2309379/150‑1 (MQ=255)
            ccGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGt                                                                                                                                         <  1:3162551/150‑1 (MQ=255)
               gTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCAACCGCGAAGACGAGGGTGAACCCACCGACGCAGTCGCGGGTgta                                                                                                                                      >  2:4363272/1‑148 (MQ=255)
               gTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTAACCCGCGAAAACGAGGGTGAAACCACCGACGCATTAGCGGGTGaa                                                                                                                                      >  2:4363171/1‑150 (MQ=255)
                ttGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAAc                                                                                                                                     <  1:4705585/150‑1 (MQ=25)
                ttGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAAc                                                                                                                                     <  2:4431576/150‑1 (MQ=255)
                 tGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAAcc                                                                                                                                    <  1:2121703/150‑1 (MQ=255)
                 tGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAAcc                                                                                                                                    <  1:2121438/150‑1 (MQ=255)
                 tGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAAcc                                                                                                                                    <  1:2121421/150‑1 (MQ=255)
                  gACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACACGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGTGGGTGAACCCACCGACGCATTCGCGGGTGAAccc                                                                                                                                   <  2:1834839/150‑1 (MQ=255)
                    ccAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGc                                                                                                                                 <  2:1055539/150‑1 (MQ=255)
                    ccAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGc                                                                                                                                 <  2:1055605/150‑1 (MQ=255)
                    ccAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGc                                                                                                                                 <  2:4329663/150‑1 (MQ=255)
                        gCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTGAACCCCCCAGAACAGGGTGTAAACCCCCACGCAATCGCGGGGGTAccccccccg                                                                                                                             >  1:122817/1‑144 (MQ=255)
                             gTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTAAACCCCCAAAAAGAGGGTTAACCCACCGAAGCCATAGCGGGTGTAACCGCGCCTGCAgg                                                                                                                        >  2:2284461/1‑150 (MQ=255)
                                cGACACCATCGCCGCCACCTCCCGAAAGACCCCGCACACGGGCGAAGTGCCGTTGGCTCCTTGATACGGGGCCGGCAGTCGCGGTTTCACCCGCGAATACGCTGGTTAACCCACCGCCGCTTTCGCGGGTGAACCCGCTCCTGCAGGGtt                                                                                                                     <  2:4896055/150‑1 (MQ=255)
                                      catcGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCCAAAACGAGGGTGTACCCACCGACGCCATCGCGGGTGAAACCGCGCCTGCAGGGTTCGCCGt                                                                                                               >  1:5046807/1‑150 (MQ=255)
                                      catcGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCAACCCCGAAAACGAGGGTGAACCCACCGACCCAGTCGCGGGTGTACCCCCTCCAGCCCGGTTAGcccg                                                                                                               >  2:3134635/1‑148 (MQ=255)
                                      catcGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTCACCCCCGAAGAACAGGGTGGAACCCACCACCCATGAGGGGGTGTACCCCCACCAGCAGGGGTGGcccg                                                                                                               >  1:333441/1‑148 (MQ=255)
                                       atcGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCCCGAATAACAGGGTGAACACACCGACGCAGTAGCGGGTGCACCCCCACCAGCAGGGGTGGCCGTc                                                                                                              >  2:2791356/1‑150 (MQ=255)
                                                cccTACCGCACGACCCAGCACACGCGCGAACTGCCGTTCGCTCCTTGATACGTGGCCTGCAGGAGCGGGTTCACCCGCGAATAGGAGGGTGAACCCACCGCCGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGa                                                                                                     <  2:1233931/149‑1 (MQ=38)
                                                aCCTACCGCAAGACCCAGCACACGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCTAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGa                                                                                                     <  2:1233951/150‑1 (MQ=255)
                                                   tcccGCAAGACCCCGCACAAGCGCGAACTGCCGTTCGCTCCTTGTTAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCTTTGGCGGGTTAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTa                                                                                                  <  2:3079120/148‑1 (MQ=39)
                                                   tACCGCACGACCCCGCACAAGCGCGAACTGCCGTTCGCTCCTTTATAGGGGGCCTGCTGGAGGGTGTTCACCCGCGAATCCGGGGGTGAACCCACCGACGCTTTCGGGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTa                                                                                                  <  2:4307335/150‑1 (MQ=255)
                                                   tACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTGACCCCCGAAAACGATGGTGAACCCACCCACGCAAACGCGGGTGAACCCGCGCCAGCACGGTTCGCCGTCACCCATTGAATa                                                                                                  >  2:3039070/1‑150 (MQ=255)
                                                      cGCAAGACCCAGCACACGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCTAATACGCGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACtt                                                                                               <  2:4060490/150‑1 (MQ=255)
                                                      cGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTAAACCCCGAAAAAGAGGGTGTAAACACCCACGCATTCGCGGGTGAAACCCCGCCGGCAGGGGTCGCCGTCCGCCGATGATTAAtt                                                                                               >  2:2178227/1‑150 (MQ=255)
                                                       gCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTGACCCCCGAATACGAGGGTGGACCCACCCACCCAGTCGCGGGTGAACCCCCGCCGGCAGGGTTAGCCGGCAGCCGATTATTACGTg                                                                                              >  1:2984279/1‑150 (MQ=255)
                                                                aGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATAAGAGGGTGTACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCATTGATTAATTTACGATGAcc                                                                                     >  2:5106472/1‑150 (MQ=39)
                                                                  cacaCGCGCGAACTGCCGTTCGCTCCTTGATACGTGGCCTGCCGGCGCGGGTTCACCCGCGATTACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCtt                                                                                   <  2:3607448/150‑1 (MQ=39)
                                                                          cgAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGTGTTCACCCGCGAATAGGAGGGTGAACCCACCTACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtca                                                                           <  2:3169358/150‑1 (MQ=255)
                                                                           gAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCCGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtcat                                                                          >  2:2085372/1‑150 (MQ=255)
                                                                               ggCCGTTCGCTCCTTGATATGTGGCCTGCTGGCGGGGGTTCACCCGCGAATCCGCGGGTGAACCCACCGACGCTTTCGCGGGTTAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCatcatc                                                                      <  2:2158018/149‑1 (MQ=255)
                                                                                 ccGTTCGCTCCTTTATACGGGGCCTGCAGGAGGGGTTTCACCAGCTAATAGGTGGGTGAACCCACCGACGCATTCGCGGTTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCga                                                                    <  2:1798598/150‑1 (MQ=255)
                                                                                    ttCGCTCTTTTAAACGTGGCCTGCTGGAGGGGGTTCACCCGCTAATACGTGGGTGAACCCACCGACGCTTTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGAt                                                                 <  1:2825439/150‑1 (MQ=38)
                                                                                       gCTCCTTGATATGTGGCCTGCTGGAGCGGGTTCACCCGCTAATACGAGGGTGAACCCACCGACGCTTTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGAtgtg                                                              <  1:1250604/150‑1 (MQ=255)
                                                                                        cTCCTTGATACGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGACCCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTgg                                                             <  2:1522367/150‑1 (MQ=255)
                                                                                          ccTTGATCAGGGGCCTGCCGGCGGGGGTTCACCCGCTAATCCGCGGGTGACCCCCCCGACGCTTTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGcc                                                           <  2:597301/150‑1 (MQ=255)
                                                                                            ttGATAAGTGGCCTGCAGGAGCGGGTTGAGCCGCGAATACGTTGGTGTAGCCACCGACGCATTGGCGGGTGAACCCGGTGGTGCGGGGTTGGGCGTGATGCTTTGTTTTCTTGACGTTGACCTTGGGTTTCTTGATGGAGTTGTggggcg                                                         >  2:4685531/1‑146 (MQ=255)
                                                                                                      gCCCGCTGGAGCGGGTCCACCCGCGACATCTTGGGTGACCCCCCCGACGCCGTGGGGGGGGCCCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTGGACGATGCCCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCa                                               <  1:3689180/150‑1 (MQ=255)
                                                                                                        cTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCaga                                             <  2:2984279/150‑1 (MQ=255)
                                                                                                             ttcgCGGGTGAACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGaagaag                                        >  1:4910538/4‑150 (MQ=255)
                                                                                                             ggtgggggATCACCCGCTAATCCGCGGGTGACCCCACCGACGCATTCGGGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGaagaag                                        <  1:3134635/145‑1 (MQ=255)
                                                                                                             ggtgCGGGTTCCCCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGaagaag                                        <  2:4486909/147‑1 (MQ=255)
                                                                                                             ggcgCGGTTTCACCCGCTAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGaagaag                                        <  2:4477222/147‑1 (MQ=255)
                                                                                                                    gTTCACCCGCTAATACGCGGTTGAACCCACCGACGCTTTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGt                                 <  2:681087/150‑1 (MQ=255)
                                                                                                                    gTTCACCCGCTAATACGAGGGTGAACCCACCGATGCATTCGCGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGt                                 <  1:1144035/150‑1 (MQ=255)
                                                                                                                    gTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGt                                 >  2:4801032/1‑150 (MQ=255)
                                                                                                                     ttCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGAtg                                                                <  1:687155/118‑1 (MQ=255)
                                                                                                                     ttCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGAtg                                                                >  2:687155/1‑118 (MQ=255)
                                                                                                                     ttCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGAGTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCCGAAGAAGCTGTAGTc                                >  2:2543331/1‑150 (MQ=255)
                                                                                                                      tCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAg                                                    <  2:6926/129‑1 (MQ=255)
                                                                                                                      tCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAg                                                    >  1:6926/1‑129 (MQ=255)
                                                                                                                          ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                           >  2:3425010/1‑150 (MQ=255)
                                                                                                                          ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                           >  2:1966774/1‑150 (MQ=255)
                                                                                                                          ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                           >  2:4978525/1‑150 (MQ=255)
                                                                                                                          ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                           >  1:181059/1‑150 (MQ=255)
                                                                                                                          ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                           >  1:2780957/1‑150 (MQ=255)
                                                                                                                          acgcgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                           >  2:578163/2‑150 (MQ=255)
                                                                                                                           cgcgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCaccacc                          >  2:2655176/1‑150 (MQ=255)
                                                                                                                               aaTACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCGCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCt                      >  2:1754376/1‑150 (MQ=255)
                                                                                                                               aaTACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCt                      >  1:2377819/1‑150 (MQ=255)
                                                                                                                                 tACGCGGGTTAACCCACCGACGCATTCGGGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTc                    <  1:4075989/150‑1 (MQ=255)
                                                                                                                                 tACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTc                    >  1:1992614/1‑150 (MQ=255)
                                                                                                                                 tACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTc                    >  1:4396939/1‑150 (MQ=255)
                                                                                                                                    gAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAg                 >  1:2380467/1‑150 (MQ=255)
                                                                                                                                    gAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAg                 >  1:461933/1‑150 (MQ=255)
                                                                                                                                      gggTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGATATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCt               >  2:2679828/1‑150 (MQ=255)
                                                                                                                                       ggTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCtt              >  2:3550118/1‑150 (MQ=255)
                                                                                                                                            aaCCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGACGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGaa          <  2:1874070/150‑1 (MQ=255)
                                                                                                                                              cccACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAgt        >  1:3793238/1‑150 (MQ=255)
                                                                                                                                              cccACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAgt        >  1:3793308/1‑150 (MQ=255)
                                                                                                                                               ccACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAgtg       >  2:1821250/1‑150 (MQ=255)
                                                                                                                                               ccACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAgtg       >  1:1531593/1‑150 (MQ=255)
                                                                                                                                                cACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAgtgt      >  1:3109861/1‑150 (MQ=255)
                                                                                                                                                 aCCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTc     >  2:4751025/1‑150 (MQ=255)
                                                                                                                                                 aCCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACAACCGGCTTCCAGCTTGGAAGTGTc     >  2:4767136/1‑150 (MQ=255)
                                                                                                                                                 aCCCACGCATTCGCGGTTCAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTc     <  1:1615859/150‑1 (MQ=255)
                                                                                                                                                  ccGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTCg    >  2:3903713/1‑150 (MQ=255)
                                                                                                                                                    gACGCTTTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCTCCACCGGCTTCCAGCTTGGAAGTGTCGaa  <  2:1872564/150‑1 (MQ=255)
                                                                                                                                                    gACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTCGaa  >  2:4733616/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTCGAA  >  NC_002947/5538898‑5539194

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: