Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,337,242 C→G 6.1% intergenic (+28/‑21) dksA → / → gluQ RNA polymerase‑binding transcription factor DksA/glutamyl‑Q tRNA(Asp) synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,337,2420CG 6.1% 1208.8 / 48.3 356intergenic (+28/‑21)dksA/gluQRNA polymerase‑binding transcription factor DksA/glutamyl‑Q tRNA(Asp) synthetase
Reads supporting (aligned to +/- strand):  ref base C (178/155);  new base G (14/8);  total (192/164)
Fisher's exact test for biased strand distribution p-value = 3.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 369 total aligned reads displayed.

GGTGCGAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGCCCGCGCCGTCAACGGC  >  NC_002947/5337093‑5337391
                                                                                                                                                     |                                                                                                                                                     
    cGAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAATAAGGAAAAATCGGTCGGGAATGGGTGTTCTTTCTTCCACCTGTTCGGGGCGGATCATgg                                                                                                                                                   >  2:6186493/1‑149 (MQ=255)
    cGAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATgc                                                                                                                                                   <  2:6683213/150‑1 (MQ=255)
     gAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATgcg                                                                                                                                                  >  2:6740850/1‑150 (MQ=255)
      aGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAATCGGTCGGCAAAGGCTGATCTTTCTTCCACCTGATCGGGGCGCATCATgcgc                                                                                                                                                 >  2:5823073/1‑150 (MQ=255)
      aGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATgcgc                                                                                                                                                 <  1:3510817/150‑1 (MQ=255)
         cGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGcccc                                                                                                                                              <  1:1117475/150‑1 (MQ=255)
          gTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATGATGCGCCCCg                                                                                                                                             <  1:4839857/150‑1 (MQ=255)
          gTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCg                                                                                                                                             <  1:236579/150‑1 (MQ=255)
           tGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGt                                                                                                                                            <  2:4673749/150‑1 (MQ=255)
            gCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGtt                                                                                                                                           <  1:6118585/150‑1 (MQ=255)
                cATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCAt                                                                                                                                       >  1:4458148/1‑150 (MQ=255)
                 aTCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCAtt                                                                                                                                      <  1:2062271/150‑1 (MQ=255)
                 aTCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCAtt                                                                                                                                      <  2:2686324/150‑1 (MQ=255)
                 aTCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCAtt                                                                                                                                      >  2:176918/1‑150 (MQ=255)
                     agaTCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTtata                                                                                                                                  >  2:326739/1‑150 (MQ=255)
                        tCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAATAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGc                                                                                                                               <  2:3755381/150‑1 (MQ=255)
                         cGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGcc                                                                                                                              <  2:1431660/150‑1 (MQ=255)
                           gCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCa                                                                                                                            <  1:208458/150‑1 (MQ=255)
                                  ccgccTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGAc                                                                                                                     >  1:4708865/1‑150 (MQ=255)
                                  ccgccTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGAc                                                                                                                     <  2:6344600/150‑1 (MQ=255)
                                  ccgccTGGAAGCCCGCCCTACCGCAGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGAGCCCCGTTTCATTTATATGACCAGCCTGAc                                                                                                                     >  1:6744294/1‑150 (MQ=255)
                                   cgccTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGAcc                                                                                                                    <  2:2841448/150‑1 (MQ=255)
                                    gccTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCa                                                                                                                   <  2:404157/150‑1 (MQ=255)
                                       tGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGa                                                                                                                >  1:3748127/1‑150 (MQ=255)
                                       tGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGa                                                                                                                <  1:6876437/150‑1 (MQ=255)
                                        ggAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGAc                                                                                                               <  2:4146806/150‑1 (MQ=255)
                                           aGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGa                                                                                                            >  2:1114879/1‑150 (MQ=255)
                                           aGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGa                                                                                                            <  1:7978/150‑1 (MQ=255)
                                                gcccTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCa                                                                                                       <  1:3832922/150‑1 (MQ=255)
                                                    tACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTa                                                                                                   <  1:1810755/150‑1 (MQ=255)
                                                      ccgccgACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGTCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCAtt                                                                                                                                      >  1:1512230/1‑113 (MQ=255)
                                                       cgccgACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACAt                                                                                                >  2:5511349/1‑150 (MQ=255)
                                                        gccgACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATc                                                                                               <  1:3066075/150‑1 (MQ=255)
                                                          cgACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCgg                                                                                             <  2:1941496/150‑1 (MQ=255)
                                                              cTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGgcgc                                                                                         <  2:2088810/150‑1 (MQ=255)
                                                              cTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGgcgc                                                                                         >  2:5964557/1‑150 (MQ=255)
                                                              cTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGgcgc                                                                                         >  1:2792261/1‑150 (MQ=255)
                                                               tgtgCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCt                                                                                        >  1:3438425/1‑150 (MQ=255)
                                                                gtgCTTCGACTGCAAGGAAATTGCCGATAAAAAGGAAAAAACTGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCAGTTTCTTTTATATGCCCAGCCTGACAATGACCGTCTCCAGCTACATCGGGCGCtt                                                                                       <  2:6119194/150‑1 (MQ=255)
                                                                   gtcaggcTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGc                                                                                    >  2:2120979/7‑150 (MQ=255)
                                                                    ttCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGcc                                                                                   >  2:2480205/1‑150 (MQ=255)
                                                                     tCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGccc                                                                                  >  1:4141332/1‑150 (MQ=255)
                                                                      cGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGcccc                                                                                 <  2:2629559/150‑1 (MQ=255)
                                                                          tGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGcccccacc                                                                             <  2:4007321/150‑1 (MQ=255)
                                                                             aaGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGGCCGACTCCAGCTACATCGGGCGCTTCGcccccaccccc                                                                          >  2:1381029/1‑150 (MQ=255)
                                                                             aaGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGcccccaccccc                                                                          >  1:2609549/1‑150 (MQ=255)
                                                                             aaGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGcccccaccccc                                                                          >  2:6930511/1‑150 (MQ=255)
                                                                               ggAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAg                                                                        >  2:6432333/1‑150 (MQ=255)
                                                                               ggAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAg                                                                        >  2:1262665/1‑150 (MQ=255)
                                                                                    ttGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCt                                                                   <  1:2480205/150‑1 (MQ=255)
                                                                                      gCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTc                                                                 <  1:6726044/150‑1 (MQ=255)
                                                                                      gCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTc                                                                 <  2:196141/150‑1 (MQ=255)
                                                                                        cGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCt                                                               <  2:7408294/150‑1 (MQ=255)
                                                                                        cGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCt                                                               >  2:2954698/1‑150 (MQ=255)
                                                                                         gagaAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTcc                                                                                                        >  2:3026813/1‑108 (MQ=255)
                                                                                         gagaAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTg                                                              >  2:5290751/1‑150 (MQ=255)
                                                                                         gagaAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTg                                                              <  1:3347392/150‑1 (MQ=255)
                                                                                         gagaAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTg                                                              >  1:5572302/1‑150 (MQ=255)
                                                                                           gaAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCt                                                                                                    >  2:4100335/1‑110 (MQ=255)
                                                                                            aaaaaGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCAc                                                           >  2:3159725/1‑150 (MQ=255)
                                                                                            aaaaaGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCAc                                                           >  1:225113/1‑150 (MQ=255)
                                                                                               aaGGAACAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCATTTc                                                        <  1:5742615/150‑1 (MQ=255)
                                                                                                aGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCg                                                       <  1:5363651/150‑1 (MQ=255)
                                                                                                 ggAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATGATGCGcccc                                                                                                                                              >  2:4697544/1‑62 (MQ=255)
                                                                                                 ggAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCgg                                                      >  1:3198726/1‑150 (MQ=255)
                                                                                                  gAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGc                                                     >  1:1282156/1‑150 (MQ=255)
                                                                                                   aaaaaaCGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCt                                                    >  2:2656778/1‑150 (MQ=255)
                                                                                                     aaaaCGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCg                                                  >  1:7353514/1‑150 (MQ=255)
                                                                                                        aCGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTg                                               <  2:1039294/150‑1 (MQ=255)
                                                                                                         cGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTgg                                              <  1:1440087/150‑1 (MQ=255)
                                                                                                           gTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTc                                            <  1:401255/150‑1 (MQ=255)
                                                                                                                  aaGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGccc                                     <  1:2985697/150‑1 (MQ=255)
                                                                                                                    ggCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTc                                   <  1:6699303/150‑1 (MQ=255)
                                                                                                                     gCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATGTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCg                                  >  2:7002959/1‑150 (MQ=255)
                                                                                                                      cTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACATCGGCTCGCTGGTCGCCGCCCTCGc                                 >  1:5109483/1‑150 (MQ=255)
                                                                                                                       tGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGcc                                >  2:4429985/1‑150 (MQ=255)
                                                                                                                       tGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGcc                                >  1:3880832/1‑150 (MQ=255)
                                                                                                                       tGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGcc                                >  2:5679491/1‑150 (MQ=255)
                                                                                                                          tctttcttCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGcctcc                             <  1:4265117/150‑1 (MQ=255)
                                                                                                                            tttcttCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTg                           >  2:3215095/1‑150 (MQ=255)
                                                                                                                               cttcCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCCGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCt                        <  2:4313352/150‑1 (MQ=255)
                                                                                                                                ttcCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTc                       <  1:7124761/150‑1 (MQ=255)
                                                                                                                                ttcCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTc                       >  2:2265592/1‑150 (MQ=255)
                                                                                                                                ttcCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGAACGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTc                       >  1:7209388/1‑150 (MQ=255)
                                                                                                                                 tcCACCTGAACGGGGCGCATCCTGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCGACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCg                      <  2:5610042/150‑1 (MQ=255)
                                                                                                                                   cACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGAt                    <  2:4533512/150‑1 (MQ=255)
                                                                                                                                   cACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGAt                    <  1:4502142/150‑1 (MQ=255)
                                                                                                                                      tgaaaCGGGGCGCATGATGCGCCCCGTTCCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGcc                 >  2:4868220/4‑150 (MQ=255)
                                                                                                                                      tgaaaCGGGGCGCATGATGCGCCCCGTTCCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGcc                 >  2:4514615/4‑150 (MQ=255)
                                                                                                                                      cTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGcc                 >  2:7243788/1‑150 (MQ=255)
                                                                                                                                       tGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGccc                <  2:3880832/150‑1 (MQ=255)
                                                                                                                                        gAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGCCcg               >  1:6574201/1‑150 (MQ=255)
                                                                                                                                        aaaCGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGCCcg               >  1:1235760/2‑150 (MQ=255)
                                                                                                                                          aCGGGGCGCATGATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGCCcgcg             >  1:6514727/1‑150 (MQ=255)
                                                                                                                                          aCGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGCCcgcg             >  2:5089670/1‑150 (MQ=255)
                                                                                                                                            ggggCGCATGATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGCCCGCGcc           >  2:4757073/1‑150 (MQ=255)
                                                                                                                                             gggCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGCCCGCGCCg          <  1:5205262/150‑1 (MQ=255)
                                                                                                                                              ggCGCATCATGCGCCCCGTTACATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGCCCGCGCCGt         <  2:6656226/150‑1 (MQ=255)
                                                                                                                                                   atcatGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGTGCGCTACGACCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGCCCGCGCCGTCAACg    >  2:3692267/1‑150 (MQ=255)
                                                                                                                                                    tcatGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGCCCGCGCCGTCAACgg   >  2:2030729/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGTGCGAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGCCCGCGCCGTCAACGGC  >  NC_002947/5337093‑5337391

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: