Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,334,215 C→G 5.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,334,2150CG 5.1% 1058.5 / 31.9 311intergenic (+28/‑21)dksA/gluQRNA polymerase‑binding transcription factor DksA/glutamyl‑Q tRNA(Asp) synthetase
Reads supporting (aligned to +/- strand):  ref base C (159/136);  new base G (6/10);  total (165/146)
Fisher's exact test for biased strand distribution p-value = 2.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 333 total aligned reads displayed.

GGTGCGAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGCCCGCGCCGTCAACGGC  >  NC_002947/5334066‑5334364
                                                                                                                                                     |                                                                                                                                                     
   gCGAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATg                                                                                                                                                    >  2:3256697/1‑150 (MQ=255)
    cGAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATgc                                                                                                                                                   <  1:6476052/150‑1 (MQ=255)
    cGAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATgc                                                                                                                                                   >  2:1561702/1‑150 (MQ=255)
     gAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATgcg                                                                                                                                                  >  1:1047735/1‑150 (MQ=255)
      aGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATgcgc                                                                                                                                                 <  2:8271755/150‑1 (MQ=255)
        tCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATGATGCGccc                                                                                                                                               <  2:544233/150‑1 (MQ=255)
         cGTGCGGCATCGAGATCGGCCTGCGCCGCCTGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCg                                                                                                                                             <  2:7414061/150‑1 (MQ=255)
         cGTGCGGCATCGAGATCGGCCGGCGCCGCCTGGAAGCCCGCCCGCCCGCCGACCTGTGCTTCGACTGCAAGGACATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGcccc                                                                                                                                              <  1:2680770/150‑1 (MQ=255)
             cGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGttt                                                                                                                                          <  1:1984019/150‑1 (MQ=255)
                cATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCAt                                                                                                                                       <  1:8982036/150‑1 (MQ=255)
                cATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCAt                                                                                                                                       <  1:3761147/150‑1 (MQ=255)
                 aTCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCAtt                                                                                                                                      <  1:1653838/150‑1 (MQ=255)
                 aTCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCAtt                                                                                                                                      <  2:3159264/150‑1 (MQ=255)
                 aTCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCAtt                                                                                                                                      >  1:6119214/1‑150 (MQ=255)
                  tCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCAttt                                                                                                                                     <  2:1742653/150‑1 (MQ=255)
                     agaTCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTtata                                                                                                                                  >  2:2752624/1‑150 (MQ=255)
                     agaTCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTtata                                                                                                                                  >  2:813855/1‑150 (MQ=255)
                          ggCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGccc                                                                                                                             <  2:2534216/150‑1 (MQ=255)
                           gCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCa                                                                                                                            >  1:1840280/1‑150 (MQ=255)
                           gCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCa                                                                                                                            <  1:2554647/150‑1 (MQ=255)
                            ccTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAg                                                                                                                           >  1:1288094/1‑150 (MQ=255)
                                 gccgccTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGa                                                                                                                      <  2:7284669/150‑1 (MQ=255)
                                   cgccTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGAcc                                                                                                                    >  2:5958306/1‑150 (MQ=255)
                                    gccTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCa                                                                                                                   >  1:7490118/1‑150 (MQ=255)
                                       tGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGa                                                                                                                <  1:5727413/150‑1 (MQ=255)
                                        ggAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGAc                                                                                                               >  2:6482294/1‑150 (MQ=255)
                                         gAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGAcc                                                                                                              <  2:8527211/150‑1 (MQ=255)
                                            gcccgcccTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGAc                                                                                                           <  1:2796021/150‑1 (MQ=255)
                                             cccgcccTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACt                                                                                                          >  2:674696/1‑150 (MQ=255)
                                               cgcccTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTcc                                                                                                        >  1:7359949/1‑150 (MQ=255)
                                               cgcccTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTcc                                                                                                        >  1:71110/1‑150 (MQ=255)
                                                gcccTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCa                                                                                                       <  2:352330/150‑1 (MQ=255)
                                                     aCCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTAc                                                                                                  >  1:30301/1‑150 (MQ=255)
                                                        gccgACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATc                                                                                               >  2:2411248/1‑150 (MQ=255)
                                                              cTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGgcgc                                                                                         <  1:3342203/150‑1 (MQ=255)
                                                                 tgCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTc                                                                                      >  2:8784807/1‑150 (MQ=255)
                                                                 tgCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTc                                                                                      <  1:2386973/150‑1 (MQ=255)
                                                                      cGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCCGCTACATCGGGCGCTTCGcccc                                                                                 >  1:6282843/1‑150 (MQ=255)
                                                                      cGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGcccc                                                                                 >  1:6974496/1‑150 (MQ=255)
                                                                      cGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGcccc                                                                                 <  1:1470530/150‑1 (MQ=255)
                                                                      cGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGcccc                                                                                 <  1:8632864/150‑1 (MQ=255)
                                                                        aCTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGccccca                                                                               <  2:629878/150‑1 (MQ=255)
                                                                           gCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGcccccaccc                                                                            >  2:4037152/1‑150 (MQ=255)
                                                                            cAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTc                                                                                      <  1:7569294/139‑1 (MQ=255)
                                                                             aaGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGcccccaccccc                                                                          >  2:1887001/1‑150 (MQ=255)
                                                                             aaGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGcccccaccccc                                                                          >  1:658312/1‑150 (MQ=255)
                                                                             aaGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGcccccaccccc                                                                          >  1:6825322/1‑150 (MQ=255)
                                                                             aaGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGcccccaccccc                                                                          >  2:4266560/1‑150 (MQ=255)
                                                                              aGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGcccccaccccca                                                                         >  2:5791612/1‑150 (MQ=255)
                                                                                        cGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCt                                                               <  1:7354609/150‑1 (MQ=255)
                                                                                        cGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCt                                                               >  2:5322472/1‑150 (MQ=255)
                                                                                         gagaAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTg                                                              >  1:5059978/1‑150 (MQ=255)
                                                                                         gagaAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTg                                                              <  1:5907830/150‑1 (MQ=255)
                                                                                         gagaAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTg                                                              <  1:6181600/150‑1 (MQ=255)
                                                                                         gagaAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGATACATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTg                                                              <  2:7662467/150‑1 (MQ=255)
                                                                                          agaAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCAACAACAAAAACAAAACCCAACAACCGCCCCGTAACAATTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGc                                                             >  1:3987461/1‑150 (MQ=255)
                                                                                           gaAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCa                                                            <  1:4576763/150‑1 (MQ=255)
                                                                                            aaaaaGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCAc                                                           <  1:3611647/150‑1 (MQ=255)
                                                                                               aaGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTc                                                        >  2:7079478/1‑150 (MQ=255)
                                                                                                aGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCg                                                       >  2:1223836/1‑150 (MQ=255)
                                                                                                  gAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGc                                                     >  1:6051594/1‑150 (MQ=255)
                                                                                                      aaaCGGTCGGCAAAGGCTGATCTTTCGTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGc                                                 <  1:8724107/150‑1 (MQ=255)
                                                                                                          ggTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGt                                             >  2:2934476/1‑150 (MQ=255)
                                                                                                           gTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTc                                            <  1:5271389/150‑1 (MQ=255)
                                                                                                             cGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTcgc                                          <  2:7966813/150‑1 (MQ=255)
                                                                                                                  aaGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGccc                                     <  1:3887709/150‑1 (MQ=255)
                                                                                                                    ggCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTc                                   <  2:1362142/150‑1 (MQ=255)
                                                                                                                    ggCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTc                                   <  2:1047735/150‑1 (MQ=255)
                                                                                                                     gCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCg                                  >  1:8819215/1‑150 (MQ=255)
                                                                                                                      cTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGc                                 >  2:4161878/1‑150 (MQ=255)
                                                                                                                       tGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGcc                                <  1:4386308/150‑1 (MQ=255)
                                                                                                                           ctttcttCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGcctcct                            <  1:8067436/150‑1 (MQ=255)
                                                                                                                           ctttcttCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGcctcct                            >  2:2241780/1‑150 (MQ=255)
                                                                                                                           ctttcttCCACCTGAACGGGGCGCATCATGCGCACCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGcctcct                            >  1:4961210/1‑150 (MQ=255)
                                                                                                                            tttcttCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTc                                                                                      <  2:4689369/91‑1 (MQ=255)
                                                                                                                            tttcttCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTc                                                                                      >  1:4689369/1‑91 (MQ=255)
                                                                                                                            tttcttCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTg                           >  1:1424790/1‑150 (MQ=255)
                                                                                                                            tttcttCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTg                           >  1:547524/1‑150 (MQ=255)
                                                                                                                               cttcCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCt                        <  2:3757726/150‑1 (MQ=255)
                                                                                                                               cttcCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCt                        >  1:3544672/1‑150 (MQ=255)
                                                                                                                               cttcCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATATATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCt                        >  2:907966/1‑150 (MQ=255)
                                                                                                                                ttcCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGTGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTc                       >  2:4968779/1‑150 (MQ=255)
                                                                                                                                 tcCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCg                      <  2:8663024/150‑1 (MQ=255)
                                                                                                                                   cACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGAt                    <  2:8928069/150‑1 (MQ=255)
                                                                                                                                   cACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGAt                    <  2:7846620/150‑1 (MQ=255)
                                                                                                                                    aCCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTc                                                                 <  2:4093147/104‑1 (MQ=255)
                                                                                                                                     ccTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGc                  <  1:2300457/150‑1 (MQ=255)
                                                                                                                                      cTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGcc                 >  2:8447359/1‑150 (MQ=255)
                                                                                                                                      cTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGcc                 >  1:8226499/1‑150 (MQ=255)
                                                                                                                                       tGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGccc                >  2:6014424/1‑150 (MQ=255)
                                                                                                                                       tGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGccc                <  1:5240226/150‑1 (MQ=255)
                                                                                                                                       gaaaCGGGGCGCATGATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGccc                >  2:8083649/3‑150 (MQ=255)
                                                                                                                                       gaaaCGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGccc                >  1:1709060/3‑150 (MQ=255)
                                                                                                                                       agAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGccc                <  1:5565417/149‑1 (MQ=255)
                                                                                                                                          aCGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGCCcgcg             <  2:5487470/150‑1 (MQ=255)
                                                                                                                                           cGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGCCCgcgc            <  2:1701383/150‑1 (MQ=255)
                                                                                                                                              ggCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGCCCGCGCCGt         >  1:2733054/1‑150 (MQ=255)
                                                                                                                                                 gcATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGCCCGCGCCGTCaa      >  2:2111086/1‑150 (MQ=255)
                                                                                                                                                   atcatGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCggg                                                                                            <  1:4191846/62‑1 (MQ=255)
                                                                                                                                                    tcatGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGCCCGCGCCGTCAACgg   <  1:8376895/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGTGCGAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGCCCGCGCCGTCAACGGC  >  NC_002947/5334066‑5334364

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