Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,843,370 G→T 54.9% T475N (ACC→AAC)  gacS ← sensor protein GacS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,843,3700GT54.9% 12.6 / 84.7 62T475N (ACC→AAC) gacSsensor protein GacS
Reads supporting (aligned to +/- strand):  ref base G (12/16);  new base T (18/16);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

CTGCGGTGGCAGCGCTTCGGCATCGTCATGAGCCTTGGGCAGGTTGAGGCTGATCCAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGC  >  NC_002947/1843222‑1843518
                                                                                                                                                    |                                                                                                                                                     
cTGCGGTGGCAGCGCTTCGGCATCGTCATGAGCCTTGGGCAGGTTGAGGCTGATCCAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCgg                                                                                                                                                       >  5:185992/1‑149 (MQ=255)
        gCAGCGCTTCGGCATCGTCATGAGCCTTGGGCAGGTTGAGGCTGATCCAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTgccgcc                                                                                                                                                <  8:131633/148‑1 (MQ=255)
         cAGCGCTTCGGCATCGGCATGAGCCTTGGGCAGGTTGAGGCTGATCCACAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGTTGCCGCCgg                                                                                                                                              <  7:207295/149‑1 (MQ=255)
               ttCGGCATCGTCATGAGCCTTGGGCAGGTTGAGGCTGATCCAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGc                                                                                                                                        >  7:8992/1‑149 (MQ=255)
                 cGGCATCGTCATGAGCCTTGGGCAGGTTGAGGCTGATCCAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGTTGCCGCCGGGCTGGc                                                                                                                                        <  6:190405/147‑1 (MQ=255)
                            tGAGCCTTGGGCAGGTTGAGGCTGATCCAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGtt                                                                                                                           <  5:43139/149‑1 (MQ=255)
                               gCCTTGGGCAGGTTGAGGCTGATCCAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGTTGCCGCCGGGCTGGCGGGCCAGTGAGTTGt                                                                                                                         >  5:47919/1‑148 (MQ=255)
                                ccTTGGGCAGGTTGAGGCTGATCCAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGt                                                                                                                         >  7:270775/1‑147 (MQ=255)
                                ccTTGGGCAGGTTGAGGCTGATCCAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCACGCGCTTGGAAATCACCAGCCCCAGACCGTTGCCGCCGGGCTGGCGGGCCAGTGAGTTGt                                                                                                                         >  7:244508/1‑147 (MQ=255)
                                      gCAGGTTGAGGCTGATCCAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCCATCTCGCCGCCCATCTGTTCGATCAGGCGCTAGGAAATCACCAGCCCCAGACCGGTGCCGCCGCGCTGGCGGGCCAGTTAGTTGTCAGCCTgc                                                                                                                 >  8:162744/1‑148 (MQ=255)
                                           ttGAGGCTGATCCAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGTTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGaa                                                                                                            <  2:167315/149‑1 (MQ=255)
                                              aGGCTGATCCAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGc                                                                                                         >  6:52358/1‑149 (MQ=255)
                                                 cTGATCCAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTg                                                                                                      <  6:185992/149‑1 (MQ=255)
                                                    aTCCAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGTTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGaa                                                                                                   >  5:252704/1‑149 (MQ=255)
                                                    aTCCAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGaa                                                                                                   <  7:162744/149‑1 (MQ=255)
                                                      ccAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACa                                                                                                 <  1:153077/149‑1 (MQ=255)
                                                       cAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGTCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACa                                                                                                 <  1:50194/148‑1 (MQ=255)
                                                       cAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACa                                                                                                 <  5:52358/148‑1 (MQ=255)
                                                       cAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAg                                                                                                >  6:125096/1‑149 (MQ=255)
                                                        aGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGTTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAgt                                                                                               >  2:146848/1‑149 (MQ=255)
                                                        aGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGTTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAgt                                                                                               <  3:72105/149‑1 (MQ=255)
                                                          aaCTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAgtgt                                                                                             <  2:176197/149‑1 (MQ=255)
                                                               cGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTgcgc                                                                                         >  8:286291/1‑148 (MQ=255)
                                                                  aCCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGTTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACAt                                                                                     >  2:203377/1‑149 (MQ=255)
                                                                     ttCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCCGGCCACTGAGTTGTCAGCCTGGCTGAAGGCCTGGATCCGTGTGCGCACCTccg                                                                                  >  8:305203/1‑148 (MQ=255)
                                                                      tCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGTTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCt                                                                                 >  4:257693/1‑148 (MQ=255)
                                                                        gCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGTTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTg                                                                                >  3:249995/1‑147 (MQ=255)
                                                                        gCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGTTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAg                                                                              <  6:47919/149‑1 (MQ=255)
                                                                        gCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGa                                                                               <  8:8992/148‑1 (MQ=255)
                                                                        gCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGa                                                                               <  7:286291/148‑1 (MQ=255)
                                                                        gCCCGGGGTACTGTCGACACCGAGCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGTTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTg                                                                                <  4:249995/147‑1 (MQ=255)
                                                                         cccGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGTTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAgg                                                                             <  1:146848/149‑1 (MQ=255)
                                                                          ccGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGt                                                                            >  3:254316/1‑149 (MQ=255)
                                                                               gTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGTTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAg                                                                       <  3:53441/149‑1 (MQ=255)
                                                                                          acCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGTTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCgg                                                            >  4:198779/1‑149 (MQ=255)
                                                                                            cGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGTTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGgtg                                                          >  5:16078/1‑149 (MQ=255)
                                                                                             gATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGTTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTAAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGgtgt                                                         >  2:254674/1‑149 (MQ=255)
                                                                                                   gccgccCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGTTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTcc                                                       <  3:257693/145‑1 (MQ=255)
                                                                                                                  gATCAGGCGCTTGGAAATCACCAGCCCCAGACCGTTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACt                                    >  5:54652/1‑149 (MQ=255)
                                                                                                                    tCAGGCGCTTGGAAATCACCAGCCCCAGACCGTTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTgc                                  >  1:258807/1‑149 (MQ=255)
                                                                                                                      aGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTgcgc                                <  1:134440/149‑1 (MQ=255)
                                                                                                                        gcgcTTGGAAATCACCAGCCCCAGACCGTTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTg                                                     >  4:87434/1‑126 (MQ=255)
                                                                                                                        gcgcTTGGAAATCACCAGCCCCAGACCGTTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTg                                                     <  3:87434/126‑1 (MQ=255)
                                                                                                                        gcgcTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTgcgcg                               >  5:65878/1‑148 (MQ=255)
                                                                                                                          gcTTGGAAATCACCAGCCCCAGACCGTTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCa                             <  1:254674/149‑1 (MQ=255)
                                                                                                                          gcTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAg                            >  1:252693/1‑149 (MQ=255)
                                                                                                                               gAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTc                       <  6:65878/149‑1 (MQ=255)
                                                                                                                                  aTCACCAGCCCCAGACCGTTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTg                    >  1:120544/1‑149 (MQ=255)
                                                                                                                                  aTCACCAGCCCCAGACCGTTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTg                    >  1:188779/1‑149 (MQ=255)
                                                                                                                                  aTCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTTCGCGCAGTCTTCGTg                    <  8:270775/149‑1 (MQ=255)
                                                                                                                                     aCCAGCCCCAGACCGTTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTc                 >  3:124353/1‑149 (MQ=255)
                                                                                                                                     aCCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTc                 >  2:145095/1‑149 (MQ=255)
                                                                                                                                        aGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTc              >  5:199900/1‑149 (MQ=255)
                                                                                                                                           cccAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTc           <  4:254316/149‑1 (MQ=255)
                                                                                                                                           cccAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTc           <  4:136748/149‑1 (MQ=255)
                                                                                                                                            ccAGACCGTTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTc                 >  3:279176/1‑142 (MQ=255)
                                                                                                                                            ccAGACCGTTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTc                 <  4:279176/142‑1 (MQ=255)
                                                                                                                                            ccAGACCGTTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTcc          <  6:252704/149‑1 (MQ=255)
                                                                                                                                             cAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCa         <  5:125096/149‑1 (MQ=255)
                                                                                                                                               gACCGTTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGc       <  6:116247/149‑1 (MQ=255)
                                                                                                                                                aCCGTTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCGGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTGCGTGCTCGGCTTCCAGc       <  8:244508/148‑1 (MQ=255)
                                                                                                                                                 ccGTTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTc                 >  5:12068/1‑137 (MQ=255)
                                                                                                                                                 ccGTTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTc                 <  6:12068/137‑1 (MQ=255)
                                                                                                                                                    gTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGgc  <  5:8402/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTGCGGTGGCAGCGCTTCGGCATCGTCATGAGCCTTGGGCAGGTTGAGGCTGATCCAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGC  >  NC_002947/1843222‑1843518

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