Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,145,849 T→A 7.3% intergenic (+36/‑87) PP_3647 → / → PP_3648 oxidoreductase/carboxymuconolactone decarboxylase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,145,8490TA 7.3% 427.3 / 18.3 137intergenic (+36/‑87)PP_3647/PP_3648oxidoreductase/carboxymuconolactone decarboxylase family protein
Reads supporting (aligned to +/- strand):  ref base T (63/63);  new base A (4/6);  total (67/70)
Fisher's exact test for biased strand distribution p-value = 7.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

Only 100 of 148 total aligned reads displayed.

GGTGACCACGCGATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGAGGACTAC  >  NC_002947/4145700‑4145998
                                                                                                                                                     |                                                                                                                                                     
  tGACCACGCGATCATTGTCGGCGCGGTGAAGAGCCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCtt                                                                                                                                                       <  1:4942451/148‑1 (MQ=255)
    aCCACGCGATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGccc                                                                                                                                                   <  1:3989819/150‑1 (MQ=255)
       aCGCGATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCcgcg                                                                                                                                                >  1:4955219/1‑150 (MQ=255)
       aCGCGATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCcgcg                                                                                                                                                >  2:3257319/1‑150 (MQ=255)
         gcgATCATTGTCGGCGCGGTGAAGGACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCAAGCCCGCGaa                                                                                                                                              <  1:4582332/150‑1 (MQ=255)
         gcgATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGaa                                                                                                                                              >  2:3221925/1‑150 (MQ=255)
         gcgATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGaa                                                                                                                                              >  1:2021049/1‑150 (MQ=255)
         gcgATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGaa                                                                                                                                              <  1:1947967/150‑1 (MQ=255)
         gcgATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCAAGCCCGCGaa                                                                                                                                              <  1:4799200/150‑1 (MQ=255)
         gcgATCAGTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCAAGCCCGCGaa                                                                                                                                              <  1:4582363/150‑1 (MQ=255)
           gATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGa                                                                                                                                            >  2:1220347/1‑150 (MQ=255)
               aTTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGcc                                                                                                                                        <  2:2396477/150‑1 (MQ=255)
                ttGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGccc                                                                                                                                       <  2:792821/150‑1 (MQ=255)
                ttGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGccc                                                                                                                                       <  2:792711/150‑1 (MQ=255)
                ttGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGccc                                                                                                                                       <  2:1619297/150‑1 (MQ=255)
                   tCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTg                                                                                                                                    >  2:4770928/1‑150 (MQ=255)
                          ggTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAg                                                                                                                             <  2:2021049/150‑1 (MQ=255)
                             gaagaaCCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAgggg                                                                                                                          >  2:3022805/1‑150 (MQ=255)
                                gaaCCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAgg                                                                                                                            >  1:1563249/1‑145 (MQ=255)
                                         gtTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCg                                                                                                              >  1:4025523/1‑150 (MQ=255)
                                          ttGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGt                                                                                                             <  2:1259418/150‑1 (MQ=255)
                                          ttGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGt                                                                                                             >  1:646905/1‑150 (MQ=255)
                                          ttGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGt                                                                                                             <  2:2524044/150‑1 (MQ=255)
                                            gCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGtct                                                                                                           <  2:4534817/150‑1 (MQ=255)
                                            gCCTGAGCAGCCAATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGtct                                                                                                           >  2:4305599/1‑150 (MQ=255)
                                             ccTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGtctc                                                                                                          >  1:3694224/1‑150 (MQ=255)
                                             ccTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTGGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGtctc                                                                                                          <  2:1498397/150‑1 (MQ=255)
                                                 agcagcCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCAccc                                                                                                      >  1:4232127/1‑150 (MQ=255)
                                                   cagcCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTc                                                                                                    >  1:1442707/1‑150 (MQ=255)
                                                   cagcCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTc                                                                                                    <  2:4079468/150‑1 (MQ=255)
                                                     gcctttGACCCCGTTGGTGTACCACAAGGGCAAGCTGGGCGCCCTGCCGCCACTTGGCTTTTAAATCTTGAGTCGATTTTGGCATCTTCGGGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCgg                                                                                                  <  1:715995/145‑1 (MQ=255)
                                                      ccGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGt                                                                                                 >  1:168371/1‑150 (MQ=255)
                                                        gATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTGTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTg                                                                                               >  2:3991094/1‑150 (MQ=255)
                                                              cccGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGccc                                                                                         >  1:2236724/1‑150 (MQ=255)
                                                              cccGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGccc                                                                                         <  1:5288953/150‑1 (MQ=255)
                                                              cccGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCGGATACAGCTTGAGGCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTGTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGccc                                                                                         <  2:4557123/150‑1 (MQ=255)
                                                                cGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGtc                                                                                                            <  1:3927806/129‑1 (MQ=255)
                                                                cGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCt                                                                                       <  2:2236724/150‑1 (MQ=255)
                                                                    ggTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCttttt                                                                                   <  2:1073856/150‑1 (MQ=255)
                                                                    ggTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCttttt                                                                                   <  1:5346678/150‑1 (MQ=255)
                                                                      tgtACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCg                                                                                             <  1:3874502/138‑1 (MQ=255)
                                                                      tgtACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCg                                                                                             >  2:3874502/1‑138 (MQ=255)
                                                                        tACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATg                                                                               >  1:1536867/1‑150 (MQ=255)
                                                                         aCCACACGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGc                                                                              <  1:2479103/150‑1 (MQ=255)
                                                                         aCCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGc                                                                              <  1:3465841/150‑1 (MQ=255)
                                                                         aCCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGc                                                                              <  1:2929032/150‑1 (MQ=255)
                                                                           cacaAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCtt                                                                                                                                                       <  2:5354743/75‑1 (MQ=255)
                                                                           cacaAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCtt                                                                                                                                                       >  1:5354743/1‑75 (MQ=255)
                                                                           cacaAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCttt                                                                                     >  1:318502/1‑141 (MQ=255)
                                                                           cacaAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCg                                                                            <  2:1935979/150‑1 (MQ=255)
                                                                             caAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGa                                                                          <  2:2534767/150‑1 (MQ=255)
                                                                              aaGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAg                                                                         <  1:4693813/150‑1 (MQ=255)
                                                                                 ggCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTAc                                                                      <  1:1293011/150‑1 (MQ=255)
                                                                                   cAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTg                                                                    >  1:2897569/1‑150 (MQ=255)
                                                                                    aaGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGc                                                                   <  1:2243373/150‑1 (MQ=255)
                                                                                      gCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGcac                                                                 >  1:1224625/1‑150 (MQ=255)
                                                                                       cTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGcaca                                                                >  2:3542612/1‑150 (MQ=255)
                                                                                          ggCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGa                                                             >  1:4258403/1‑150 (MQ=255)
                                                                                          ggCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGa                                                             <  2:1280541/150‑1 (MQ=255)
                                                                                          ggCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGa                                                             <  2:1358145/150‑1 (MQ=255)
                                                                                          ggCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGa                                                             <  1:1892692/150‑1 (MQ=255)
                                                                                           gcgcCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAg                                                            <  1:4322442/150‑1 (MQ=255)
                                                                                              cccTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGTGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCaa                                                         >  1:1438188/1‑150 (MQ=255)
                                                                                                cTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCCATg                                                       >  2:4695537/1‑150 (MQ=255)
                                                                                                  gccgccACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGaa                                                     <  1:3924390/150‑1 (MQ=255)
                                                                                                   ccgccACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGaaa                                                    >  2:2472223/1‑150 (MQ=255)
                                                                                                     gccACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGa                                                                                                                                               <  2:5349809/57‑1 (MQ=255)
                                                                                                        aCTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTAc                                               >  2:1609803/1‑150 (MQ=255)
                                                                                                        aCTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTAc                                               >  2:2681309/1‑150 (MQ=255)
                                                                                                        aCTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTAc                                               <  1:3039735/150‑1 (MQ=255)
                                                                                                          ttGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGa                                             >  2:3226882/1‑150 (MQ=255)
                                                                                                           tGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGaa                                            <  2:4955219/150‑1 (MQ=255)
                                                                                                                gATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACt                                                                     <  2:3092547/120‑1 (MQ=255)
                                                                                                                 aTACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGc                                      <  2:2650130/150‑1 (MQ=255)
                                                                                                                  tACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGcc                                     <  1:3188887/150‑1 (MQ=255)
                                                                                                                     aTCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATc                                                                                                               >  1:3674267/1‑73 (MQ=255)
                                                                                                                     aTCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATc                                                                                                               <  2:3674267/73‑1 (MQ=255)
                                                                                                                     aTCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGa                                  >  1:1225484/1‑150 (MQ=255)
                                                                                                                       cTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAgg                                                                                                                            >  2:3889196/1‑58 (MQ=255)
                                                                                                                           aGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATcc                            <  1:1609803/150‑1 (MQ=255)
                                                                                                                           aGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATcc                            <  2:397975/150‑1 (MQ=255)
                                                                                                                               gATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAACAGGGCCTGTAAATCCGCAc                        >  2:4862211/1‑150 (MQ=255)
                                                                                                                               gATGTTGGCATCTTAGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTAGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCAGCaa                        >  1:3759792/1‑149 (MQ=255)
                                                                                                                                aTGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCAcc                       >  2:2923673/1‑150 (MQ=255)
                                                                                                                                      gCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTa                 <  1:4087476/150‑1 (MQ=255)
                                                                                                                                      gCATCTTCGCGGGCAAGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTa                 >  2:3079861/1‑150 (MQ=255)
                                                                                                                                       cATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTAc                >  1:341418/1‑150 (MQ=255)
                                                                                                                                       cATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGAGCCTGAAAATCCGCACCCAGGTAc                >  2:5131138/1‑150 (MQ=255)
                                                                                                                                        aTCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACt               >  2:4693683/1‑150 (MQ=255)
                                                                                                                                          cTTCGCGGGCAAGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAAt                                                        <  1:248094/107‑1 (MQ=255)
                                                                                                                                          cTTCGCGGGCAAGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAAt                                                        >  2:248094/1‑107 (MQ=255)
                                                                                                                                              gcgGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGa         <  1:1818852/150‑1 (MQ=255)
                                                                                                                                              gcgGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGa         <  1:3069835/150‑1 (MQ=255)
                                                                                                                                              gcgGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCACTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGa         <  1:2719169/150‑1 (MQ=255)
                                                                                                                                               cgGGCAAGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGAg        <  1:3079861/150‑1 (MQ=255)
                                                                                                                                                  gCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGAGGAc     <  1:3221925/150‑1 (MQ=255)
                                                                                                                                                    ttGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGAGGACTa   <  2:4258403/150‑1 (MQ=255)
                                                                                                                                                     tGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTTGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGAGGACTAc  <  1:3022805/150‑1 (MQ=255)
                                                                                                                                                     tGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGAGGACTAc  <  2:4232127/150‑1 (MQ=255)
                                                                                                                                                     tGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGAGGACTAc  <  1:2749516/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGTGACCACGCGATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGAGGACTAC  >  NC_002947/4145700‑4145998

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