Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 406,026 C→G 6.7% V11L (GTG→CTG)  PP_0337 ← sensory box protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947406,0260CG 6.7% 180.3 / 3.4 60V11L (GTG→CTG) PP_0337sensory box protein
Reads supporting (aligned to +/- strand):  ref base C (26/30);  new base G (2/2);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.26e-01

CACACAGCTCGCCGGTTTTCAGGCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCTGTGGACAGGTTATCGGCT  >  NC_002947/405879‑406174
                                                                                                                                                   |                                                                                                                                                    
cacacaGCTCGCCGGTTTTCAGGCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCcaca                                                                                                                                                     >  3:83057/1‑149 (MQ=255)
   acaGCTCGCCGGTTTTCAGGCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACgg                                                                                                                                                  >  6:29776/1‑149 (MQ=255)
       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                  <  6:297037/149‑1 (MQ=255)
                                                                                                                                   aCGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCtgt                  <  1:231560/149‑1 (MQ=255)
                                                                                                                                   aCGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCtgt                  <  4:259093/149‑1 (MQ=255)
                                                                                                                                    cGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAg                                               >  7:217211/1‑119 (MQ=255)
                                                                                                                                    cGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAg                                               <  8:217211/119‑1 (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)
                                                                                                                                                   |                                                                                                                                                    
CACACAGCTCGCCGGTTTTCAGGCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACGACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAAGACTCGCTTGAAGAGGGAACGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCTGTGGACAGGTTATCGGCT  >  NC_002947/405879‑406174

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