Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,542,072 G→C 10.7% 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,542,0720GC10.7% 569.2 / 31.1 210intergenic (+52/+51)nadE/PP_4870NH(3)‑dependent NAD(+) synthetase/azurin
Reads supporting (aligned to +/- strand):  ref base G (90/91);  new base C (14/11);  total (104/106)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 215 total aligned reads displayed.

TCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTCGAA  >  NC_002947/5541923‑5542221
                                                                                                                                                     |                                                                                                                                                     
tCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCg                                                                                                                                                       <  2:6212207/150‑1 (MQ=255)
    gCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGc                                                                                                                                                   <  2:5938384/150‑1 (MQ=255)
    gCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCAACCCCGAATACGAGGGTGAACCCACCCACGc                                                                                                                                                   >  1:1056419/1‑150 (MQ=255)
       tGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCCCcgccgccat                                                                                                                                                >  2:6082852/1‑147 (MQ=255)
            aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                           <  2:3798770/150‑1 (MQ=25)
            aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                           <  1:8342/150‑1 (MQ=25)
            aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                           <  2:4752113/150‑1 (MQ=255)
            aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                           <  1:4468096/150‑1 (MQ=255)
            aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                           <  1:6061267/150‑1 (MQ=255)
            aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                           <  1:3396585/150‑1 (MQ=255)
            aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                           <  1:3396633/150‑1 (MQ=255)
            aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCAACCGCCAAAACGTGGGGGTAACCAACCACGCAATCGCgg                                                                                                                                           >  2:5803313/1‑150 (MQ=255)
              ccGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGt                                                                                                                                         <  2:3051561/150‑1 (MQ=25)
              ccGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGt                                                                                                                                         <  2:2491303/150‑1 (MQ=25)
              ccGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGt                                                                                                                                         <  1:5777961/150‑1 (MQ=255)
              ccGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGt                                                                                                                                         <  1:2311310/150‑1 (MQ=255)
                 gTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGaa                                                                                                                                      >  2:4518811/1‑150 (MQ=255)
                  ttGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAAc                                                                                                                                     <  2:1191625/150‑1 (MQ=25)
                     aCCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACACGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCg                                                                                                                                  <  1:111791/150‑1 (MQ=255)
                             gCGTTCGACATCAGCGTCGCCACCTACCGCAAGACCCAGCAAACGCGCGAACTGCCGGTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGGGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCa                                                                                                                          <  2:6690293/150‑1 (MQ=38)
                             gCGTTCGACAACATCGTCGCCACCTACCGCAAGACCCAGCAAACGCGCGAACTGCCGTTCGCTCCTTTATACGTGGCCTGCAGGAGCGGGTTCACCCGCTAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCa                                                                                                                          <  2:3875413/150‑1 (MQ=38)
                                    aCATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCAACGACGCATTCGCGGGTGAACCCGCTCCAGCAGGGTTCg                                                                                                                   >  2:5677271/1‑150 (MQ=255)
                                       tcatcGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGATCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTCACCCGCGAATACGAGGGTGGACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCg                                                                                                                >  2:6571235/1‑150 (MQ=255)
                                          tcgtcgCCACCTACCGCAAGACCCAGCACACGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCa                                                                                                             <  2:6220512/150‑1 (MQ=255)
                                          tcgtcgCCACCTACCGCAAGACCCAGCACACCCGCGAACTGCCGTGCGCTCCTTGATACGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCa                                                                                                             <  2:5572583/150‑1 (MQ=255)
                                          tcgtcgCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGGCACCCGCGAATACGAGGGTGAACCCACCGACGCATGCGCGGGGGAACCCGCGCCTGCCGGGGTCGCCggcc                                                                                                             >  2:5071119/1‑147 (MQ=255)
                                            gtcgCCACCTACCGCAAGACCCAGCACACGCGCGAACTGCCGTTCGCTCCTTGATACGGGGCCTGCTGGAGCGGGTTCACCCGCGAATATGAGGGTGAACCCACCGACGCATACGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGc                                                                                                           <  1:3950672/150‑1 (MQ=39)
                                                ccACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTGCCCCGCGAATACGAGGGTGGACCCACCCACGCAATAGCGGGTGAACCCGCGCCAGCAGGGTTAGCCGGCCGCCgtt                                                                                                       >  1:3909916/1‑150 (MQ=255)
                                                 cACCTCCCGCAAGACCCAGCACCCGCGCGAACTGCCGTTCGCTCCTTGATAAGGGGCCTGCAGGAGCGGGTTCACCCGCGAATACGCGGGTGAACCCACCGCCGCATTCGGGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTg                                                                                                      <  1:6285121/150‑1 (MQ=38)
                                                  aCCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGa                                                                                                     <  2:6738631/150‑1 (MQ=255)
                                                      cccGCACGACCCACCACACGCGCGAACTGCCGTTCGCTCCTTGATCAGGGGCCTGCAGGCGCGGGTTCACCCGCGAATACGTGGGTGAACCCACCGACGCATTCGCGGTTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTAc                                                                                                 <  1:5054575/149‑1 (MQ=255)
                                                        cGCAAGACCCCGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGTGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGGGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACtt                                                                                               <  2:5076866/150‑1 (MQ=255)
                                                         gCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATTCGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGGTCGCCGTCAGCCTTTGATTACTTg                                                                                              >  2:4445621/1‑150 (MQ=255)
                                                          cAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCCCGAAAACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCACCTGCAGGGGTCGCCGGCAGCCGTTGATTAATTGa                                                                                             >  1:6863595/1‑150 (MQ=38)
                                                           aaGACCCAGCACACGCGCGAACTGCGGTTCGATCCTTTATACGTGGCCTGCAGGAGCGGGTTCACCCGCTAATCCGAGGGTGAACCAACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGAc                                                                                            <  2:5199587/150‑1 (MQ=38)
                                                              aCCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTCACCCCCCAATACCCGGGTGTAACCACCGCCGCCAGCGCGGGTGTAACCGCACCCGCACGGGTAGCCGGCCGCCGTTGAGTACGTGAccag                                                                                         >  2:5387257/1‑147 (MQ=255)
                                                                cccgCACAAGCGCTACCTGCGGTTCGCCCCTTTATATGGGGCCTGCAGGCGCGGGTCCACCCGCTAATACGTGGTTTACCCCACCGACGATTTCGGGGGTTAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGa                                                                                       <  1:4932555/147‑1 (MQ=255)
                                                                ccAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCCCCCAAAAGAGGGTGAACCCACCCACGCATTCGCGGGTGGACCCGCGCCAGCAGGGTTCGCCCTCAGCCGTTGGTTACATGACGATGa                                                                                       >  2:3558693/1‑150 (MQ=255)
                                                                  aGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGGATACGAGGGTGAAACCACCGACGCCTTCGCGGGTGTACCCGCGCCGGCAGGGTTAGCCGGCCGCCGTTGATTACTGGTAGATGAcc                                                                                     >  1:6284685/1‑150 (MQ=255)
                                                                      caAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAAACCACCGACGCATTCGCGGGTGAACCCGCTCCAGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGc                                                                                 >  2:2731504/1‑150 (MQ=255)
                                                                           gcgAACTGCGTTTCGCTCCTTTTTAAGTGGCCTGCTGGAGGGGTTTCACCCGCTATTAGGTGGTTTAACCCACCGACGCTTTCGGGGTTTAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtc                                                                            <  2:1922493/150‑1 (MQ=255)
                                                                           gcgAACTGCCGTTCGCTCGTTTATACGGGGCCTGCTGGCGCGGGTCCACCCGCTAATACGAGGGTGAACCCACCGACGATTTCGGGGGTTAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtc                                                                            <  1:1108022/150‑1 (MQ=255)
                                                                           gcgAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCTGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtc                                                                            <  2:3903214/150‑1 (MQ=255)
                                                                            cgAACGGCCGTTCGCTCATTGATATGTGGCCTGCTGGAGCGGTTTAACCCGCTAATCCGAGGGTGAACCCACCGACGCTTTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtca                                                                           <  1:2820015/150‑1 (MQ=255)
                                                                                 tGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCatcatc                                                                      <  2:6713082/150‑1 (MQ=255)
                                                                                 ggCCGTGCGGTCCTTGATCCGGGGCCTGCTGGAGGGGGTTCACCGGCTAATACGCGGGTGAACCCACCGACGCATTGGCGGGTTAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCatcatc                                                                      <  1:6205152/149‑1 (MQ=255)
                                                                                  gCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGTTTCACCCGCTAATACTAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCg                                                                     <  1:5675110/150‑1 (MQ=255)
                                                                                     gTTCGCTCCTTTATATGTGGCCTGCAGGAGGGGGTTCACCCGCTAATACGAGGTTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCgaga                                                                  <  2:1779803/150‑1 (MQ=38)
                                                                                     gTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGTTTCACCCGCTAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCACCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCgaga                                                                  <  2:702771/150‑1 (MQ=255)
                                                                                         gCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAAACCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCGTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGAtgtg                                                              >  2:5695196/1‑150 (MQ=255)
                                                                                           tCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAAAACGAGGGTGTACCCACCCACCCATTCGCGGGTGAACCCGCTCCTGCAGGGTTAGCCGTCAGCCATTGAGTAAGTTAAGGTGACCGTGCCGTTCATCATCGAGATGTGGc                                                            >  2:4381295/1‑150 (MQ=25)
                                                                                           tCCTTGATAAGTGGCCTGCAGGAGCGGGGTCACCCCCGCAAACGTGGGGGCACCCACCCACGCAGTCGCGGGGGGCCCCCCACCAGCACGGGTGGCCGGCCGCCGTTGATTTCGTGAAGATGGCCCTGGCGTTCATCAGCGAGAAGTggg                                                            >  2:3648665/1‑149 (MQ=255)
                                                                                             cTTGATACGGGGCCGGCTGGGGGGGTTCCCCCCGCGAAGCCGCGGGTCACCCCCCCGCCGCATTCGGGGGCGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGccc                                                          <  1:86557/150‑1 (MQ=255)
                                                                                             cTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCCCGAAAAAGACGGTGAACCCCCCCCCCCCATCGCGGGTGAACCCCCGCCTGCAGGGTTCGCCGGCCCCCGTTGATTACATGACGCTTACCATGCCGTTACTCAGAGAGATGTGGccc                                                          >  2:4912311/1‑150 (MQ=255)
                                                                                                    cggggccggccGGCGCGGGTGCCCCCGCGAATACGGGGGTGACCCCACCGCCGCTTGCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAgg                                                   <  1:2170393/147‑1 (MQ=255)
                                                                                                      gggCCCGCTGGCGGGGGTCCACCGGCTAATACGCGGGTGACCCCACCGACGCATGGGCGGGTCCCCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTg                                                 <  1:5982236/149‑1 (MQ=255)
                                                                                                        gCCGGCAGGCGCGGGGTCCCCCGCGAATACGTGGGGGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTGGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCa                                               <  2:4260059/150‑1 (MQ=255)
                                                                                                         ccTGCAGGAGCGGGTTCACCCGCGAATAAGAGGGTGAAACCACCGACGCATTCGCGGGTGAACCCGCGCCGGCCGGGTTCGCCGGCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCag                                              >  2:2218915/1‑150 (MQ=38)
                                                                                                            gCAGGAGCGGGTTCACCCGCGAATACGAGGGTGACCCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCagaag                                           <  2:5665909/150‑1 (MQ=255)
                                                                                                              tggcgCGGGTTCACCGGCGAATACGCGGGTGAACCCACCGCCGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTGGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAgaagaa                                         <  2:4111217/146‑1 (MQ=39)
                                                                                                               tcggtgggTTCACCCTCGTATTCGCGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTAGTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGaagaag                                        >  2:2586375/6‑150 (MQ=255)
                                                                                                                  gCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTg                                     >  1:4398595/1‑150 (MQ=255)
                                                                                                                   gtgtttCCCCGGCTAATAAGTGGGTTAACCAACCTACTCATTCGCGGGTTAACCCGCTCCTGCAGTGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGt                                    <  1:4561407/146‑1 (MQ=255)
                                                                                                                    ggtttCCCCCGCGAATAGGGGGTTCAACCCACCGACGCTTTCGGGTGTTAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTa                                   <  1:1831042/147‑1 (MQ=255)
                                                                                                                    gggTTCACCCGCTATTACGAGGGATAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTa                                   <  2:1282875/150‑1 (MQ=255)
                                                                                                                    gggTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCTAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTa                                   <  2:1058858/150‑1 (MQ=255)
                                                                                                                      gTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAGGGTGCAGAAGAAGCTGTAGt                                 >  1:3137892/1‑150 (MQ=255)
                                                                                                                      gTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGt                                 >  1:6598325/1‑150 (MQ=255)
                                                                                                                      gTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGt                                 >  1:5926506/1‑150 (MQ=255)
                                                                                                                        tCACCCGCGAATACGAGGGTGAACCCACCGACTCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCCTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTca                               <  1:4518811/150‑1 (MQ=255)
                                                                                                                        tCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTca                               >  2:2895598/1‑150 (MQ=255)
                                                                                                                          aCCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcacc                             >  1:679237/1‑150 (MQ=255)
                                                                                                                           cccGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCGGAAGAAGCTGTAGTcacca                            >  2:834048/1‑150 (MQ=25)
                                                                                                                            ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTg                                                                                                                            >  2:2601349/1‑53 (MQ=255)
                                                                                                                            ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTg                                                                                                                            <  1:2601349/53‑1 (MQ=255)
                                                                                                                            ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGACGGTGCAGAAGAAGCTGTAGTcaccac                           >  2:3891319/1‑150 (MQ=255)
                                                                                                                            ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                           >  1:3434075/1‑150 (MQ=255)
                                                                                                                            ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                           >  1:4107092/1‑150 (MQ=255)
                                                                                                                             cgcgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCCTCGAGATGTGGCCCGGGACGGTGCAGAAGAAGCTGTAGTCaccacc                          >  2:6013363/1‑150 (MQ=255)
                                                                                                                              gagAATACGAGGGTGAACCCACCGACGCATTCGCGGTTGAACACGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCATTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGATGTAGTCACCACCg                         >  2:3744285/3‑150 (MQ=255)
                                                                                                                                      gAGGGTGAACCCACCGACTCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAg                 >  2:6388235/1‑150 (MQ=255)
                                                                                                                                      gAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAg                 >  1:5177366/1‑150 (MQ=255)
                                                                                                                                      gAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAg                 >  1:5543962/1‑150 (MQ=255)
                                                                                                                                      gAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAg                 >  1:5664124/1‑150 (MQ=255)
                                                                                                                                         ggTGAACCCCCCGCCGCATGCGCGGGTGCACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCtt              <  2:3679454/150‑1 (MQ=255)
                                                                                                                                         ggTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCtt              >  1:4340708/1‑150 (MQ=255)
                                                                                                                                          gTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTg             >  1:3189547/1‑150 (MQ=255)
                                                                                                                                          gTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTg             >  2:6265080/1‑150 (MQ=255)
                                                                                                                                             accccACCGACGCTTTCGCGGGTCACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGaa          <  1:3612954/148‑1 (MQ=255)
                                                                                                                                             aaCCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGaa          >  2:3011068/1‑150 (MQ=255)
                                                                                                                                               cccACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAgt        <  2:3252265/150‑1 (MQ=255)
                                                                                                                                                ccACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAgtg       >  2:281466/1‑150 (MQ=255)
                                                                                                                                                  aCCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCagaaga                                          <  2:5022619/113‑1 (MQ=255)
                                                                                                                                                  aCCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTc     >  2:4621513/1‑150 (MQ=255)
                                                                                                                                                  aCCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTc     <  2:5361927/150‑1 (MQ=255)
                                                                                                                                                  aCCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTc     >  1:1456012/1‑150 (MQ=255)
                                                                                                                                                   ccGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGCAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGACGTGTCg    >  1:4729965/1‑150 (MQ=255)
                                                                                                                                                    cGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTCGa   >  1:5448659/1‑150 (MQ=255)
                                                                                                                                                    cGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTCGa   >  2:1590689/1‑150 (MQ=255)
                                                                                                                                                     gACGCTTTGGCGGGTTAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTCGaa  <  1:3720992/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTCGAA  >  NC_002947/5541923‑5542221

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: