Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,997,649 C→A 8.6% R167L (CGG→CTG)  PP_4405 ← sensory box protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,997,6490CA 8.6% 247.7 / 10.1 81R167L (CGG→CTG) PP_4405sensory box protein
Reads supporting (aligned to +/- strand):  ref base C (33/41);  new base A (5/2);  total (38/43)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.28e-01

TCCATGTCCATGGCAATCAGCGACAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATGCGCACCCAGCATGCGTGCCAC  >  NC_002947/4997502‑4997791
                                                                                                                                                   |                                                                                                                                              
tccatgtccatgGCAATCAGCGACAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGc                                                                                                                                              >  2:1806393/1‑150 (MQ=255)
tccatgtccatgGCAATCAGCGACAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGc                                                                                                                                              >  2:1850336/1‑150 (MQ=255)
   atgtccatgGCAATCAGCGACAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGc                                                                                                                                           >  2:1723677/1‑150 (MQ=255)
   atgtccatgGCAATCAGCGACAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGc                                                                                                                                           >  1:813678/1‑150 (MQ=255)
          tgGCAATCAGCGACAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAActgct                                                                                                                                    <  2:1984259/150‑1 (MQ=255)
          tgGCAATCAGCGACAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCAGCGGCAActgga                                                                                                                                    <  2:618209/150‑3 (MQ=255)
           ggCAATCAGCGACAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAActgctg                                                                                                                                   >  1:173938/1‑150 (MQ=255)
           ggCAATCAGCGAAAACGGATGACGGACGCGCTGAGCGGGTTGTCCTCCCAAGCGAAATCCGTTCTGCATCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAActgctg                                                                                                                                   <  2:2329634/150‑1 (MQ=255)
            gCAATCAGCGACAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTgg                                                                                                                                  <  2:813678/150‑1 (MQ=255)
            gCAATCAGCGACAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTgg                                                                                                                                  <  1:2058160/150‑1 (MQ=255)
            gCAATCAGCGACAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTgg                                                                                                                                  <  1:1896418/150‑1 (MQ=255)
                   gCGACAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCAccc                                                                                                                           <  1:2045574/150‑1 (MQ=255)
                     gACAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAg                                                                                                                         >  2:1018734/1‑150 (MQ=255)
                     gACAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAg                                                                                                                         <  2:2198427/150‑1 (MQ=255)
                      aCAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAgg                                                                                                                        <  1:1291787/150‑1 (MQ=255)
                       cAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGa                                                                                                                       <  2:85417/150‑1 (MQ=255)
                         aCGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAg                                                                                                                         <  2:2115918/146‑1 (MQ=255)
                         aCGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAg                                                                                                                         >  1:2115918/1‑146 (MQ=255)
                          cGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACaa                                                                                                                    <  1:2083453/150‑1 (MQ=255)
                            gCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAActgctg                                                                                                                                   >  1:1243892/1‑133 (MQ=255)
                            gCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAActgctg                                                                                                                                   <  2:1243892/133‑1 (MQ=255)
                               ggCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCAc                                                                                                               <  2:1982832/150‑1 (MQ=255)
                                gCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACg                                                                                                              >  1:1765145/1‑150 (MQ=255)
                                gCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACg                                                                                                              <  2:2404722/150‑1 (MQ=255)
                                        gcTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGtt                                                                                                      >  2:2241043/1‑150 (MQ=255)
                                         cTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTCACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTc                                                                                                     <  2:1739053/150‑1 (MQ=255)
                                            ggCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGc                                                                                                  >  1:574575/1‑150 (MQ=255)
                                               gcgcTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCAc                                                                                               <  1:1648697/150‑1 (MQ=255)
                                                     gCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGaaa                                                                                         <  2:1765145/150‑1 (MQ=255)
                                                       aCTGCAGCTGCAACCGGCTCTCCAGCCGGTAACGGTTGGCGACACTGTGCAACGCCACCTTTGTCGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACt                                                                                       >  2:1289622/1‑150 (MQ=255)
                                                       aCTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACt                                                                                       >  1:2085482/1‑150 (MQ=255)
                                                        cTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTTGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTc                                                                                      >  2:1880200/1‑150 (MQ=255)
                                                        cTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTCCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTc                                                                                      >  1:1475555/1‑150 (MQ=255)
                                                          ccAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGTTCAACGCATCGCGTTCGGCCAATTCAAGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCa                                                                                    <  1:771714/150‑1 (MQ=255)
                                                          ccAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCa                                                                                    <  1:1723677/150‑1 (MQ=255)
                                                           cAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCATCTGaaa                                                                                         <  1:842939/144‑1 (MQ=255)
                                                           cAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGaaa                                                                                         >  2:842939/1‑144 (MQ=255)
                                                               tGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGAc                                                                               >  1:2392027/1‑150 (MQ=255)
                                                                 cAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTg                                                                             <  1:2241043/150‑1 (MQ=255)
                                                                   aCGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGt                                                                           >  1:1489981/1‑150 (MQ=255)
                                                                   aCGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGt                                                                           >  1:1475503/1‑150 (MQ=255)
                                                                       gcTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGttctt                                                                       <  2:2085482/150‑1 (MQ=255)
                                                                       gcTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGttctt                                                                       >  2:497874/1‑150 (MQ=255)
                                                                       gcTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGttctt                                                                       <  1:2088743/150‑1 (MQ=255)
                                                                       gcTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGttctt                                                                       <  2:173938/150‑1 (MQ=255)
                                                                           tcCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGt                                                                   <  2:279251/150‑1 (MQ=255)
                                                                           tcCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGt                                                                   >  1:279251/1‑150 (MQ=255)
                                                                             cAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTc                                                                 <  2:108417/150‑1 (MQ=255)
                                                                             cAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTc                                                                 >  1:108417/1‑150 (MQ=255)
                                                                               gCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCa                                                               <  1:497874/150‑1 (MQ=255)
                                                                                    tAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTga                                                          <  2:1475555/150‑1 (MQ=255)
                                                                                    tAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTga                                                          <  1:1018734/150‑1 (MQ=255)
                                                                                    tAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTga                                                          >  1:106807/1‑150 (MQ=255)
                                                                                    tAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTGCCAGTTga                                                          <  1:885615/150‑1 (MQ=255)
                                                                                    tAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGAGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTga                                                          <  2:342671/150‑1 (MQ=255)
                                                                                                 aaccggTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCt                                             <  1:37146/146‑1 (MQ=255)
                                                                                                   cGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTg                                           <  1:1598206/150‑1 (MQ=255)
                                                                                                     ggTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGcc                                         <  2:1054318/150‑1 (MQ=255)
                                                                                                      gTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCa                                        <  2:574575/150‑1 (MQ=255)
                                                                                                          aCGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGt                                    >  2:1361208/1‑150 (MQ=255)
                                                                                                                      cGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTTGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGAt                        <  1:1880200/150‑1 (MQ=255)
                                                                                                                      cGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGAt                        <  1:1850336/150‑1 (MQ=255)
                                                                                                                       ggCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTg                                                                             <  2:1780208/96‑1 (MQ=255)
                                                                                                                       ggCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTg                                                                             >  1:1780208/1‑96 (MQ=255)
                                                                                                                       ggCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTGTGTTCCAGTTGAGACACCAAACGCTTGCCGGCGTCGATATTGCGATg                       >  1:998204/1‑150 (MQ=255)
                                                                                                                       ggCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGacac                                                     <  1:1323918/120‑1 (MQ=255)
                                                                                                                       ggCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGacac                                                     >  2:1323918/1‑120 (MQ=255)
                                                                                                                        gCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATgc                      >  2:1681551/1‑150 (MQ=255)
                                                                                                                           aaTTCATGCTGCTCGTCCAGTTACCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATGCGCa                   <  2:998204/150‑1 (MQ=255)
                                                                                                                                  ggTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGTTATCTGCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTAGATCTTTTGATGCACACccacca            >  1:150302/1‑150 (MQ=255)
                                                                                                                                     gaaccagcAGTTGCAGCGGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATGCGCACCCAGCATGc         >  2:554445/8‑150 (MQ=255)
                                                                                                                                       tcGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATGCGCACCCAGCATGCGt       >  1:426547/1‑150 (MQ=255)
                                                                                                                                         gTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATGCGCACCCAGCATGCGTGc     >  1:1566078/1‑150 (MQ=255)
                                                                                                                                         cagcAGTTGCAGCGGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTg                                                                             <  1:412851/75‑1 (MQ=39)
                                                                                                                                         cagcAGTTGCAGCGGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTg                                                                             >  2:412851/4‑78 (MQ=39)
                                                                                                                                           gcAGTTGCAGCGGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGGAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATGCGCACCCAGCATGCGTGCCa   >  2:1609571/2‑150 (MQ=255)
                                                                                                                                           gcAGTTGCAGCGGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATGCGCACCCAGCATGCGTGCCa   >  1:589335/2‑150 (MQ=255)
                                                                                                                                           gcAGTTGCAGCGGCAACTGCTGGTTCACCCAGGACAACTCACGGGCGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATGCGCACCCAGCATGCGTGccc   >  2:363690/2‑149 (MQ=255)
                                                                                                                                            cAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATGCGCACCCAGCATGCGTGCCAc  >  1:140422/1‑150 (MQ=255)
                                                                                                                                              gTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATGCGCACCCa               >  1:1557379/1‑135 (MQ=255)
                                                                                                                                              gTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTGCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATGCGCACCCa               <  2:1557379/135‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
TCCATGTCCATGGCAATCAGCGACAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATGCGCACCCAGCATGCGTGCCAC  >  NC_002947/4997502‑4997791

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