Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,988,826 C→T 24.4% intergenic (+120/‑426) kefB‑III → / → PP_mr64 glutathione‑regulated potassium/H+ antiporter/C4 ncRNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,988,8260CT24.4% 95.0 / 27.5 55intergenic (+120/‑426)kefB‑III/PP_mr64glutathione‑regulated potassium/H+ antiporter/C4 ncRNA
Reads supporting (aligned to +/- strand):  ref base C (15/25);  new base T (6/7);  total (23/32)
Fisher's exact test for biased strand distribution p-value = 7.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.26e-01

CGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGG  >  NC_002947/5988680‑5988974
                                                                                                                                                        |                                                                                                                                                    
cGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAGTTATTCCTGTGGGAcgcg                                                                                                                                                     >  6:116761/1‑148 (MQ=255)
      tttGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGa                                                                                                                                              <  8:142428/149‑1 (MQ=255)
        tGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAg                                                                                                                                             <  5:116761/148‑1 (MQ=255)
         gATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCC‑‑‑GCTCATACGT‑‑TTGTTTCGGGCAAGTTATTCCTGTGGGACGAGCACGCGAGcc                                                                                                                                           >  8:103797/1‑149 (MQ=255)
                    tttcagcgcctgggaggtgttcgccgggagagtttcagcgcctgggagatcgagcgccgcccgcgCGGCGCTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGcat                                                                                                                                >  6:95135‑M2/66‑148 (MQ=255)
                       cagcgcctgggaggtgttcggcgggagagtttcagcgcctgggagatcgagggccgcccgcgCGGCGCTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGggcatgt                                                                                                                              <  5:52476‑M2/86‑4 (MQ=255)
                       cagcgcctgggaggtgttcgccgggagagtttcagcgcctgggagatcgagcgccgcccgcgCGGCGCTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGgggcatgtc                                                                                                                             >  7:215317‑M2/63‑145 (MQ=255)
                        tGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCt                                                                                                                            >  3:77571/1‑149 (MQ=255)
                        agcgcctgggaggtgttcgccgggagagtttcagcgcctgggagatcgagcgccgcccgcgCGGCGCTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGgggcatgtc                                                                                                                             <  8:215317‑M2/87‑5 (MQ=255)
                          aGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTc                                                                                                                           >  7:72836/1‑148 (MQ=255)
                          aGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTc                                                                                                                           >  3:48643/1‑148 (MQ=255)
                           ggCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGa                                                                                                                         <  4:76036/149‑1 (MQ=255)
                           ggCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGACAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGa                                                                                                                         >  2:185356/1‑149 (MQ=255)
                                 gaggtgttcgccgggagagtttcagcgcctgggagatcgagcgccgcccgcgCGGCGCTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATcgc                                                                                                                   <  5:95135‑M2/97‑1 (MQ=255)
                                  gggccatcgccgggagagtttcagcgcctgggagatcgagcgccgcccgcgCGGCGCTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATcgcg                                                                                                                  <  4:33330‑M2/98‑1 (MQ=255)
                                    gtgttcggcgggagagtttcagcgcctgggagatcgagcgccgcccgcgCGGCGCTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGcgtc                                                                                                                <  6:102326‑M2/100‑1 (MQ=255)
                                    gtgttcgccgggagagtttcagcgcctgggagagcgagcgccgcccgcgCGGCGCTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGcgt                                                                                                                 <  1:263808‑M2/99‑1 (MQ=255)
                                    gcgcGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGcgtc                                                                                                                <  4:77571/149‑1 (MQ=255)
                                     cgcgTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGcgtcg                                                                                                               >  3:258922/1‑149 (MQ=255)
                                      gcgTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGgc                                                                                                                                    >  2:258125/1‑127 (MQ=255)
                                      gcgTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGgc                                                                                                                                    <  1:258125/127‑1 (MQ=255)
                                          gCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTaca                                                                                                           >  7:33643/1‑148 (MQ=255)
                                           ccgggagagtttcagcgcctgggagatcgagcgccgcccgcgCGGCGCTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTacaaa                                                                                                         <  7:37460‑M2/107‑1 (MQ=255)
                                               gAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTacaaaca                                                                                                       <  8:72836/147‑1 (MQ=255)
                                                 gggTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGc                                                                                                   <  3:134500/149‑1 (MQ=255)
                                                  ggTTGTAGAGCGTGTGAGATTGTGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCCCATCTCCTACGTCCTCATTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTacaaa                                                                                                         <  1:84046/148‑1 (MQ=25)
                                                     tGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTc                                                                                               <  6:75656/149‑1 (MQ=255)
                                                      gTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTc                                                                                               <  4:258922/148‑1 (MQ=255)
                                                         gTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATc                                                                                                                     >  4:174209/1‑123 (MQ=255)
                                                         gTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTcgcgc                                                                                           >  1:31127/1‑149 (MQ=255)
                                                                    aTTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGcac                                                                                <  4:48643/149‑1 (MQ=255)
                                                                         gcgccgcccgcgCGGCGCTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGcgtc                                                                                                                >  6:121405‑M2/13‑112 (MQ=255)
                                                                         gcgccgcccgcgCGGCGCTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGcgtc                                                                                                                <  5:121405‑M2/100‑1 (MQ=255)
                                                                         gcgcCGCCCGCGCGGCGCTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCaggag                                                                           >  5:247235/1‑149 (MQ=14)
                                                                                ccgcgCGGCGCTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAg                                                                                                     >  4:249155‑M2/6‑116 (MQ=255)
                                                                                ccgcgCGGCGCTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAg                                                                                                     <  3:249155‑M2/111‑1 (MQ=255)
                                                                                 cgcgcgGCGCTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGTCTGGc                                                                   >  3:3442/1‑149 (MQ=14)
                                                                                 cgcgCGGCGATTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACTCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGc                                                                   <  1:185356‑M2/145‑1 (MQ=255)
                                                                                       gcgcTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGaaa                                                             <  7:103797/149‑1 (MQ=255)
                                                                                           ttCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGaaacaaa                                                         <  8:33643/149‑1 (MQ=255)
                                                                                           ttCGCGGGTAAACCC‑‑‑GCTCCTACCT‑‑TTGTTTCGGGCCAGTTATTACTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTAGCGCGCGCCTGGCACAGGAGGGATTGGCCAGaaacaaa                                                         >  6:240718/1‑149 (MQ=255)
                                                                                             cgcgGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGaaacaaac                                                        <  2:53419/148‑1 (MQ=255)
                                                                                               cgGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTa                                                     <  2:31127/149‑1 (MQ=255)
                                                                                                      aCCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCgg                                              >  6:245598/1‑149 (MQ=255)
                                                                                                         c‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGtt                                           <  6:72185/149‑1 (MQ=255)
                                                                                                                 cTACGT‑‑TTGTTTCGGGCCAATTATTCCTGCGGGATTAGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGaaacaaac                                                        >  3:201883/1‑130 (MQ=2)
                                                                                                                     gt‑‑ttgtttCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgcg                                   <  8:64077‑M2/148‑7 (MQ=255)
                                                                                                                           tgtttCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGTACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgcgaagc                               <  8:124034‑M2/149‑11 (MQ=255)
                                                                                                                             tttCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGAACTGGCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGcgcc                            >  7:139244/1‑149 (MQ=14)
                                                                                                                               tCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACACGGCGGTCGCGCGCGCCTGGCACTGGAGGGATTGGCCATAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgcc                            >  8:56030‑M2/1‑135 (MQ=255)
                                                                                                                                   gCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcg                       <  3:24111‑M2/149‑19 (MQ=255)
                                                                                                                                   gCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCgcgg                      >  3:104530/1‑149 (MQ=14)
                                                                                                                                              cTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctc             <  4:15247‑M2/148‑29 (MQ=255)
                                                                                                                                              cTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctc             <  1:124560‑M2/148‑29 (MQ=255)
                                                                                                                                                 tgGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcga           >  4:24111‑M2/1‑117 (MQ=255)
                                                                                                                                                   ggACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAATCAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatct        <  2:191966‑M2/148‑34 (MQ=255)
                                                                                                                                                      cgcgCACGCGAGCCCCCTGGCGGACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATg    <  3:108838/149‑1 (MQ=255)
                                                                                                                                                       gcgcACGCGAGCCCCCTGGCGGACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATg    >  4:108838/1‑148 (MQ=255)
                                                                                                                                                        cgcACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgcccaatctcacagg  >  7:64077‑M2/1‑110 (MQ=255)
                                                                                                                                                        |                                                                                                                                                    
CGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGG  >  NC_002947/5988680‑5988974

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: