Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,539,045 G→C 13.6% 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,0450GC13.6% 455.5 / 38.8 175intergenic (+52/+51)nadE/PP_4870NH(3)‑dependent NAD(+) synthetase/azurin
Reads supporting (aligned to +/- strand):  ref base G (66/83);  new base C (13/12);  total (79/96)
Fisher's exact test for biased strand distribution p-value = 5.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 177 total aligned reads displayed.

CTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTCGAA  >  NC_002947/5538898‑5539194
                                                                                                                                                   |                                                                                                                                                     
 tGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGCAGACGCGGGTGAACCCACCGACg                                                                                                                                                    >  2:5447339/1‑150 (MQ=255)
  gCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACCTCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGc                                                                                                                                                   <  2:547816/150‑1 (MQ=255)
  gCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGc                                                                                                                                                   <  1:6682951/150‑1 (MQ=255)
        gCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGCAAACGAGGGTGAACCCACCCACGCCTTcgc                                                                                                                                             >  2:7941206/1‑150 (MQ=255)
          aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                           <  2:4864544/150‑1 (MQ=25)
          aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                           <  2:5903046/150‑1 (MQ=25)
          aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                           <  1:4752405/150‑1 (MQ=25)
          aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                           <  1:7589429/150‑1 (MQ=25)
          aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                           <  2:4351501/150‑1 (MQ=25)
          aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                           <  2:8769025/150‑1 (MQ=25)
          aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                           <  2:1391026/150‑1 (MQ=25)
          aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                           <  1:712816/150‑1 (MQ=255)
          aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                           <  1:712770/150‑1 (MQ=255)
          aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                           <  1:4675078/150‑1 (MQ=255)
          aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                           <  1:8682035/150‑1 (MQ=255)
          aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                           <  2:3312285/150‑1 (MQ=255)
          aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                           <  2:2237217/150‑1 (MQ=255)
          aGCCGGTGGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                           <  1:3487603/150‑1 (MQ=255)
           gCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCggg                                                                                                                                          <  1:1185038/150‑1 (MQ=255)
           gCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTAAACCGCGAATACGAGGGTGAACCCACCGACGCATTCGCggg                                                                                                                                          >  2:3539560/1‑150 (MQ=255)
             cGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTg                                                                                                                                        <  2:7387714/150‑1 (MQ=25)
             cGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTGACCCCCCCAAAAGAGGGTTGAACCACCCCCCCAATCGCggggg                                                                                                                                        >  1:2746584/1‑148 (MQ=255)
                ttGACCAGGCAGCGTTCGAGATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAAc                                                                                                                                     <  2:3154669/150‑1 (MQ=255)
                    ccAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGc                                                                                                                                 <  2:189846/150‑1 (MQ=255)
                     cAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCCCGAATACCAGGGTGAACCCACCGACGCATTCGCGGGTGAAACCGcg                                                                                                                                >  2:8441706/1‑149 (MQ=255)
                        gCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCCCCAATAAGAGGGTGTAACCACCCACCCCATAGCGGGTGGACACCCTccg                                                                                                                             >  1:692254/1‑149 (MQ=255)
                        gCAGCGTTCGACATCAGCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCGTTGATAAGTGGCCTGCAGGAGCGGGTACACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCt                                                                                                                             <  2:8420325/150‑1 (MQ=255)
                            cGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTAACCCGCCCAAAAGAGGGTGAACCCACCCACGCATTCGCGGGTGAACCCCCGCcagcag                                                                                                                         >  1:6708695/1‑150 (MQ=255)
                             gTTCGACATCATCGTCGCCACCTCCCGCAAGACCCCGCCCACGCGCGACCTGCCGTTCGCTCCTTGATACGTGGCCTGCAGGAGGGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGACCCCGCTCCTGCAgg                                                                                                                        <  2:6584078/150‑1 (MQ=255)
                             gTTCGACATCATCGTCGCCACCTCCCGCAAGACCCCGCCAAGGGGCGAAGGGCCGTTCGGTCCTTTATAAGGGGCCTGCAGGCGGGGTTTCACCCGCGAATACGTGGTTGAACCCACCGACGATTTCGGGGTTGAACCCGCTCCTGCAgg                                                                                                                        <  2:1843483/150‑1 (MQ=255)
                             gTTCGACATCATCGTCGCCACCTACCGCAACCCCCAGCACACGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGTTTCACCGGCTAATACGTGGTTGAACCCACCGACGAATGCGCGGGTGAACCCGCTCCTGCAgg                                                                                                                        <  2:6442561/150‑1 (MQ=255)
                              tgcGACATCTTCTTCGCCACCTCCCGCACGACCCAGCACACGCGCGAACTGCCGTTCGCTCCGTTATAAGTGGCCTGCAGGCGCGGGGTCACCCGCGAACCCGCGGGTGAACCCACCGACGCATTCGCGGGCGAACCCGCTCCTGCAggg                                                                                                                       <  2:5670546/148‑1 (MQ=255)
                                       atcGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCCAATAAGAGGGTGAAACCACCCACGCAGTCGCGGGTGAACCCGCGCCGGCAGGGTTCGCCggc                                                                                                              >  2:890507/1‑148 (MQ=255)
                                             gCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATACGTGGCCTGCAGGCGGGGTTTCACCCGCTAATACGAGGGTGAACCCACCGACGCATTCGCGGGTTAACCCGCTCCTGCAGGGTTCGCCGTCAGCCtt                                                                                                        <  2:3590984/150‑1 (MQ=255)
                                                cccTCCCGCAACACCCCGCACACGCGCGAACTGCCGTTCGCTCCTTGATACGGGGCCTGCAGGAGCGGTTTCACCCGCGAATACGCGGGTGAACCCACCGACGCATTCGGGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGa                                                                                                     <  2:3470970/149‑1 (MQ=38)
                                                  cacccGCAAGCCCCCGCACACGCGCCAACTGCCGTTCGCTCCTTGATAAGGGGCCTGCAGGCGCGGGTTCACCGGCTAATACGTGGGTTAACCCACCGACGAATTGGCGGGTTAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGAtt                                                                                                   <  1:6595045/147‑1 (MQ=255)
                                                  cTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGTGCGGGGTCACCCGCGCAAACGCGGGTGTACCCACCGACGCCATCGCGGGGGAACCCGCGCCGGGAGGGTTCGCCGTAAGCCATTGAtt                                                                                                   >  2:2668481/1‑150 (MQ=255)
                                                  cTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGTACCCACCGACGCAGTCGCGGGTGAACCCCCGCCTGCCGGGTTCGCCGTCAGCCGTTGAtt                                                                                                   >  1:7531344/1‑150 (MQ=255)
                                                     ccGCAAGACCCAGCAAACGCGCTAACTGCCGTGCGCTCCTTTATATGTGGCCGGCTGGAGCGGGTTCACCCGCGAATACGCGGGTGAACCCACCGACGCTTTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACt                                                                                                <  1:6019357/150‑1 (MQ=38)
                                                      cGCAAGACCCCGCACAAGCGCGAACTGCCGTTCGCTCCTTGATACGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACtt                                                                                               <  2:7422247/150‑1 (MQ=255)
                                                        cAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGa                                                                                             <  2:6803746/150‑1 (MQ=255)
                                                        cAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATAAGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGGTTACGTGa                                                                                             >  2:6816512/1‑150 (MQ=255)
                                                            aCCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTGGGCCCGCGAATAAGAGGGGGTAACCACCCACCCCATCGCGGGTGAACCCCCGCCGGCAGGGGTGGCCGGCAGCCGTTGGTTACTTGGCGAt                                                                                         >  1:7563725/1‑150 (MQ=255)
                                                              cccgCACACGCGCTAACTGCCGTGCGCTCCTTGATATGTGGCCTGCTGGAGCGGGTTCACCCGCTAATAGGTGGGTGAACCCACCGACGCTTTCGCGGTTTAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGa                                                                                       <  1:1316083/147‑1 (MQ=255)
                                                                aGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGCATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCCGGGTTCGCCGTCAGCCTTTGATTACGTGACGATGAcc                                                                                     >  2:6069496/1‑150 (MQ=255)
                                                                     aggagcgAACGGCAGTTCGCGCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTTACGATGACCTTGcc                                                                                >  2:5949331/5‑150 (MQ=38)
                                                                      cgcgcGAACTGCCGTTCGCTCCTGGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCt                                                                               <  2:524293/149‑1 (MQ=255)
                                                                              cTGCCGTGCGCTCCTTTATACGGGGCCTGCAGGAGCGGGTACACCCGCGATTACGAGGGTGAACCCACCGCCGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTcatcat                                                                       <  2:3849508/150‑1 (MQ=38)
                                                                                   gTTCGCTCCTTGATATGTGGCCTGCAGGTGGGGGTTCACCCGCGAATACGAGGGTTAACCCACCGCCGCTTTCGGGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCgaga                                                                  <  2:948958/150‑1 (MQ=255)
                                                                                    ttCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAAAAAGAGGGTGAACCCCCCGACGCATTAGCGGGTGAACCCGCGACTGCAGGGTTCGCCGTCCGCCTTTGAGTACTTGACGGTGACCATTCCTTTCATCATCgagct                                                                 >  2:5887024/1‑148 (MQ=255)
                                                                                    ttCGCTCCTTGATAAGTGCCCTGCAGGAGCGTGTACACCCGCGAATACGTGGGTCAACCCACCGACGAATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCATTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGAt                                                                 <  2:1906342/150‑1 (MQ=38)
                                                                                      cGCCCCTTTATATGTGGCCTGCAGGCGGGGGTACACCCGCTAATACTCGGGTGACCCCACCGCCGCTTGCGGGGGGCAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGAtgt                                                               <  1:5540042/150‑1 (MQ=255)
                                                                                                 atggggCCTTCAGGCGGGGGTCCACCCGCTAACCCGTGGGTGACCCCCCCGCCGCTTTGGGGGGTTCCCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAg                                                    <  1:1913134/146‑1 (MQ=255)
                                                                                                  cgTGGCCTGCTGGAGCGGGTTCACCCGCGAATACGCGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAgg                                                   <  2:3983984/149‑1 (MQ=255)
                                                                                                       ccTGCGGGTGGGGTTGCACCGGCGAATACGCGGGTCACCCCACCGACGCATTGGGGGGTCAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCag                                              <  1:5017223/150‑1 (MQ=255)
                                                                                                         ggccggcgCGGGGTCCCCCGCGAAGCCGCGGGTGAACCCACCGCCGCATGCGCGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCagaa                                            <  2:7296033/143‑1 (MQ=255)
                                                                                                          gCAGGAGCGGGTTCACCCGCGACTACGAGGGTGAACCCACCTACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGTTGTGGCCCGTGAAGGTGCagaag                                           <  2:8208449/150‑1 (MQ=255)
                                                                                                              gAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAACCCGCTCCTGCGGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGc                                       >  2:7070096/1‑150 (MQ=21)
                                                                                                               aGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCt                                      <  2:3691294/150‑1 (MQ=255)
                                                                                                                gCGGGTTCACCCGCGAATACGAGGTTGAACCCACCGACGAATTCGCGGGTGAACACGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACATTGCCTTTCATCATCGAGATTTGGCCAGGGAAGGTGCAGAAGAAGCTg                                     >  1:739543/1‑150 (MQ=255)
                                                                                                                  gggTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTa                                   <  2:4347650/150‑1 (MQ=255)
                                                                                                                   ggTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTGGCCATTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAg                                  <  2:7538176/150‑1 (MQ=255)
                                                                                                                    gTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGt                                 >  2:5523497/1‑150 (MQ=255)
                                                                                                                    gTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGt                                 >  1:6736379/1‑150 (MQ=255)
                                                                                                                        aCCCGCGAATACGAGGGTGACCCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcacc                             <  2:2884844/150‑1 (MQ=255)
                                                                                                                        aCCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcacc                             >  1:4817453/1‑150 (MQ=255)
                                                                                                                         cccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcacca                            >  1:4320331/1‑150 (MQ=255)
                                                                                                                          ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                           >  2:5095699/1‑150 (MQ=255)
                                                                                                                          ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                           >  1:2860671/1‑150 (MQ=255)
                                                                                                                          ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                           >  2:2488770/1‑150 (MQ=255)
                                                                                                                          ccGCGAATACGAGGGTGAACCCACCGACGCAGTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                           >  2:844908/1‑150 (MQ=255)
                                                                                                                          ccGCGAAGACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCCTCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                           >  2:4568590/1‑150 (MQ=255)
                                                                                                                           cgcgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCaccacc                          >  2:806702/1‑150 (MQ=255)
                                                                                                                            gcgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCg                         >  1:7708810/1‑150 (MQ=255)
                                                                                                                            gcgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCg                         >  2:611753/1‑150 (MQ=255)
                                                                                                                              gAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGc                       >  1:5183869/1‑150 (MQ=255)
                                                                                                                                 tACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTc                    >  1:1934815/1‑150 (MQ=255)
                                                                                                                                 tACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTc                    >  1:2422468/1‑150 (MQ=255)
                                                                                                                                 caggcgggTCAACCCACCGCGGCTTCCGCGGGGCACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTc                    <  1:3150322/145‑1 (MQ=255)
                                                                                                                                    gAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAg                 >  2:5754963/1‑150 (MQ=255)
                                                                                                                                       ggTGAACCCACCGACGCATTCGCGGTTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCtt              <  2:4301470/150‑1 (MQ=255)
                                                                                                                                       ggTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGACGCGGTAGTCACCACCGGCTTCAAGCtt              >  2:1586943/1‑150 (MQ=255)
                                                                                                                                        ctcaaCCCACCCACTAATGCGCGGGACAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTg             <  1:2503261/147‑1 (MQ=255)
                                                                                                                                          gAACCCACCGACGCTTTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGa           <  2:1377215/150‑1 (MQ=255)
                                                                                                                                           accccACCGACGCATTCGGGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGaa          <  1:8883665/148‑1 (MQ=255)
                                                                                                                                           accccACCGACGCATGCGGGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGaa          <  1:8441706/148‑1 (MQ=255)
                                                                                                                                           aaCCCACCGACTCAGTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGaa          <  2:175837/150‑1 (MQ=255)
                                                                                                                                            aCCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACtt                                                                                               >  2:3421884/1‑64 (MQ=255)
                                                                                                                                            aCCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAg         >  2:213394/1‑150 (MQ=255)
                                                                                                                                            aCCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAg         >  2:4417787/1‑150 (MQ=255)
                                                                                                                                            aCCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGGTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAg         >  1:8184342/1‑150 (MQ=255)
                                                                                                                                              ccACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAgtg       >  1:9057034/1‑150 (MQ=255)
                                                                                                                                               cACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAgtgt      >  1:2248629/1‑150 (MQ=255)
                                                                                                                                               cACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAgtgt      >  1:2248482/1‑150 (MQ=255)
                                                                                                                                                 ccGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCa                  >  2:5663507/1‑136 (MQ=255)
                                                                                                                                                 ccGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCa                  <  1:5663507/136‑1 (MQ=255)
                                                                                                                                                  cGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTCGa   >  2:771954/1‑150 (MQ=255)
                                                                                                                                                  cGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTCGa   >  2:4644256/1‑150 (MQ=255)
                                                                                                                                                  cGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTCGa   >  2:3696200/1‑150 (MQ=255)
                                                                                                                                                   gACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTCGaa  <  2:7558496/150‑1 (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: