Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,988,829 A→G 24.1% intergenic (+123/‑423) 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,8290AG24.1% 105.5 / 31.4 58intergenic (+123/‑423)kefB‑III/PP_mr64glutathione‑regulated potassium/H+ antiporter/C4 ncRNA
Reads supporting (aligned to +/- strand):  ref base A (18/26);  new base G (7/7);  total (25/33)
Fisher's exact test for biased strand distribution p-value = 7.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

TTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGG  >  NC_002947/5988686‑5988974
                                                                                                                                                     |                                                                                                                                                 
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             <  1:124560‑M2/148‑29 (MQ=255)
                                                                                                                                        cTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctc             <  4:15247‑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)
                                                                                                                                                     |                                                                                                                                                 
TTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCC‑‑‑GCTCCTACGT‑‑TTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGG  >  NC_002947/5988686‑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: