Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 467,959 G→T 5.6% Q1000K (CAG→AAG)  PP_0386 ← sensory box protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947467,9590GT 5.6% 252.9 / 3.5 81Q1000K (CAG→AAG) PP_0386sensory box protein
Reads supporting (aligned to +/- strand):  ref base G (39/37);  new base T (2/3);  total (41/40)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.11e-01

CCAGTCGCCGAGCTGCAGGATGGTGCCGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGCGGGTGTAGGCTTCGACGCGGTTGCGGCCCTTGGCCTTGGAGCGGTACATGGCGGCG  >  NC_002947/467813‑468104
                                                                                                                                                  |                                                                                                                                                 
ccAGTCGCCGAGCTGCAGGATGGTGCCGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTgg                                                                                                                                                  >  8:290610/1‑148 (MQ=255)
 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: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                                                                                                                              <  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                                                                                                  >  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)
                                                                            gTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGc                                                                     >  5:123683/1‑149 (MQ=255)
                                                                               ggggggCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGGCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGt                                                                  <  7:219736/149‑1 (MQ=255)
                                                                                ggTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTg                                                                 <  2:65704/149‑1 (MQ=255)
                                                                                ggTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTCGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTg                                                                 <  7:277651/149‑1 (MQ=255)
                                                                                ggTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTCGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTg                                                                 <  7:281869/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        <  6:103549/142‑1 (MQ=255)
                                                                                                                                                cTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGCGGGTGTAGGCTTCGACGCGGTTGCGGCCCTTGGCCTTGGAGCGGTACATg        >  5:103549/1‑142 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CCAGTCGCCGAGCTGCAGGATGGTGCCGTTCTCTTCCGCCAAGTGGATGAAGTGTTCCGGCGGTACCTCGCCAAAGGTCGGGTGGCTCCAGCGGATCAGCGCTTCGGCACCCACCAGGCTCTGGGTCTTCAGGCTGAGCTTGGGCTGGAAGCACAGGCTCATTTCGTTGCGTTCGACGGCCCGGCGCAGCTCATGCTCCAGGGCGATGCGCTCGCTGGCCTGGGCGGTGAGGTCGCGGGTGTAGGCTTCGACGCGGTTGCGGCCCTTGGCCTTGGAGCGGTACATGGCGGCG  >  NC_002947/467813‑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: ‑