Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 467,961 A→C 5.3% F999C (TTC→TGC)  PP_0386 ← sensory box protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947467,9610AC 5.3% 242.0 / 3.1 76F999C (TTC→TGC) PP_0386sensory box protein
Reads supporting (aligned to +/- strand):  ref base A (37/35);  new base C (2/2);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

CAGTCGCCGAGCTGCAGGATGGTGCCGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGCGGGTGTAGGCTTCGACGCGGTTGCGGCCCTTGGCCTTGGAGCGGTACATGGCGGCG  >  NC_002947/467814‑468104
                                                                                                                                                   |                                                                                                                                               
cAGTCGCCGAGCTGCAGGATGGTGCCGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGaa                                                                                                                                                <  5:59441/149‑1 (MQ=255)
      ccGAGCTGCAGGATGGTGCCGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAg                                                                                                                                          <  7:81488/149‑1 (MQ=255)
      ccGAGCTGCAGGATGGTGCCGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAg                                                                                                                                          <  7:132823/149‑1 (MQ=255)
             gCAGGATGGTGCCGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCAtt                                                                                                                                   >  5:304927/1‑149 (MQ=255)
               aGGATGGTGCCGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTc                                                                                                                                 >  3:133543/1‑149 (MQ=255)
                 gATGGTGCCGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGa                                                                                                                                                 <  4:130855/131‑1 (MQ=255)
                 gATGGTGCCGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGa                                                                                                                                                 <  4:125996/131‑1 (MQ=255)
                 gATGGTGCCGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGa                                                                                                                                                 >  3:130855/1‑131 (MQ=255)
                 gATGGTGCCGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGa                                                                                                                                                 >  3:125996/1‑131 (MQ=255)
                 gATGGTGCCGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGt                                                                                                                               >  3:199651/1‑149 (MQ=255)
                  aTGGTGCCGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGTTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCGTGGGCTGGAAGCACAGGCTCATTTCGtt                                                                                                                              <  1:20721/149‑1 (MQ=255)
                        ccGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTc                                                                                                                        <  1:175730/149‑1 (MQ=255)
                           ttCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCCTGTGCTTCCAGCCCAAGCTCagcctgaa                                                                                                                              <  6:239345/140‑8 (MQ=38)
                           ttCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCCTGTGCTTCCAGCCCAAGCTCagcctgaa                                                                                                                              >  5:239345/1‑133 (MQ=38)
                                ttCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGccc                                                                                                                >  8:91938/1‑149 (MQ=255)
                                   cGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGgc                                                                                                             <  7:108785/149‑1 (MQ=255)
                                    gCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCCTGTGCTTCCAGCCCAAGCTCagcc                                                                                                                                  >  8:69484/1‑124 (MQ=255)
                                    gCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCCTGTGCTTCCAGCCCAAGCTCagcc                                                                                                                                  <  7:69484/127‑4 (MQ=255)
                                       aaGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGtt                                                                                                                              <  8:127808/128‑1 (MQ=255)
                                       aaGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGtt                                                                                                                              >  7:127808/1‑128 (MQ=255)
                                         gTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCt                                                                                                       >  4:102031/1‑149 (MQ=255)
                                          tGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCt                                                                                                       >  1:126610/1‑148 (MQ=255)
                                            gATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCAt                                                                                                    >  7:24300/1‑149 (MQ=255)
                                              tGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGc                                                                                                  <  2:262261/149‑1 (MQ=255)
                                              tGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGc                                                                                                  >  3:61055/1‑149 (MQ=255)
                                                aacggtTCCGGCGGTACCTCCCCAAAGGTCGGGGGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTGCGACGGCCCGGCGCAGCTCATGc                                                                                                  <  7:273986/143‑1 (MQ=255)
                                                aaGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTc                                                                                                <  5:114787/149‑1 (MQ=255)
                                                 aGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTCGCGTTCGACGGACCGGCGCAGCTCAttatc                                                                                                >  8:277651/1‑144 (MQ=255)
                                                 aGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTAGCGTTCGAAGGCCCGGCGCAGCTCATGCta                                                                                                >  8:281869/1‑147 (MQ=255)
                                                   tgtTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAg                                                                                             <  6:304927/149‑1 (MQ=255)
                                                   tgtTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAg                                                                                             <  4:133543/149‑1 (MQ=255)
                                                     ttCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAggg                                                                                           <  2:219853/149‑1 (MQ=255)
                                                      tCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGc                                                                                          >  4:289889/1‑149 (MQ=255)
                                                       ccgccgGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCAACCACCATGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGc                                                                                          >  5:225245/1‑148 (MQ=255)
                                                           cggTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGTGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATgc                                                                                     <  6:74778/149‑1 (MQ=255)
                                                             gTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATgcgc                                                                                   <  1:1489/149‑1 (MQ=255)
                                                                ccTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGcgctc                                                                                 <  7:290610/148‑1 (MQ=255)
                                                                 cTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCCCCCCCCAGGCTCGGGGGCGTCCGGCGGAGCTTGGGCTGGAAGCACAGGCGCAGTGCGGGGCGTTCGACGGCCCGGCGCCGCTCATCCTCCAGGGCGAGGcgctc                                                                                 >  8:285605/1‑147 (MQ=255)
                                                                    gCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTgg                                                                            >  4:144342/1‑149 (MQ=255)
                                                                           gTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGc                                                                     >  5:123683/1‑149 (MQ=255)
                                                                              ggggggCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGGCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGt                                                                  <  7:219736/149‑1 (MQ=255)
                                                                               ggTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTg                                                                 <  2:65704/149‑1 (MQ=255)
                                                                               ggTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTCGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTg                                                                 <  7:281869/149‑1 (MQ=255)
                                                                               ggTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTCGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTg                                                                 <  7:277651/149‑1 (MQ=255)
                                                                                   gCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGt                                                             >  5:65964/1‑149 (MQ=255)
                                                                                       cAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTcgcg                                                         >  8:58361/1‑149 (MQ=255)
                                                                                        aGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGCgg                                                        <  3:289889/149‑1 (MQ=255)
                                                                                                 gcgcTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGCGGGTGTAGGCt                                               >  7:212088/1‑149 (MQ=255)
                                                                                                     ttCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGCGGGTGTAGGCTTCGa                                           >  1:197363/1‑149 (MQ=255)
                                                                                                        ggCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTggg                                                                      <  1:199574/119‑1 (MQ=255)
                                                                                                        ggCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTggg                                                                      >  2:199574/1‑119 (MQ=255)
                                                                                                         gCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGCGGGTGTAGGCTTCGAcgcg                                       >  5:70715/1‑149 (MQ=255)
                                                                                                            cccaccAGGCTCTGGGTCTTCAGGATGAGCGTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTGCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGCGGGTGTAGGCTTCGACGCGGtt                                    <  2:219887/149‑1 (MQ=255)
                                                                                                              caccaGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGCGGGTGTAGGCTTCGACGCGGTTGc                                  >  4:62009/1‑149 (MQ=255)
                                                                                                               accaGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGCGGGTGTAGGCTTCGAcgccgccat                                  >  1:42251/1‑142 (MQ=255)
                                                                                                                ccaGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGCGGGTGTAGGCTTCGACGCGGTTGCgg                                <  4:61055/149‑1 (MQ=255)
                                                                                                                      tctGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGCGGGTGTAGGCTTCGACGCGGTTGCGGCccttg                          >  5:212488/1‑149 (MQ=255)
                                                                                                                      tctGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGCGGGTGTAGGCTTCGACGCGGTTGCGGCccttg                          >  5:128593/1‑149 (MQ=255)
                                                                                                                       ctGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGCGGGTGTAGGCTTCGACGCGGTTGCGGCccttgg                         >  3:311292/1‑149 (MQ=255)
                                                                                                                        tGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGTGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGAGGGTGTAGCCTTCGACGCGGTTGCGGCccttggc                        >  7:229827/1‑149 (MQ=255)
                                                                                                                        tGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGCGGGTGTAGGCTTCGACGCGGTTGCGGCccttg                          <  8:55429/147‑1 (MQ=255)
                                                                                                                        tGGGGCTGCCGGCGGAGCTTGGGCTGGACTCACAGGCTCATTTCGTGGCGTTCGCCGGCCCGGCGCAGCTCATGCTCCAGGGCGAGGCGCTCGCTGGCCGGGGCGGGGAGGTCGCGGGTGTAGGCGTCGACGCGGTTGCGGCccttggc                        <  7:285605/149‑1 (MQ=255)
                                                                                                                          ggTCTTCAGGCTGAGCTTGGGCTGGAAGCAGAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGCGGGTGTAGGCTTCGACGCGGTTGCGGCccttggcc                       <  8:212088/148‑1 (MQ=255)
                                                                                                                          ggTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGATGTCGCGGGTGTAGGCTTCGACGCGGTTGCGGCccttggc                        <  8:24300/147‑1 (MQ=255)
                                                                                                                              ttCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGCGGGTGTAGGCTTCGACGCGGTTGCGGCccttggccttg                    >  5:198144/1‑147 (MQ=255)
                                                                                                                              ttCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGCGGGTGTAGGCTTCGACGCGGTTGCGGCccttggccttg                    >  7:56744/1‑147 (MQ=255)
                                                                                                                              ttCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGCGGGTGTAGGCTTCGACGCGGTTGCGGCccttg                          <  4:117023/141‑1 (MQ=255)
                                                                                                                              ttCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGCGGGTGTAGGCTTCGACGCGGTTGCGGCccttg                          >  3:117023/1‑141 (MQ=255)
                                                                                                                                cAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGCGGGTGTAGGCTTCGACGCGGTTGCGGCccttggccttgg                   <  7:58361/146‑1 (MQ=255)
                                                                                                                                 aGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGCGGGTGTAGGCTTCGACGCGGTTGCGGCCCTTGGCCTTGGAGc                <  3:144342/148‑1 (MQ=255)
                                                                                                                                  ggCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGCGGGTGTAGGCTTCGACGCGGTTGCGGCccttggcct                      >  8:116097/1‑141 (MQ=255)
                                                                                                                                  ggCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGCGGGTGTAGGCTTCGACGCGGTTGCGGCccttggcct                      <  7:116097/141‑1 (MQ=255)
                                                                                                                                  ggCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGCGGGTGTAGGCTTCGACGCGGTTGCGGCCCTTGGCCTTGGAGCgg              >  3:147088/1‑149 (MQ=255)
                                                                                                                                  ggCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGATCCAGGGCGATGCGCTAGCTGGCCAGGGCGGTGAGGTCGCGGGTGTAGGCTTCGACGCGGTTGCGGCCCTTGGCCTTGGAGCgg              >  1:260172/1‑149 (MQ=255)
                                                                                                                                   gCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGCGGGTGTAGGCTTCGACGCGGTTGCGGCCCTTGGCCTTGGAGCGGt             <  6:128593/149‑1 (MQ=255)
                                                                                                                                   gCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCGCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGCGGGTGTAGGCTTCGACGCGGTTGCGGCCCTTGGCCTTGGAGCGGt             <  4:199651/149‑1 (MQ=255)
                                                                                                                                            gggCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATTCGCTCGCTGTCATGGGCGGTGAGGAAGCGTGTGTAGCCTTCGACGCGGTTGCGGCCATTGGCCTTGGAGCGGTACATggcgg    >  5:327727/1‑149 (MQ=255)
                                                                                                                                             ggCTTGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGCGGGTGTAGGCTTCGACGCGGTTGCGGCCCTTGGCCTTGGAGCGGTACATggcggc   <  6:123683/149‑1 (MQ=255)
                                                                                                                                              gCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGCGGGTGTAGGCTTCGACGCGGTTGCGGCCCTTGGCCTTGGAGCGGTACATGgcggcg  >  5:67761/1‑149 (MQ=255)
                                                                                                                                               cTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGCGGGTGTAGGCTTCGACGCGGTTGCGGCCCTTGGCCTTGGAGCGGTACATg        >  5:103549/1‑142 (MQ=255)
                                                                                                                                               cTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGCGGGTGTAGGCTTCGACGCGGTTGCGGCCCTTGGCCTTGGAGCGGTACATg        <  6:103549/142‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
CAGTCGCCGAGCTGCAGGATGGTGCCGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGCGGGTGTAGGCTTCGACGCGGTTGCGGCCCTTGGCCTTGGAGCGGTACATGGCGGCG  >  NC_002947/467814‑468104

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: