Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,960,565 A→T 10.0% intergenic (‑137/‑192) 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,960,5650AT10.0% 196.5 / 13.0 70intergenic (‑137/‑192)mnxG/PP_3491manganese‑oxidizing multicopper oxidase/enoly‑coenzyme A hydratase/isomerase family protein
Reads supporting (aligned to +/- strand):  ref base A (31/32);  new base T (4/3);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.47e-01

GCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGA  >  NC_002947/3960419‑3960709
                                                                                                                                                  |                                                                                                                                                
gcgTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGATTCg                                                                                                                                               >  1:1003358/1‑150 (MQ=255)
gcgTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCg                                                                                                                                               >  1:53794/1‑150 (MQ=255)
gcgTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCg                                                                                                                                               >  1:838761/1‑150 (MQ=255)
      aCATGCTGCACTTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACTATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGATTCGGGGGaa                                                                                                                                         <  2:1003358/150‑2 (MQ=255)
       cATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACg                                                                                                                                        <  1:286328/150‑1 (MQ=255)
         tgctgcACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGa                                                                                                                                      >  2:303505/1‑150 (MQ=255)
         tgctgcACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGa                                                                                                                                      >  2:1314500/1‑150 (MQ=255)
          gctgcACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGaa                                                                                                                                     <  1:1314500/150‑1 (MQ=255)
             gcACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATcc                                                                                                                                  >  1:632303/1‑150 (MQ=255)
              cACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATccc                                                                                                                                 <  2:1770775/150‑1 (MQ=255)
                  tCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTgg                                                                                                                             <  1:1344485/150‑1 (MQ=255)
                  tCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTgg                                                                                                                             <  1:139401/150‑1 (MQ=255)
                    cGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTgggg                                                                                                                           <  2:632303/150‑1 (MQ=255)
                     gggACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTggggg                                                                                                                          <  2:1397412/150‑1 (MQ=255)
                         cacaGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGccac                                                                                                                      >  2:1779640/1‑147 (MQ=255)
                          acaGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAg                                                                                                                     >  1:999544/1‑150 (MQ=255)
                          acaGGATCTCTGCGGGGAAAGGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTACTCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAg                                                                                                                     <  1:1532818/150‑1 (MQ=255)
                            aGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGgc                                                                                                                   >  1:229160/1‑150 (MQ=255)
                              gATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGATTACCCCGAATTGGGGGTCGGAATCCCGTGGGGGCAAAGgctg                                                                                                                 >  2:1589287/1‑150 (MQ=255)
                                     gCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGTACGGAATCCCGTGGGGGCAAAGGCTGCTCcaca                                                                                                          <  1:1589287/150‑1 (MQ=255)
                                        gggAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAgc                                                                                                       <  2:806049/150‑1 (MQ=255)
                                        gggAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAACCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAgc                                                                                                       <  2:1381661/150‑1 (MQ=255)
                                               ggAGCAAGGTGCGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTCCGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCACCCGAATCGGGGGACGGAATCCCGTGGCGGCAAAGGCTGCTCCACAAGCGCGCGCt                                                                                                <  2:1000551/150‑1 (MQ=255)
                                                gAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTg                                                                                               <  1:622475/150‑1 (MQ=255)
                                                   cAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGcc                                                                                            >  1:1224054/1‑150 (MQ=255)
                                                          tgtgCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGATGGCCTcacaca                                                                                     >  1:1238860/1‑150 (MQ=255)
                                                            tgCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAg                                                                                   <  2:1109853/150‑1 (MQ=255)
                                                            tgCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAg                                                                                   <  2:999544/150‑1 (MQ=255)
                                                            tgCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAg                                                                                   >  1:408532/1‑150 (MQ=255)
                                                             gCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAgg                                                                                  <  1:1189773/150‑1 (MQ=255)
                                                              ccAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGt                                                                                 <  2:1570421/150‑1 (MQ=255)
                                                                agaagaAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTaa                                                                               <  2:838761/150‑1 (MQ=255)
                                                                          ttCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAg                                                                                                                     >  2:375324/1‑102 (MQ=255)
                                                                          ttCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAg                                                                                                                     >  2:375361/1‑102 (MQ=255)
                                                                          ttCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAg                                                                                                                     <  1:375361/102‑1 (MQ=255)
                                                                          ttCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAg                                                                                                                     <  1:375324/102‑1 (MQ=255)
                                                                             aaTATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACCGTAc                                                                  <  2:1462199/150‑1 (MQ=255)
                                                                                    aTCCATAACTTGTACGAGCAACCGGAACTTAGTGAGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGca                                                           <  2:408532/150‑1 (MQ=255)
                                                                                     tCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGcac                                                          >  2:1106186/1‑150 (MQ=255)
                                                                                     tCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGcac                                                          >  1:1384763/1‑150 (MQ=255)
                                                                                       aaTAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGTTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCtt                                                              <  2:247004/144‑1 (MQ=255)
                                                                                       aaTAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGTTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCtt                                                              >  1:247004/1‑144 (MQ=255)
                                                                                       aaTAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTATAAACGCACGGTACCCTTGCACaa                                                        <  2:945580/150‑1 (MQ=255)
                                                                                         tAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAAtt                                                      >  1:1668082/1‑150 (MQ=255)
                                                                                         tAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAAtt                                                      <  1:1618507/150‑1 (MQ=255)
                                                                                             ttGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTAccc                                                                >  2:180176/1‑136 (MQ=255)
                                                                                             ttGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTAccc                                                                <  1:180176/136‑1 (MQ=255)
                                                                                              tGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAATAACGCACGGTACCCTTGCACAATTGTCGa                                                 >  2:28915/1‑150 (MQ=255)
                                                                                              tGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGa                                                 <  2:53794/150‑1 (MQ=255)
                                                                                               gTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGATTCGGGGGaa                                                                                                                                         <  1:491203/61‑2 (MQ=255)
                                                                                               gTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGATTCGGGGGaa                                                                                                                                         >  2:491203/1‑60 (MQ=255)
                                                                                                 aCGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATtg                                              >  2:162936/1‑150 (MQ=255)
                                                                                                      cAACCGCAACTGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACCCATGGTAACCTTGCACAATTGTCGATTGTGCaa                                         <  1:250591/150‑1 (MQ=255)
                                                                                                       aaCCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAg                                        <  2:229160/150‑1 (MQ=255)
                                                                                                        aCCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCACGAATTGGGGTACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAgg                                       >  1:1774550/1‑150 (MQ=255)
                                                                                                               cGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTc                                <  1:303505/150‑1 (MQ=255)
                                                                                                                         aggaTCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAgg                      >  2:1085939/1‑150 (MQ=255)
                                                                                                                                cccAGCAGGTTCCCCCGAATCGGGGGACTGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGc               >  1:730831/1‑150 (MQ=255)
                                                                                                                                cccAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGc               <  2:1224054/150‑1 (MQ=255)
                                                                                                                                cccAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAAATTATTATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGc               >  1:99069/1‑150 (MQ=255)
                                                                                                                                 ccAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCg              >  1:744211/1‑150 (MQ=255)
                                                                                                                                  cagcagGTTCCCCCGAATCGGGGGACGGAATCACGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGa             >  2:1071053/1‑150 (MQ=255)
                                                                                                                                   agcagGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCa                             <  2:42665/133‑1 (MQ=255)
                                                                                                                                   agcagGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCa                             >  1:42665/1‑133 (MQ=255)
                                                                                                                                   agcagGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGAt            >  1:606697/1‑150 (MQ=255)
                                                                                                                                         gtCCCCCGATTCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCa      >  2:1376297/2‑150 (MQ=255)
                                                                                                                                          tCCCCCGATTCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATc                   >  2:343432/1‑136 (MQ=255)
                                                                                                                                          tCCCCCGATTCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATc                   <  1:343432/136‑1 (MQ=255)
                                                                                                                                             cccGGATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGa  >  1:1651819/1‑150 (MQ=255)
                                                                                                                                             cccGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGa  >  1:1150241/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGA  >  NC_002947/3960419‑3960709

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