Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 406,031 C→G 6.3% S9T (AGC→ACC)  PP_0337 ← sensory box protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947406,0310CG 6.3% 202.1 / 3.3 63S9T (AGC→ACC) PP_0337sensory box protein
Reads supporting (aligned to +/- strand):  ref base C (27/32);  new base G (2/2);  total (29/34)
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.20e-01

CTCGCCGGTTTTCAGGCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCTGTGGACAGGTTATCGGCTGA  >  NC_002947/405886‑406176
                                                                                                                                                 |                                                                                                                                                 
cTCGCCGGTTTTCAGGCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTgg                                                                                                                                                <  2:32103/149‑1 (MQ=255)
cTCGCCGGTTTTCAGGCAGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTgg                                                                                                                                                <  2:32107/149‑1 (MQ=255)
        ttttCAGGCCGAACTTGCTCTCACAGCCTTGGTCAAGGTAGATCATGGGACAGGTGGCTTCTTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCAGCCGTGTggcggccgag                                                                                                                                        <  3:53592/149‑5 (MQ=255)
            cAGGCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGtt                                                                                                                                    >  5:297037/1‑149 (MQ=255)
             aGGCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGtt                                                                                                                                    <  7:167130/148‑1 (MQ=255)
               gCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGc                                                                                                                                 <  2:52889/149‑1 (MQ=255)
                ccGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCt                                                                                                                                >  5:26285/1‑149 (MQ=255)
                  gAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCttt                                                                                                                              <  5:292352/149‑1 (MQ=255)
                        gcgcTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCAt                                                                                                                          <  6:26285/147‑1 (MQ=255)
                        gcgcTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCg                                                                                                                        <  7:66063/149‑1 (MQ=255)
                          gcTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGaa                                                                                                                      <  5:29776/149‑1 (MQ=255)
                             cacaGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGAc                                                                                                                   >  5:51709/1‑149 (MQ=255)
                              acaGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACt                                                                                                                  >  8:4156/1‑149 (MQ=255)
                              acaGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACt                                                                                                                  >  6:120285/1‑149 (MQ=255)
                               caGGCTGTGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTc                                                                                                                 <  6:51709/149‑1 (MQ=255)
                                  gCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCt                                                                                                              >  2:136387/1‑149 (MQ=255)
                                   cTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGGCACACGGCTGGCGGCATCGGGTTTGCTTTTCATCGAAGACTCGCtt                                                                                                             <  5:20164/149‑1 (MQ=255)
                                       gTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATAGGGTTGGCTTTTCATCGATGACTCGCTTg                                                                                                            >  8:116347/1‑146 (MQ=255)
                                           aGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAgg                                                                                                      >  3:259093/1‑148 (MQ=255)
                                            ggTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAggg                                                                                                     >  8:236491/1‑148 (MQ=255)
                                             gTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATAGGGTTGGCTGTTCATCGAAGACTCGCTTGAAGAGGGaa                                                                                                   >  8:273758/1‑149 (MQ=255)
                                                       gCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTc                                                                                                                 <  3:176526/125‑1 (MQ=255)
                                                       gCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTc                                                                                                                 >  4:176526/1‑125 (MQ=255)
                                                       gCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTtg                                                                                         >  3:231008/1‑149 (MQ=255)
                                                        ccaaGGTTGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTtgt                                                                                        <  7:112791/146‑1 (MQ=255)
                                                                gCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTggc                                                                                                                                               <  6:280108/86‑1 (MQ=255)
                                                                gCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTggc                                                                                                                                               >  5:280108/1‑86 (MQ=255)
                                                                gCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGGTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTgcgg                                                                                <  7:273758/149‑1 (MQ=255)
                                                                     gTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGa                                                                            <  6:298204/148‑1 (MQ=255)
                                                                      ttGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACTCGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGAtt                                                                          <  1:136387/149‑1 (MQ=255)
                                                                                cGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAag                                                                                                         <  1:172027/108‑1 (MQ=255)
                                                                                cGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAag                                                                                                         >  2:172027/1‑108 (MQ=255)
                                                                                 gAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGtt                                                                <  5:120285/148‑1 (MQ=255)
                                                                                  aaaCGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCGTAAAGAGGGAAACCGGCCGTGAGGCTGCGGCCGATTTGGAAAGAtta                                                               >  7:191780/1‑147 (MQ=255)
                                                                                          caTGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTc                                                      <  7:4156/149‑1 (MQ=255)
                                                                                               cGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAg                                                 >  2:231560/1‑149 (MQ=255)
                                                                                                gAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAgg                                                >  6:64721/1‑149 (MQ=255)
                                                                                                 agacgcgAGGGCACGGGTAACTTTGTCACGACCTCGGCCGACACACGGCTGGCGGCATCGGGTTTGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAggg                                               <  8:120803/146‑1 (MQ=255)
                                                                                                    cgcgAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCat                                                            <  8:139619/133‑1 (MQ=255)
                                                                                                    cgcgAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCat                                                            >  7:139619/1‑133 (MQ=255)
                                                                                                    cgcgAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTg                                            <  8:66488/149‑1 (MQ=255)
                                                                                                      cgAGGGCACGTGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCCTCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTgcgg                                                                                >  2:41811/1‑111 (MQ=255)
                                                                                                            cACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCgg                                    <  5:64721/149‑1 (MQ=255)
                                                                                                                 gTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGACCGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATg                               >  7:285463/1‑149 (MQ=255)
                                                                                                                          tCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCt                      >  2:34864/1‑149 (MQ=255)
                                                                                                                          tCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCt                      >  1:288970/1‑149 (MQ=255)
                                                                                                                            aCGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCtgt                    <  4:259093/149‑1 (MQ=255)
                                                                                                                            aCGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCtgt                    <  1:231560/149‑1 (MQ=255)
                                                                                                                            aCGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCtgt                    <  6:297037/149‑1 (MQ=255)
                                                                                                                             cGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAg                                                 <  8:217211/119‑1 (MQ=255)
                                                                                                                             cGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAg                                                 >  7:217211/1‑119 (MQ=255)
                                                                                                                              gACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAg                                                                   >  4:12039/1‑100 (MQ=255)
                                                                                                                              gACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAg                                                                   <  3:12039/100‑1 (MQ=255)
                                                                                                                               cccgatgccgccAGCCGTGTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCTGTGGa                 >  7:232958/7‑149 (MQ=255)
                                                                                                                                ccgatgccgccAGCCGTGTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGc                                      <  4:210193/122‑1 (MQ=255)
                                                                                                                                   atgccgccAGCCGTGTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCTGTGGACAgg             >  1:145810/3‑149 (MQ=255)
                                                                                                                                    ggCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCTGTGGACAg              <  4:231008/147‑1 (MQ=255)
                                                                                                                                    ggCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCTGTGGACAg              <  4:83057/147‑1 (MQ=255)
                                                                                                                                        gccACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCTGTGGACAGGTTATc        >  4:38859/1‑149 (MQ=255)
                                                                                                                                          cacaCGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAg                                                 >  2:117940/1‑106 (MQ=255)
                                                                                                                                          cacaCGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAg                                                 <  1:117940/106‑1 (MQ=255)
                                                                                                                                            cacGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCTGTGGACAGGTTATCGGCt    >  5:296344/1‑149 (MQ=255)
                                                                                                                                              cGGCTGGCGGCATCGGGTTGGCTTTTCATAGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCTGTGGACAGGTTATCGGCTGa  <  7:236491/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CTCGCCGGTTTTCAGGCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCTGTGGACAGGTTATCGGCTGA  >  NC_002947/405886‑406176

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