Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,425,661 G→A 9.3% intergenic (+63/‑442) 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,6610GA 9.3% 205.0 / 9.7 75intergenic (+63/‑442)PP_1249/PP_1250DUF4223 domain‑containing protein/group II intron‑encoding maturase
Reads supporting (aligned to +/- strand):  ref base G (33/35);  new base A (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 = 5.39e-01

TACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGT  >  NC_002947/1425518‑1425802
                                                                                                                                               |                                                                                                                                             
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                                                                                                                            >  5:102790/1‑149 (MQ=255)
              ttCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTgcg                                                                                                                            >  3:146620/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                                                                                                                                          >  2:121969/1‑121 (MQ=255)
                            ccAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAAACCCCCGGCCATGTGCATGACCggg                                                                                                                                          <  1:121969/121‑1 (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)
                                                                                                                                               |                                                                                                                                             
TACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGT  >  NC_002947/1425518‑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: ‑