Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 467,952 T→C 5.0% K1002R (AAG→AGG)  PP_0386 ← sensory box protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947467,9520TC 5.0% 238.5 / 2.1 82K1002R (AAG→AGG) PP_0386sensory box protein
Reads supporting (aligned to +/- strand):  ref base T (41/35);  new base C (2/2);  total (43/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.99e-01

TCCAGTACCCAGTCGCCGAGCTGCAGGATGGTGCCGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGCGGGTGTAGGCTTCGACGCGGTTGCGGCCCTTGGCCTTGGAGCGGT  >  NC_002947/467805‑468093
                                                                                                                                                   |                                                                                                                                             
tCCAGTACCCAGTCGCCGAGCTGCAGGATGGTGCCGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCtt                                                                                                                                              >  7:55429/1‑149 (MQ=255)
 ccAGTACCCAGTCGCCGAGCTGCAGGATGGTGCCGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTg                                                                                                                                             >  1:65704/1‑149 (MQ=255)
  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                                                                                                                                      <  4:125996/131‑1 (MQ=255)
                          gATGGTGCCGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGa                                                                                                                                      <  4:130855/131‑1 (MQ=255)
                          gATGGTGCCGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGa                                                                                                                                      >  3:125996/1‑131 (MQ=255)
                          gATGGTGCCGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGa                                                                                                                                      >  3:130855/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                                                                                                                   >  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                                                                                  <  4:133543/149‑1 (MQ=255)
                                                            tgtTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAg                                                                                  <  6:304927/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                                                                                                                                           >  1:233932/1‑74 (MQ=255)
                                                                              ccAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTggg                                                                                                                                           <  2:233932/74‑1 (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)
                                                                                                                                                   |                                                                                                                                             
TCCAGTACCCAGTCGCCGAGCTGCAGGATGGTGCCGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGCGGGTGTAGGCTTCGACGCGGTTGCGGCCCTTGGCCTTGGAGCGGT  >  NC_002947/467805‑468093

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