Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,425,655 C→T 9.3% intergenic (+57/‑448) PP_1249 → / → PP_1250 DUF4223 domain‑containing protein/group II intron‑encoding maturase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,425,6550CT 9.3% 206.4 / 8.8 75intergenic (+57/‑448)PP_1249/PP_1250DUF4223 domain‑containing protein/group II intron‑encoding maturase
Reads supporting (aligned to +/- strand):  ref base C (33/35);  new base T (3/4);  total (36/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 = 2.34e-01

CACCTACGACTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGT  >  NC_002947/1425508‑1425802
                                                                                                                                                   |                                                                                                                                                   
cACCTACGACTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGca                                                                                                                                                    >  3:212960/1‑149 (MQ=255)
      cGACTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATg                                                                                                                                                <  4:277880/147‑1 (MQ=255)
          tACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCggg                                                                                                                                          >  6:319710/1‑149 (MQ=255)
           aCCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCgggg                                                                                                                                         <  1:83071/149‑1 (MQ=255)
            ccTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGTCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGCCATGTGCATGACCgggg                                                                                                                                         <  5:247301/148‑1 (MQ=255)
            ccTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCggggg                                                                                                                                        >  2:59368/1‑149 (MQ=255)
                ttCCACCCTTCGGTTTCCATCTCCATGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGCCATGTGCATGACCgggggcctt                                                                                                                                    <  3:260311/149‑5 (MQ=255)
                  ccACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGttttt                                                                                                                                   <  4:188642/148‑1 (MQ=255)
                  ccACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGttttt                                                                                                                                   <  1:59368/148‑1 (MQ=255)
                   cACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTtgt                                                                                                                                 >  2:54646/1‑149 (MQ=255)
                       cTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTgc                                                                                                                             >  4:275715/1‑149 (MQ=255)
                        ttCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTgcg                                                                                                                            >  3:146620/1‑149 (MQ=255)
                        ttCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTgcg                                                                                                                            >  5:102790/1‑149 (MQ=255)
                            gTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAg                                                                                                                        >  5:306037/1‑149 (MQ=255)
                                 cATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATg                                                                                                                   >  7:240302/1‑149 (MQ=255)
                                     tcCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGgcgc                                                                                                               >  5:37634/1‑149 (MQ=255)
                                      ccAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAAACCCCCGGCCATGTGCATGACCggg                                                                                                                                          <  1:121969/121‑1 (MQ=255)
                                      ccAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAAACCCCCGGCCATGTGCATGACCggg                                                                                                                                          >  2:121969/1‑121 (MQ=255)
                                       cAAGATCATCGGTGGCGGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCATGCCTGgcgcac                                                                                                             <  3:82964/149‑3 (MQ=255)
                                        aaGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCATGCCTGgcgcaca                                                                                                            <  5:319710/149‑4 (MQ=255)
                                         aGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTg                                                                                                           <  8:240302/149‑1 (MQ=255)
                                          gATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTgc                                                                                                          <  4:212960/149‑1 (MQ=255)
                                               tcGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTa                                                                                                     <  4:146620/149‑1 (MQ=255)
                                                 ggTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAt                                                                                                    <  6:306037/148‑1 (MQ=255)
                                                         cGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAAc                                                                                             <  6:44506/147‑1 (MQ=255)
                                                              cgatcgatAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTg                                                                                      <  6:41705/149‑1 (MQ=255)
                                                               gatcgatAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTgg                                                                                     >  2:76162/1‑149 (MQ=255)
                                                                  aagattaCTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTtg                                                                                  >  7:114691/7‑149 (MQ=255)
                                                                   gatAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGttttt                                                                                                                                   <  2:20392/99‑1 (MQ=255)
                                                                   gatAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGttttt                                                                                                                                   >  1:20392/1‑99 (MQ=255)
                                                                   gatAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTtgt                                                                                 <  3:275715/149‑1 (MQ=255)
                                                                    atAAACTGCCTCAGCAGCAGTAATGTTGGCGGTACCCTGATCGCGACGGTTGGGGCGTGATCGAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGCGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTtgtg                                                                                <  3:185607/149‑1 (MQ=37)
                                                                        aCTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTtgtg                                                                                <  1:144746/145‑1 (MQ=255)
                                                                        aCTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTtgtg                                                                                >  2:144746/1‑145 (MQ=255)
                                                                             cTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTtgtg                                                                                >  2:159687/1‑140 (MQ=255)
                                                                             cTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTtgtg                                                                                <  1:159687/140‑1 (MQ=255)
                                                                                agcagcagTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTgg                                                                    >  5:238777/1‑149 (MQ=255)
                                                                                agcagcagTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTgg                                                                    <  1:54646/149‑1 (MQ=255)
                                                                                    gcagTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGa                                                                >  6:323801/1‑149 (MQ=255)
                                                                                      agTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACt                                                              >  7:250739/1‑149 (MQ=255)
                                                                                      agTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACt                                                              >  1:222838/1‑149 (MQ=255)
                                                                                        tAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTg                                                            <  8:250739/149‑1 (MQ=255)
                                                                                          aTCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGCCATGTGCATGACCg                                                                                                                                            <  4:42769/67‑1 (MQ=37)
                                                                                               ggCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAg                                                         <  6:102790/145‑1 (MQ=255)
                                                                                                gCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGa                                                     <  1:144835/148‑1 (MQ=255)
                                                                                                 cGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGaaa                                                   <  8:114691/149‑1 (MQ=255)
                                                                                                    tACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTtgtgt                                                                                                                               <  3:14819/70‑1 (MQ=255)
                                                                                                    tACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTtgtgt                                                                                                                               >  4:14819/1‑70 (MQ=255)
                                                                                                     aCCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTcc                                               >  4:36708/1‑149 (MQ=255)
                                                                                                      cccTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCct                                              <  6:37634/149‑1 (MQ=255)
                                                                                                      cccTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCct                                              <  6:238777/149‑1 (MQ=255)
                                                                                                      cccTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCct                                              <  5:323801/149‑1 (MQ=255)
                                                                                                       ccTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCctc                                             >  4:201953/1‑149 (MQ=255)
                                                                                                        cTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCctc                                             >  2:222021/1‑148 (MQ=255)
                                                                                                            tCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTacac                                        <  3:36708/149‑1 (MQ=255)
                                                                                                               cgACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTgc                                                                                                          <  4:175443/80‑1 (MQ=255)
                                                                                                               cgACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTgc                                                                                                          >  3:175443/1‑80 (MQ=255)
                                                                                                                gACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAg                                                         >  8:309027/1‑128 (MQ=255)
                                                                                                                gACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAg                                                         <  7:309027/128‑1 (MQ=255)
                                                                                                                gACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGACAAGTGGGTGACTTAGAGGTGAAAGTCCTCTacacaca                                     >  4:82239/1‑147 (MQ=255)
                                                                                                                   gCTTGCGGCGTTATCCAAGGCCCCCGGTCATGCACATGGACGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGc                                 <  3:82239/149‑1 (MQ=255)
                                                                                                                    cTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGc                                 >  7:2523/1‑148 (MQ=255)
                                                                                                                     ttGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGTTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACACGGc                                 >  5:65741/1‑147 (MQ=255)
                                                                                                                          gcgTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAgggg                           >  8:42236/1‑148 (MQ=255)
                                                                                                                            gTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCACCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGGGACTTGGAGGGGAAAGTCCTCTACACACCCGGCAAGGGGAAg                        <  6:285930/149‑1 (MQ=255)
                                                                                                                                tCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTc                                                <  1:315293/121‑1 (MQ=255)
                                                                                                                                tCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTc                                                >  2:315293/1‑121 (MQ=255)
                                                                                                                                 ccAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAgtgt                     >  3:151914/1‑147 (MQ=255)
                                                                                                                                 ccAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAgtgt                     >  5:333066/1‑147 (MQ=255)
                                                                                                                                        ccccGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCctc                                             >  4:266434/1‑116 (MQ=255)
                                                                                                                                        ccccGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCctc                                             <  3:266434/116‑1 (MQ=255)
                                                                                                                                         cccGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAgg           >  8:291906/1‑149 (MQ=255)
                                                                                                                                         cccGGCCATGTGCATGACCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAgg           >  2:132981/1‑149 (MQ=255)
                                                                                                                                          ccGGCCATGTGCATGACCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAgg           >  4:65545/1‑148 (MQ=255)
                                                                                                                                                aTGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAACGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGgt    <  1:76162/149‑1 (MQ=255)
                                                                                                                                                aTGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGgt    <  2:222838/149‑1 (MQ=255)
                                                                                                                                                   cacaTGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGgtgt  <  4:151914/148‑1 (MQ=39)
                                                                                                                                                   |                                                                                                                                                   
CACCTACGACTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGT  >  NC_002947/1425508‑1425802

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: