Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 467,955 G→T 5.1% P1001H (CCC→CAC)  PP_0386 ← sensory box protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947467,9550GT 5.1% 241.3 / 3.2 78P1001H (CCC→CAC) PP_0386sensory box protein
Reads supporting (aligned to +/- strand):  ref base G (38/36);  new base T (2/2);  total (40/38)
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.21e-01

CAGTACCCAGTCGCCGAGCTGCAGGATGGTGCCGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGCGGGTGTAGGCTTCGACGCGGTTGCGGCCCTTGGCCTTGGAGCGGTACATGGCGGC  >  NC_002947/467807‑468103
                                                                                                                                                    |                                                                                                                                                    
cAGTACCCAGTCGCCGAGCTGCAGGATGGTGCCGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTgg                                                                                                                                                      >  4:113696/1‑149 (MQ=255)
  gTACCCAGTCGCCGAGCTGCAGGATGGTGCCGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGc                                                                                                                                                    >  1:219853/1‑149 (MQ=255)
      ccAGTCGCCGAGCTGCAGGATGGTGCCGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTgg                                                                                                                                                 >  8:290610/1‑148 (MQ=255)
       cAGTCGCCGAGCTGCAGGATGGTGCCGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGaa                                                                                                                                               <  5:59441/149‑1 (MQ=255)
             ccGAGCTGCAGGATGGTGCCGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAg                                                                                                                                         <  7:132823/149‑1 (MQ=255)
             ccGAGCTGCAGGATGGTGCCGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAg                                                                                                                                         <  7:81488/149‑1 (MQ=255)
                    gCAGGATGGTGCCGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCAtt                                                                                                                                  >  5:304927/1‑149 (MQ=255)
                      aGGATGGTGCCGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTc                                                                                                                                >  3:133543/1‑149 (MQ=255)
                        gATGGTGCCGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGa                                                                                                                                                >  3:130855/1‑131 (MQ=255)
                        gATGGTGCCGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGa                                                                                                                                                >  3:125996/1‑131 (MQ=255)
                        gATGGTGCCGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGa                                                                                                                                                <  4:130855/131‑1 (MQ=255)
                        gATGGTGCCGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGa                                                                                                                                                <  4:125996/131‑1 (MQ=255)
                        gATGGTGCCGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGt                                                                                                                              >  3:199651/1‑149 (MQ=255)
                         aTGGTGCCGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGTTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCGTGGGCTGGAAGCACAGGCTCATTTCGtt                                                                                                                             <  1:20721/149‑1 (MQ=255)
                               ccGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTc                                                                                                                       <  1:175730/149‑1 (MQ=255)
                                  ttCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCCTGTGCTTCCAGCCCAAGCTCagcctgaa                                                                                                                             >  5:239345/1‑133 (MQ=38)
                                  ttCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCCTGTGCTTCCAGCCCAAGCTCagcctgaa                                                                                                                             <  6:239345/140‑8 (MQ=38)
                                       ttCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGccc                                                                                                               >  8:91938/1‑149 (MQ=255)
                                          cGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGgc                                                                                                            <  7:108785/149‑1 (MQ=255)
                                           gCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCCTGTGCTTCCAGCCCAAGCTCagcc                                                                                                                                 <  7:69484/127‑4 (MQ=255)
                                           gCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCCTGTGCTTCCAGCCCAAGCTCagcc                                                                                                                                 >  8:69484/1‑124 (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                                                                                                 >  3:61055/1‑149 (MQ=255)
                                                     tGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGc                                                                                                 <  2:262261/149‑1 (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)
                                                                            ccAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTggg                                                                                                                                                     <  2:233932/74‑1 (MQ=255)
                                                                            ccAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTggg                                                                                                                                                     >  1:233932/1‑74 (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                                                                     >  2:199574/1‑119 (MQ=255)
                                                                                                               ggCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTggg                                                                     <  1:199574/119‑1 (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:128593/1‑149 (MQ=255)
                                                                                                                             tctGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGCGGGTGTAGGCTTCGACGCGGTTGCGGCccttg                         >  5:212488/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                         >  3:117023/1‑141 (MQ=255)
                                                                                                                                     ttCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGCGGGTGTAGGCTTCGACGCGGTTGCGGCccttg                         <  4:117023/141‑1 (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)
                                                                                                                                                    |                                                                                                                                                    
CAGTACCCAGTCGCCGAGCTGCAGGATGGTGCCGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGCGGGTGTAGGCTTCGACGCGGTTGCGGCCCTTGGCCTTGGAGCGGTACATGGCGGC  >  NC_002947/467807‑468103

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