Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,997,652 T→G 8.3% Q166P (CAG→CCG)  PP_4405 ← sensory box protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,997,6520TG 8.3% 238.2 / 12.4 84Q166P (CAG→CCG) PP_4405sensory box protein
Reads supporting (aligned to +/- strand):  ref base T (33/44);  new base G (5/2);  total (38/46)
Fisher's exact test for biased strand distribution p-value = 2.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

ATGTCCATGGCAATCAGCGACAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATGCGCACCCAGCATGCGTGCCACGGAACCGTC  >  NC_002947/4997505‑4997800
                                                                                                                                                   |                                                                                                                                                    
atgtccatgGCAATCAGCGACAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGc                                                                                                                                                    >  1:813678/1‑150 (MQ=255)
atgtccatgGCAATCAGCGACAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGc                                                                                                                                                    >  2:1723677/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:2198427/150‑1 (MQ=255)
                  gACAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAg                                                                                                                                  >  2:1018734/1‑150 (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:1475503/1‑150 (MQ=255)
                                                                aCGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGt                                                                                    >  1:1489981/1‑150 (MQ=255)
                                                                    gcTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGttctt                                                                                <  2:2085482/150‑1 (MQ=255)
                                                                    gcTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGttctt                                                                                <  2:173938/150‑1 (MQ=255)
                                                                    gcTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGttctt                                                                                >  2:497874/1‑150 (MQ=255)
                                                                    gcTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGttctt                                                                                <  1:2088743/150‑1 (MQ=255)
                                                                        tcCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGt                                                                            >  1:279251/1‑150 (MQ=255)
                                                                        tcCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGt                                                                            <  2:279251/150‑1 (MQ=255)
                                                                          cAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTc                                                                          >  1:108417/1‑150 (MQ=255)
                                                                          cAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTc                                                                          <  2:108417/150‑1 (MQ=255)
                                                                            gCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCa                                                                        <  1:497874/150‑1 (MQ=255)
                                                                                 tAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTga                                                                   <  2:1475555/150‑1 (MQ=255)
                                                                                 tAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTga                                                                   >  1:106807/1‑150 (MQ=255)
                                                                                 tAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTga                                                                   <  1:1018734/150‑1 (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                                                                                      >  1:1780208/1‑96 (MQ=255)
                                                                                                                    ggCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTg                                                                                      <  2:1780208/96‑1 (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                                                                                      >  2:412851/4‑78 (MQ=39)
                                                                                                                                      cagcAGTTGCAGCGGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTg                                                                                      <  1:412851/75‑1 (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)
                                                                                                                                                 gCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCtt                                                     <  2:916824/100‑1 (MQ=255)
                                                                                                                                                 gCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCtt                                                     >  1:916824/1‑100 (MQ=255)
                                                                                                                                                 gCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGc                                               <  2:106516/106‑1 (MQ=255)
                                                                                                                                                 gCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGc                                               >  1:106516/1‑106 (MQ=255)
                                                                                                                                                  cTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATGCGCACCCAGCATGCGTGCCACGGAACCGTc  <  2:1489981/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ATGTCCATGGCAATCAGCGACAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATGCGCACCCAGCATGCGTGCCACGGAACCGTC  >  NC_002947/4997505‑4997800

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