Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,609,858 C→A 6.6% K125N (AAG→AAT) ‡ PP_3183 ← SCO1/SenC family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,609,8580CA 6.6% 190.1 / 3.7 61K125N (AAG→AAT) ‡PP_3183SCO1/SenC family protein
Reads supporting (aligned to +/- strand):  ref base C (27/30);  new base A (2/2);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.28e-01

CCTGGTTGCCGATGATGAGGCTCAGGTTGTGGTCTTTGGAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCGGTCACCGAGGATTTTCTGTACCTGGCGTAGGCGCGCAGTTTCCACTGGGCAGGAATCGGAG  >  NC_002947/3609713‑3610005
                                                                                                                                                 |                                                                                                                                                   
ccTGGTTGCCGATGATGAGGCTCAGGTTGTGGTCTTTGGAGCGCCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGCCCAGCCTGGGCCGATGCCGAACttt                                                                                                                                                  >  2:23621/1‑149 (MQ=255)
 cTGGTTGCCGATGATGAGGCTCAGGTTGTGGTCTTTGGAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACtttt                                                                                                                                                 <  2:171468/149‑1 (MQ=255)
    gTTGCCGATGATGAGGCTCAGGTTGTGGTCTTTGGAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCgg                                                                                                                                              >  4:361737/1‑149 (MQ=255)
        ccGATGATGAGGCTCAGGTTGTGGTCTTTGGAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATa                                                                                                                                          <  6:141929/149‑1 (MQ=255)
        ccGATGATGAGGCTCAGGTTGTGGTCTTTGGAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATa                                                                                                                                          >  4:200909/1‑149 (MQ=255)
        ccGATGATGAGGCTCAGGTTGTGGTCTTTGGAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAAAAGTTCGGCAt                                                                                                                                           <  4:179965/148‑1 (MQ=255)
         cGATGATGAGGCTCAGGTTGTGGTCTTTGGAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAAAAGTTCGGCAt                                                                                                                                           <  6:85192/147‑1 (MQ=255)
         cGATGATGAGGCTCAGGTTGTGGTCTTTGGAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAAAAGTTCGGCAt                                                                                                                                           >  5:85192/1‑147 (MQ=255)
                gaGGCTCAGGTTGTGGTCTTTGGAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGttt                                                                                                                                    >  4:129895/1‑147 (MQ=255)
                 aGGCTCAGGTTGTGGTCTTTGGAGCGCCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGCCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGa                                                                                                                                  <  1:23621/148‑1 (MQ=255)
                  ggCTCAGGTTGTGGTCTTTGGAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAgg                                                                                                                                <  3:129895/149‑1 (MQ=255)
                     tCAGGTTGTGGTCTTTGGAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAgg                                                                                                                                <  7:175799/146‑1 (MQ=255)
                     tCAGGTTGTGGTCTTTGGAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAgg                                                                                                                                >  8:175799/1‑146 (MQ=255)
                       aGGTTGTGGTCTTTGGAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGc                                                                                                                           >  3:364783/1‑149 (MQ=255)
                             ttgTCTTTGGAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGgt                                                                                                                     <  4:136687/147‑1 (MQ=255)
                             tgGTCTGTGGAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGgt                                                                                                                     <  4:131016/149‑1 (MQ=255)
                                tCTTTGGAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTc                                                                                                                  <  7:24344/149‑1 (MQ=255)
                                tCTTTGGAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTc                                                                                                                  <  1:19227/149‑1 (MQ=255)
                                tCTTTGGAGCGGCCGTGTTCCAGGCCATCGATCCACAGCCGCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTc                                                                                                                  <  1:13151/149‑1 (MQ=255)
                                 cTTTGGAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCg                                                                                                                 <  2:112818/149‑1 (MQ=255)
                                  tttGGAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGt                                                                                                                >  7:212393/1‑149 (MQ=255)
                                  tttGGAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGt                                                                                                                >  2:129075/1‑149 (MQ=255)
                                      gAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTtat                                                                                                             <  3:200909/148‑1 (MQ=255)
                                             cGTTCTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCgggg                                                                                                                       <  6:250091/131‑1 (MQ=255)
                                             cGTTCTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCgggg                                                                                                                       >  5:250091/1‑131 (MQ=255)
                                                tttCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGAt                                                                                                   <  1:270486/148‑1 (MQ=255)
                                                tttCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGAt                                                                                                   >  2:270486/1‑148 (MQ=255)
                                                    cAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATTGGGATCGATGCTGa                                                                                              >  4:276297/1‑149 (MQ=255)
                                                     aGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGAt                                                                                             <  4:280091/149‑1 (MQ=255)
                                                          aTCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGa                                                                                          <  2:201054/147‑1 (MQ=255)
                                                          aTCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGa                                                                                          >  1:201054/1‑147 (MQ=255)
                                                          aTCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGGCGTTATAGGGCTCGAGGCTGATg                                                                                            >  8:250816/1‑145 (MQ=255)
                                                             gATCCACATCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTaga                                                                                     <  3:361737/149‑1 (MQ=255)
                                                              aTCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCTGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTagag                                                                                    >  8:64230/1‑149 (MQ=255)
                                                              aTCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCTGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTagag                                                                                    >  2:121587/1‑149 (MQ=255)
                                                                  acaGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAgg                                                                                   >  1:88645/1‑146 (MQ=255)
                                                                        ccAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCt                                                                          >  3:168653/1‑149 (MQ=255)
                                                                             cTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGcc                                                                     <  1:129075/149‑1 (MQ=255)
                                                                                 ggcgTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGAcc                                                                 >  8:64361/1‑149 (MQ=255)
                                                                                  gcgTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGAccc                                                                <  8:212393/149‑1 (MQ=255)
                                                                                       gTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCGGTcc                                                           >  4:148805/1‑148 (MQ=255)
                                                                                           cTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCGGTCACCGa                                                       <  7:99061/149‑1 (MQ=255)
                                                                                             cGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCGGTCACCGAg                                                      >  2:293195/1‑148 (MQ=255)
                                                                                              gATGTCGTCGGGTTCGCCGGGCAGCAAGGTCCAGCCTGGGCAGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCGGTCACCGAgg                                                     <  5:292533/148‑1 (MQ=255)
                                                                                                               cGGTCAGCAAGGTCCAGCCTGGGCTGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCttt                                                                        <  1:335459/112‑1 (MQ=255)
                                                                                                               cGGTCAGCAAGGTCCAGCCTGGGCTGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCttt                                                                        >  2:335459/1‑112 (MQ=255)
                                                                                                                      cAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCGGTCACCGAGGATTTTCTGTACCTGGCGTAGgcgcg                            >  7:104401/1‑149 (MQ=255)
                                                                                                                      cAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCGGTCACCGAGGATTTTCTGTACCTGGCGTAGgcgc                             >  8:77529/1‑148 (MQ=255)
                                                                                                                          gTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCGGTCACCGAGGATTTTCTGTACCTGGCGTAgg                                <  6:111129/141‑1 (MQ=255)
                                                                                                                          gTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCGGTCACCGAGGATTTTCTGTACCTGGCGTAgg                                >  5:111129/1‑141 (MQ=255)
                                                                                                                              aGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCGGTCACCGAGGATTTTCTGTACCTGGCGTAGGCGCGCAGTTTc                     >  7:349252/1‑148 (MQ=255)
                                                                                                                                ccTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCTGTCACCGAGGATTTTCTGTACCTGGCGTAGGCGCGCAGTTTCCAc                  <  1:293195/149‑1 (MQ=255)
                                                                                                                                ccTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCGGTCACCGAGGATTTTCTGTACCTGGCGTAGGCGCGCAGTTTc                     >  1:196015/1‑146 (MQ=255)
                                                                                                                                ccTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCGGTCACCGAGGATTTTCTGTACCTGGCGTAGGCGCGCAGTTTc                     <  2:196015/146‑1 (MQ=255)
                                                                                                                                    ggCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCGGTCACCGAGGATTTTCTGTACCTGGCGTAGGCGCGCAGTTTCCACTg                <  8:104401/147‑1 (MQ=255)
                                                                                                                                       cGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCGGTCACCGAGGATTTTCTGTACCTGGCGTAGGCGCGCAGTTTCCACTGGGCAg           <  7:64361/149‑1 (MQ=255)
                                                                                                                                       cGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCGGTCACCGAGGATTTTCTGTACCTGGCGTAGGCGCGCAGTTTCCACTGGGCAg           <  4:364783/149‑1 (MQ=255)
                                                                                                                                        taTGCCGAAAAGTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCGGTCACCGAGGATTTTCTGTACCTGGCGTAGGCGCGCAGTTTCCACTGGGCAg           >  6:78077/2‑148 (MQ=255)
                                                                                                                                              gAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATTGAGTAGAGGAAGATGTCTTTGCCGACCCGGTCACCGAGGATTTTCTGTACCTGGCGTAg                                 <  3:32225/120‑1 (MQ=255)
                                                                                                                                              gAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCAGGTCACCGAGGATTTTCTGTACCTGGCGTAg                                 >  4:32225/1‑120 (MQ=255)
                                                                                                                                                aCTTTTCGTCATAGCGTTTGAGGGTGGCGGGGGTGGCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACACGGTCACCGAGGATTTTCTGTACCTGGCGTAGGCGCGCAGTTTCCACTGGGCAGGAATCGGag  <  7:250816/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CCTGGTTGCCGATGATGAGGCTCAGGTTGTGGTCTTTGGAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCAGCAAGGTCCAGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCGGTCACCGAGGATTTTCTGTACCTGGCGTAGGCGCGCAGTTTCCACTGGGCAGGAATCGGAG  >  NC_002947/3609713‑3610005

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