Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,962,581 T→A 5.5% intergenic (‑231/‑98) mnxG ← / → PP_3491 manganese‑oxidizing multicopper oxidase/enoly‑coenzyme A hydratase/isomerase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,962,5810TA 5.5% 1038.3 / 28.5 305intergenic (‑231/‑98)mnxG/PP_3491manganese‑oxidizing multicopper oxidase/enoly‑coenzyme A hydratase/isomerase family protein
Reads supporting (aligned to +/- strand):  ref base T (132/155);  new base A (10/7);  total (142/163)
Fisher's exact test for biased strand distribution p-value = 3.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 310 total aligned reads displayed.

ACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAGACCATGTACTGGCCGAGG  >  NC_002947/3962432‑3962730
                                                                                                                                                     |                                                                                                                                                     
aCTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGt                                                                                                                                                       <  1:2176814/150‑1 (MQ=255)
  ttGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCg                                                                                                                                                     >  2:2086444/1‑150 (MQ=255)
  ttGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCg                                                                                                                                                     >  1:6442828/1‑150 (MQ=255)
  ttGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCg                                                                                                                                                     >  1:7640140/1‑150 (MQ=255)
  ttGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAAGTGTCg                                                                                                                                                     >  1:2556450/1‑150 (MQ=255)
   tGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGa                                                                                                                                                    <  2:3270096/150‑1 (MQ=255)
       cGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATtgt                                                                                                                                                <  1:6414936/150‑1 (MQ=255)
       cGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATtgt                                                                                                                                                <  2:876753/150‑1 (MQ=255)
         aTCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGAt                                                                                                                                                   <  2:1852105/145‑1 (MQ=255)
            aaCCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAg                                                                                                                                           <  2:153063/150‑1 (MQ=255)
            aaCCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAg                                                                                                                                           <  2:5138184/150‑1 (MQ=255)
             aCCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAgg                                                                                                                                          <  1:2374213/150‑1 (MQ=255)
               cGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCa                                                                                                                                        >  2:8240083/1‑150 (MQ=255)
                gCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCaa                                                                                                                                            >  1:8123741/1‑145 (MQ=255)
                 cAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCaaa                                                                                                                                      <  1:8630005/150‑1 (MQ=255)
                 cAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCaaa                                                                                                                                      >  1:2826855/1‑150 (MQ=255)
                 cAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCaaa                                                                                                                                      >  1:6045142/1‑150 (MQ=255)
                       aTTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCa                                                                                                                                >  2:4097800/1‑150 (MQ=255)
                           gggAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCACATTCGCATGcc                                                                                                                            >  2:5052184/1‑150 (MQ=255)
                              aggaTCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAgg                                                                                                                         <  1:5927083/150‑1 (MQ=255)
                                     cccAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGc                                                                                                                  >  1:8594878/1‑150 (MQ=255)
                                     cccAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGc                                                                                                                  >  2:8900761/1‑150 (MQ=255)
                                      ccAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCg                                                                                                                 >  1:5533434/1‑150 (MQ=255)
                                            ggTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGAtcctc                                                                                                           <  2:7354467/150‑1 (MQ=255)
                                            ggTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGAtcctc                                                                                                           >  2:1717909/1‑150 (MQ=255)
                                              ttCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCa                                                                                                         <  2:3429546/150‑1 (MQ=255)
                                              ttCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGAGTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCa                                                                                                         <  2:7817914/150‑1 (MQ=255)
                                               tCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAg                                                                                                        <  1:1969933/150‑1 (MQ=255)
                                                cccccGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGc                                                                                                       <  1:4583291/150‑1 (MQ=255)
                                                cccccGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGc                                                                                                       <  1:5983323/150‑1 (MQ=255)
                                                cccccGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGc                                                                                                       <  2:6096673/150‑1 (MQ=255)
                                                  cccGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGa                                                                                                     >  1:6831873/1‑150 (MQ=255)
                                                   ccGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGaa                                                                                                    >  1:7861064/1‑150 (MQ=255)
                                                       aTCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTcc                                                                                                <  1:8900761/150‑1 (MQ=255)
                                                        tCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCa                                                                                               <  1:2189292/150‑1 (MQ=255)
                                                           ggggACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGAc                                                                                            >  2:6793014/1‑150 (MQ=255)
                                                             ggACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACtt                                                                                          >  1:4020024/1‑150 (MQ=255)
                                                              gACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTc                                                                                         <  1:4350034/150‑1 (MQ=255)
                                                                 ggAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACt                                                                                      <  1:3737758/150‑1 (MQ=255)
                                                                 ggAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACt                                                                                      >  2:8948030/1‑150 (MQ=255)
                                                                  gAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACtt                                                                                     <  2:8309936/150‑1 (MQ=255)
                                                                  gAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACtt                                                                                     <  1:5831908/150‑1 (MQ=255)
                                                                  gAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACtt                                                                                     <  1:972106/150‑1 (MQ=255)
                                                                   aaTCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTg                                                                                    >  2:2838233/1‑150 (MQ=255)
                                                                     tCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGa                                                                                  >  1:2338760/1‑150 (MQ=255)
                                                                       ccGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGaaa                                                                                >  2:5760032/1‑150 (MQ=255)
                                                                        cGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGaaaa                                                                               <  1:5046004/150‑1 (MQ=255)
                                                                        cGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGaaaa                                                                               <  1:3081935/150‑1 (MQ=255)
                                                                         gTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAAt                                                                              <  1:561084/150‑1 (MQ=255)
                                                                         gTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAAt                                                                              <  1:1034473/150‑1 (MQ=255)
                                                                         gTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAAt                                                                              <  2:2226628/150‑1 (MQ=255)
                                                                          tGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATc                                                                             <  1:5100813/150‑1 (MQ=255)
                                                                          tGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAAGTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATc                                                                             <  2:6609274/150‑1 (MQ=255)
                                                                            ggggCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCCCACAAGGTAAAACCGCACGGTACCCTTGCACAATTGTCGCTTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCaa                                                                           <  2:3184473/150‑1 (MQ=255)
                                                                            ggggCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCaa                                                                           <  2:7540061/150‑1 (MQ=255)
                                                                            ggggCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTAGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCaa                                                                           <  2:9048948/150‑1 (MQ=255)
                                                                             gggCAACGGCGGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAACCGCACGGTACCCTTGCACAAGTGTCGATTGTGCAAGGCAAAGTCGCATGCCAGGATCCGGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAAt                                                                          <  2:7474344/150‑1 (MQ=255)
                                                                             gggCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAAt                                                                          <  2:6397585/150‑1 (MQ=255)
                                                                             gggCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAAt                                                                          >  1:7795324/1‑150 (MQ=255)
                                                                             gggCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAAt                                                                          <  2:5439756/150‑1 (MQ=255)
                                                                               gCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATaa                                                                        >  1:5468659/1‑150 (MQ=255)
                                                                                   aGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACaaa                                                                    <  1:2228487/150‑1 (MQ=255)
                                                                                   aGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACaaa                                                                    <  2:4729863/150‑1 (MQ=255)
                                                                                      ctgctCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACCCACGCTACCCTTGCACAATTGACGATTGTGCAAGGCAAATTCGCATCCCACGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTCGAAAATCAATAACAAAGCa                                                                 <  2:6976688/150‑1 (MQ=255)
                                                                                       tgctCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAg                                                                <  2:6434088/150‑1 (MQ=255)
                                                                                         ctCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCTAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAggg                                                              <  2:2354349/150‑1 (MQ=255)
                                                                                          tCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGga                                                             <  1:1821461/150‑1 (MQ=255)
                                                                                          tCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGga                                                             <  2:2441177/150‑1 (MQ=255)
                                                                                            cacaAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGgaga                                                           >  2:7453524/1‑150 (MQ=255)
                                                                                              caAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAAc                                                         >  1:7978206/1‑150 (MQ=255)
                                                                                              caAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAAc                                                         <  1:3261195/150‑1 (MQ=255)
                                                                                                  cgcgcgcgGGCCTCACACACGCCAAAAACGCACGGTCCCCTTGCACAATTGTCGCTTGTGCAAGGCAAAGTCGCAGGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCaa                                                     <  2:6274030/150‑1 (MQ=255)
                                                                                                   gcgcgcTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAAt                                                    <  1:4227407/150‑1 (MQ=255)
                                                                                                   gcgcgcTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAAt                                                    <  2:3662295/150‑1 (MQ=255)
                                                                                                     gcgcTGGCCTCACACAAGGTAAAAACGCTCGTTACCCTTGCACAATTGTCGATTGTGCAAGTCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGa                                                  <  2:1198308/150‑1 (MQ=255)
                                                                                                         tGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGc                                              >  1:7341515/1‑150 (MQ=255)
                                                                                                             cTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCacac                                          >  2:8371011/1‑150 (MQ=255)
                                                                                                                  acaAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCa                                     >  2:8670574/1‑150 (MQ=255)
                                                                                                                  acaAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCa                                     >  2:4266293/1‑150 (MQ=255)
                                                                                                                  acaAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAAACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCa                                     >  1:4692521/1‑150 (MQ=255)
                                                                                                                     agctAAAAACGCACGGTACCCTTGCACAATTCTCGATTGAGCAAGGCAAATTCGCATGCCAGGATCCAGCGATCCTCCAGCCAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATAc                                  <  2:8508340/147‑1 (MQ=255)
                                                                                                                       ggaaaaaCGCACGGTACCCTTGCACAATTGTCGAGTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATAccc                                >  1:4922874/3‑150 (MQ=255)
                                                                                                                       gTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATAccc                                <  1:8290330/150‑1 (MQ=255)
                                                                                                                       gTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATAccc                                >  1:5410828/1‑150 (MQ=255)
                                                                                                                           aaaCGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGgc                            >  2:1232505/1‑150 (MQ=255)
                                                                                                                             aCGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGgcag                          >  1:6712623/1‑150 (MQ=255)
                                                                                                                                  cGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACa                     <  2:6015235/150‑1 (MQ=255)
                                                                                                                                   ggTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAg                    >  2:8100222/1‑150 (MQ=255)
                                                                                                                                    gTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAGa                   <  2:6770527/150‑1 (MQ=255)
                                                                                                                                    gTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAGa                   <  2:4028260/150‑1 (MQ=255)
                                                                                                                                      aCCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAGAcc                 <  2:2543176/150‑1 (MQ=255)
                                                                                                                                         cTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAGACCATg              <  2:3534629/150‑1 (MQ=255)
                                                                                                                                         cTTGCACAATCGACAATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAGACCATg              >  1:5622166/1‑150 (MQ=255)
                                                                                                                                          ttGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAGACCATGt             >  2:390018/1‑150 (MQ=255)
                                                                                                                                          ttGCACAATCGACAATTGTGCAAGACAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAGACCATGt             >  2:6947761/1‑150 (MQ=255)
                                                                                                                                              acaATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAGACCATGTACTg         >  2:4420346/1‑150 (MQ=255)
                                                                                                                                               caATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCAAGACTTCTCTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAGACCATGTACTgg        >  2:2943581/1‑150 (MQ=255)
                                                                                                                                                aaTTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAGACCATGTACTGGc       >  1:8344697/1‑150 (MQ=255)
                                                                                                                                                   tGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAGACCATGTACTGGCCGa    <  1:3015364/150‑1 (MQ=255)
                                                                                                                                                    gTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAGACCATGTACTGGCCGAg   <  2:5501313/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAGACCATGTACTGGCCGAGG  >  NC_002947/3962432‑3962730

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