Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,843,613 T→C N408S (AAC→AGC)  gacS ← sensor protein GacS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,843,6130TC98.8% 1156.5 / ‑0.3 331N408S (AAC→AGC) gacSsensor protein GacS
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base C (157/170);  total (159/172)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 331 total aligned reads displayed.

AGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGATCAGGTCCCGCAGGTTG  >  NC_002947/1843466‑1843759
                                                                                                                                                   |                                                                                                                                                  
 gtgtGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTAAACTTGATGGCGTTGCTGACCAGGCTgg                                                                                                                                                 <  1:100022/150‑1 (MQ=255)
  tgtgCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGt                                                                                                                                                >  1:4953965/1‑150 (MQ=255)
  tgtgCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGt                                                                                                                                                <  2:5909285/150‑1 (MQ=255)
  tgtgCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGt                                                                                                                                                >  1:1700351/1‑150 (MQ=255)
    tgCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGa                                                                                                                                              >  1:5852043/1‑150 (MQ=255)
       gcACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCATCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGa                                                                                                                                           >  2:6479121/1‑150 (MQ=255)
         acaTCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATc                                                                                                                                         >  1:1800310/1‑150 (MQ=255)
            tCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGc                                                                                                                                      <  1:5430615/150‑1 (MQ=255)
               tGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTg                                                                                                                                   <  1:5658025/150‑1 (MQ=255)
               tGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTg                                                                                                                                   <  2:2977365/150‑1 (MQ=255)
                gAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGa                                                                                                                                  <  2:2692044/150‑1 (MQ=255)
                gAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGa                                                                                                                                  <  2:564141/150‑1 (MQ=255)
                 aGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAg                                                                                                                                 <  2:6638796/150‑1 (MQ=255)
                 aGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAg                                                                                                                                 <  2:3721405/150‑1 (MQ=255)
                  ggTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGc                                                                                                                                <  1:4177825/150‑1 (MQ=255)
                        aGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAAc                                                                                                                          <  2:1075652/150‑1 (MQ=255)
                        aGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAAc                                                                                                                          >  1:3709359/1‑150 (MQ=255)
                        aGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAAc                                                                                                                          >  1:2070868/1‑150 (MQ=255)
                        aGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAAc                                                                                                                          <  2:5327186/150‑1 (MQ=255)
                          gCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACgg                                                                                                                        <  2:3506380/150‑1 (MQ=255)
                            cGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGt                                                                                                                      >  2:1896322/1‑150 (MQ=255)
                            cGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGt                                                                                                                      >  2:3338602/1‑150 (MQ=255)
                              aTGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCg                                                                                                                    >  2:1986214/1‑150 (MQ=255)
                              aTGCCGGTGTCCTGCACGCTGATCCGCAACTTCGCGCAGTCTTCGTCCTCGTCTTCCAGCATGGCGAGCACGACGATGGTGCCTACGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCg                                                                                                                    >  2:4675353/1‑150 (MQ=255)
                               tGCCGGGGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGc                                                                                                                   <  1:3450183/150‑1 (MQ=255)
                                 ccGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCt                                                                                                                 <  2:4880722/150‑1 (MQ=255)
                                  cGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTa                                                                                                                >  2:276019/1‑150 (MQ=255)
                                    gtgtCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTAcc                                                                                                              <  2:304940/150‑1 (MQ=255)
                                    gaagacTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTAcc                                                                                                              >  1:2706268/6‑150 (MQ=255)
                                     tgtACTGCACGCTGATGCGCAACTGCGCGCTGTCTTCGTGCTAGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGATCTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCTACGGGTCGCCTagaa                                                                                                             <  2:3112259/150‑4 (MQ=255)
                                         cTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGa                                                                                                         >  1:622341/1‑150 (MQ=255)
                                         cTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGa                                                                                                         >  1:6582948/1‑150 (MQ=255)
                                         cTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGa                                                                                                         >  1:5234187/1‑150 (MQ=255)
                                          tGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGaa                                                                                                        <  1:1488574/150‑1 (MQ=255)
                                          tGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGaa                                                                                                        <  1:2774124/150‑1 (MQ=255)
                                          tGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGaa                                                                                                        >  2:243633/1‑150 (MQ=255)
                                          tGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGaa                                                                                                        <  1:138139/150‑1 (MQ=255)
                                           gCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGacc                                                                                                       <  2:1622699/150‑3 (MQ=255)
                                            cACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGc                                                                                                      >  2:1710719/1‑150 (MQ=255)
                                                cTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGgt                                                                                                  <  2:4174546/150‑1 (MQ=255)
                                                  gATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGgtgt                                                                                                >  2:2384874/1‑150 (MQ=255)
                                                     gcgcAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTc                                                                                             >  2:2199688/1‑150 (MQ=255)
                                                     gcgcAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTc                                                                                             >  2:5730705/1‑150 (MQ=255)
                                                       gcAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTcgc                                                                                           <  2:4828572/150‑1 (MQ=255)
                                                        caccTGCGCTCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTcgcg                                                                                          <  2:5673454/147‑1 (MQ=255)
                                                         aaCTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGGAGCAGGTGTGTCGCgg                                                                                         >  1:5893110/1‑150 (MQ=255)
                                                         aaCTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCgg                                                                                         >  1:6774011/1‑150 (MQ=255)
                                                         aaCTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCgg                                                                                         >  1:2185186/1‑150 (MQ=255)
                                                            tGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAg                                                                                      <  1:5219175/150‑1 (MQ=255)
                                                               gcgcAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATc                                                                                   >  1:224118/1‑150 (MQ=255)
                                                                cgcAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCa                                                                                  <  2:2315357/150‑1 (MQ=255)
                                                                cgcAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCa                                                                                  >  1:1251032/1‑150 (MQ=255)
                                                                cgcAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCa                                                                                  <  1:3193667/150‑1 (MQ=255)
                                                                 gcAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAg                                                                                 >  1:5722249/1‑150 (MQ=255)
                                                                  cAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAgg                                                                                <  1:3305665/150‑1 (MQ=255)
                                                                      cTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCa                                                                            <  1:5588321/150‑1 (MQ=255)
                                                                      cTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCa                                                                            <  1:4483125/150‑1 (MQ=255)
                                                                          gTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCaa                                                                        >  1:2569400/1‑150 (MQ=255)
                                                                           tGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAAt                                                                       >  1:349531/1‑150 (MQ=255)
                                                                              tCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTcc                                                                    >  2:5257099/1‑150 (MQ=255)
                                                                               cGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCa                                                                   >  2:6395022/1‑150 (MQ=255)
                                                                               cGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCa                                                                   <  1:5022844/150‑1 (MQ=255)
                                                                                 tCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAgc                                                                 <  1:4681681/150‑1 (MQ=255)
                                                                                   ttCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGCGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAgctg                                                               >  1:3065155/1‑150 (MQ=255)
                                                                                      gaGCATGGCGCGCACGAGGTTGGTGCCTTCGCGGGTGATCTTGATGGCGTTGCTGTCCAGGCTGGTGAGGATCTGCTTGAGCCGCTACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCtt                                                            <  2:4464114/149‑1 (MQ=255)
                                                                                        gCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTgg                                                          <  1:335481/150‑1 (MQ=255)
                                                                                           tGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGt                                                       >  1:4513425/1‑150 (MQ=255)
                                                                                             gcgcgcACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTgg                                                     <  2:3803190/150‑1 (MQ=255)
                                                                                                 gcACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGggcgg                                                 >  1:4286747/1‑150 (MQ=255)
                                                                                                 gcACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGggcgg                                                 >  2:6192335/1‑150 (MQ=255)
                                                                                                    cgacgaTGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCg                                              <  1:2356684/150‑1 (MQ=255)
                                                                                                     gacgaTGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCgg                                             >  1:6250573/1‑150 (MQ=255)
                                                                                                     gacgaTGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCgg                                             <  1:6901977/150‑1 (MQ=255)
                                                                                                      acgaTGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGa                                            <  1:4408145/150‑1 (MQ=255)
                                                                                                         aTGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCg                                         <  2:701615/150‑1 (MQ=255)
                                                                                                          tGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGa                                        <  2:1513871/150‑1 (MQ=255)
                                                                                                          tGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGa                                        >  1:4018850/1‑150 (MQ=255)
                                                                                                          tGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGa                                        <  2:5183781/150‑1 (MQ=255)
                                                                                                          tGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGa                                        <  2:4018850/150‑1 (MQ=255)
                                                                                                           ggTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAg                                       >  1:6349454/1‑150 (MQ=255)
                                                                                                              gCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGAt                                    <  1:5510456/150‑1 (MQ=255)
                                                                                                               ccTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATg                                   <  1:1905655/150‑1 (MQ=255)
                                                                                                                 ttCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGt                                 <  1:533608/150‑1 (MQ=255)
                                                                                                                     cgGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAgg                             <  2:1258830/150‑1 (MQ=255)
                                                                                                                       ggTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGgt                           <  1:2962317/150‑1 (MQ=255)
                                                                                                                       ggTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGgt                           <  2:5128241/150‑1 (MQ=255)
                                                                                                                         tGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGgtgt                         <  1:3142310/150‑1 (MQ=255)
                                                                                                                            aCTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTCCt                      <  1:363400/150‑1 (MQ=255)
                                                                                                                            aCTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTCCt                      <  2:2789289/150‑1 (MQ=255)
                                                                                                                                gATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGAt                  <  2:4201928/150‑1 (MQ=255)
                                                                                                                                 aTGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGATc                 <  2:5982324/150‑1 (MQ=255)
                                                                                                                                  tGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGATCa                >  2:1008940/1‑150 (MQ=255)
                                                                                                                                  tGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGATCa                <  2:5440718/150‑1 (MQ=255)
                                                                                                                                      gTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGATCAGGTc            <  1:24458/150‑1 (MQ=255)
                                                                                                                                           tGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGATCAGGTCCCGCa       <  1:812332/150‑1 (MQ=255)
                                                                                                                                             aCCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGATCAGGTCCCGCAgg     >  1:4231244/1‑150 (MQ=255)
                                                                                                                                             aCCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGATCAGGTCCCGCAgg     >  1:3705195/1‑150 (MQ=255)
                                                                                                                                               cAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGATCAGGTCCCGCAGGtt   <  2:6582948/150‑1 (MQ=255)
                                                                                                                                               cAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGATCAGGTCCCGCAGGtt   <  2:6745770/150‑1 (MQ=255)
                                                                                                                                                aGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGATCAGGTCCCGCAGGTTg  >  2:4758258/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGATCAGGTCCCGCAGGTTG  >  NC_002947/1843466‑1843759

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: ‑