Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,160,662 G→C 5.7% intergenic (+28/+93) PP_3662 → / ← PP_3663 hypothetical protein/GGDEF domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,160,6620GC 5.7% 229.7 / 4.0 70intergenic (+28/+93)PP_3662/PP_3663hypothetical protein/GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (32/34);  new base C (2/2);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.69e-01

GGTGACAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCCGCTCGAAAAACCGAATGCTGTCGCGCC  >  NC_002947/4160518‑4160807
                                                                                                                                                 |                                                                                                                                                 
ggTGACAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATCCATGc                                                                                                                                               <  3:66616/149‑1 (MQ=255)
     cAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATCCATGCCCGTc                                                                                                                                          <  6:84084/149‑1 (MQ=255)
         gcCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAAc                                                                                                                                      <  5:194391/149‑1 (MQ=255)
           cGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAAcgc                                                                                                                                    >  8:161297/1‑149 (MQ=255)
               cTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCgggg                                                                                                                                <  3:419258/149‑1 (MQ=255)
               cTGGACCCGCGGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCggg                                                                                                                                 <  6:253361/148‑1 (MQ=255)
                tGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCggggg                                                                                                                               >  2:284714/1‑149 (MQ=255)
                  gACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGcc                                                                                                                             >  3:278401/1‑149 (MQ=255)
                   aCCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCg                                                                                                                            >  3:135031/1‑149 (MQ=255)
                     ccGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAAcgc                                                                                                                                    <  4:412375/139‑1 (MQ=255)
                     ccGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAAcgc                                                                                                                                    >  3:412375/1‑139 (MQ=255)
                     ccGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCggg                                                                                                                          >  8:1847/1‑149 (MQ=255)
                       gCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTCGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGt                                                                                                                         >  4:119014/1‑148 (MQ=255)
                           ggCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGTCGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACgg                                                                                                                    <  1:284714/149‑1 (MQ=255)
                           ggCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACg                                                                                                                     <  4:140217/148‑1 (MQ=255)
                               ggCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGc                                                                                                                <  4:50182/149‑1 (MQ=255)
                                 ctgctgCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGt                                                                                                              <  4:306306/149‑1 (MQ=255)
                                 ctgctgCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGt                                                                                                              <  7:161297/149‑1 (MQ=255)
                                    ctgcAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTc                                                                                                             <  4:293139/147‑1 (MQ=255)
                                      gcAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTGTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTa                                                                                                         >  2:368727/1‑149 (MQ=255)
                                       cAGGCATTCGTCGACCATCACCTGATGAATCTGCCCGGAGGCGCGGCGTATGTACCTTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGATGACGTGAATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGATCGTCTGTAc                                                                                                        <  8:234325/149‑1 (MQ=255)
                                       cAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGAGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTAc                                                                                                        <  8:240658/149‑1 (MQ=255)
                                        aGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACg                                                                                                       >  2:237707/1‑149 (MQ=255)
                                         ggCATTCGGCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACgg                                                                                                      <  5:209608/149‑1 (MQ=255)
                                          gCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGa                                                                                                     >  3:70302/1‑149 (MQ=255)
                                            aTTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGAt                                                                                                    >  5:184011/1‑148 (MQ=255)
                                              tcgtcgACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTAATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTa                                                                                                         >  1:303837/1‑142 (MQ=255)
                                              tcgtcgACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTa                                                                                                         <  2:303837/141‑1 (MQ=255)
                                               cgtcgACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCAc                                                                                                 <  1:237707/148‑1 (MQ=255)
                                                gtcgACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCgg                                                                                                                           <  1:254691/121‑1 (MQ=255)
                                                gtcgACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCgg                                                                                                                           >  2:254691/1‑121 (MQ=255)
                                                    aCCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACAt                                                                                           >  2:252329/1‑149 (MQ=255)
                                                     ccAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATg                                                                                          >  4:136800/1‑149 (MQ=255)
                                                      cAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATgg                                                                                         >  6:207633/1‑149 (MQ=255)
                                                      cAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATgg                                                                                         >  4:345702/1‑149 (MQ=255)
                                                       aGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATggg                                                                                        <  3:345702/149‑1 (MQ=255)
                                                         cACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTTTCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGt                                                                                       >  3:15286/1‑148 (MQ=255)
                                                         cACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGtt                                                                                      <  1:271246/149‑1 (MQ=255)
                                                                  aaGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGaaa                                                                              <  4:166703/148‑1 (MQ=255)
                                                                      tGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTcc                                                                         <  4:70302/149‑1 (MQ=255)
                                                                         ccGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTcccc                                                                       >  5:303599/1‑148 (MQ=255)
                                                                         ccGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTcccc                                                                       <  6:303599/148‑1 (MQ=255)
                                                                                   ggcgTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGa                                                            <  1:298423/149‑1 (MQ=255)
                                                                                      gTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTc                                                         >  4:373482/1‑149 (MQ=255)
                                                                                       tATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGGAAGCTCCCCCACCTAGGTGACTCg                                                        >  2:336842/1‑149 (MQ=255)
                                                                                            tACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTc                                                     >  4:141656/1‑148 (MQ=255)
                                                                                               cGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTg                                                 >  3:127712/1‑149 (MQ=255)
                                                                                                   cTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCg                                              <  4:278401/148‑1 (MQ=255)
                                                                                                   cTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCg                                              <  6:184011/148‑1 (MQ=255)
                                                                                                      tCATGTGGCAACCTGACCCGCTGCTTTACTTAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGt                                          <  7:35773/149‑1 (MQ=255)
                                                                                                          gTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGc                                                       >  2:28978/1‑132 (MQ=255)
                                                                                                          gTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGc                                                       <  1:28978/132‑1 (MQ=255)
                                                                                                             gCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTACGTGACTCGCTCACTGCCGGAGTTCAacca                                   >  1:396607/1‑149 (MQ=255)
                                                                                                               aaCCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTAAaccacc                                 >  5:352833/1‑149 (MQ=255)
                                                                                                                 ccTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCAccgc                               <  3:119014/149‑1 (MQ=255)
                                                                                                                       ccGCTGCTTTACTGAGGACGGGCATGGATGGCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCCGCTc                         >  1:36211/1‑149 (MQ=255)
                                                                                                                       ccGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCCGCTc                         <  5:207633/149‑1 (MQ=255)
                                                                                                                                 aCTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCCGCTCGAAAAACCg                <  3:136800/148‑1 (MQ=255)
                                                                                                                                  cTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCCGCTCGAAAAACCGaa              <  7:1847/149‑1 (MQ=255)
                                                                                                                                  cTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCCGCTCGAAAAACCGa               <  4:15286/148‑1 (MQ=255)
                                                                                                                                       gACGGGCATCCATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCCGCTCGAAAAACCGAATGCTg         >  2:145342/1‑149 (MQ=255)
                                                                                                                                       gACGGGCATCCATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCCGCTCGAAAAACCGAATGCTg         >  3:197360/1‑149 (MQ=255)
                                                                                                                                          gggCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCCGCTCGaaaaa                   <  2:105304/136‑1 (MQ=255)
                                                                                                                                          gggCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCCGCTCGaaaaa                   >  1:105304/1‑136 (MQ=255)
                                                                                                                                          gggCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCCGCTCGAAAAACCGAATGCTGTcg      >  2:365999/1‑149 (MQ=255)
                                                                                                                                           ggCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCCGCTCGAAAAACCGAATGCTGTcgc     <  4:135031/149‑1 (MQ=255)
                                                                                                                                           ggCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCCGCTCGAAAAACCGAATGCTGTcgc     <  3:373482/149‑1 (MQ=255)
                                                                                                                                            gCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCCGCTCGAAAAACCGAATGCTGTcgcg    >  6:198958/1‑149 (MQ=255)
                                                                                                                                            gCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCCGCTCGAAAAACCGAATGCTGTcgcg    >  3:318603/1‑149 (MQ=255)
                                                                                                                                              aTGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCCGCTCGAAAAACCGAATGCTGTCGCGcc  >  5:331007/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GGTGACAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCCGCTCGAAAAACCGAATGCTGTCGCGCC  >  NC_002947/4160518‑4160807

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: