Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,019,947	G→A	6.8%	G55G (GGC→GGT)	PP_4423 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,019,947	0	G	A	6.8%	174.9 / 2.1	59	G55G (GGC→GGT)	PP_4423	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (25/30);  new base A (2/2);  total (27/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 = 2.71e-01

TCGGAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCACCCCACGCAGAGTCATCACGTGAGTAAGGCAGCTTCAAGAACCCGTGCGACACCCCTTGTTGGTCGAAATCGATAGTTGGCTT  >  NC_002947/5019802‑5020093                                                                                                                                                   |                                                                                                                                                   tCGGAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTgccgc                                                                                                                                                   >  4:170574/1‑147 (MQ=255)   ggAGAGGTTCTGGAGGCTTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGt                                                                                                                                               <  3:170574/149‑1 (MQ=255)      gagGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTgag                                                                                                                                            <  1:11012/149‑1 (MQ=255)        ggTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTgagga                                                                                                                                          <  3:273187/149‑1 (MQ=255)        ggTGCTGGAGGGGTTTTTAAAAGCAGCTTAGTTCATGTAGGGCAAACTAATGATTGTCCCTCTGACTGTATCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGACGCCAGTgagga                                                                                                                                          <  6:4353/149‑1 (MQ=255)                    gTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGc                                                                                                                              >  4:15993/1‑149 (MQ=255)                     tttttAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGcc                                                                                                                             <  4:94441/149‑1 (MQ=255)                      ttttAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATAATCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGcc                                                                                                                             >  5:299709/1‑148 (MQ=255)                             aGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGt                                                                                                                     >  7:82023/1‑149 (MQ=255)                             aGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGt                                                                                                                     <  6:299709/149‑1 (MQ=255)                                aGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGtttt                                                                                                                  <  5:187863/149‑1 (MQ=255)                                  gATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTc                                                                                                                <  8:14100/149‑1 (MQ=255)                                       ttCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGAc                                                                                                           >  4:73322/1‑149 (MQ=255)                                                 gCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAgg                                                                                                  >  6:80584/1‑148 (MQ=255)                                                  cAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATa                                                                                                    >  5:12818/1‑145 (MQ=255)                                                  cAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATa                                                                                                    <  6:12818/145‑1 (MQ=255)                                                  cAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGaa                                                                                                >  5:222991/1‑149 (MQ=255)                                                  cAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGaa                                                                                                <  3:73322/149‑1 (MQ=255)                                                   aaaaTAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGaa                                                                                                >  5:304844/1‑148 (MQ=255)                                                     aataatGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATa                                                                                                    <  4:186813/142‑1 (MQ=255)                                                     aataatGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATa                                                                                                    >  3:186813/1‑142 (MQ=255)                                                      ataatGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGTCACCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTg                                                                                                       <  1:233529/138‑1 (MQ=255)                                                      ataatGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGTCACCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTg                                                                                                       <  1:229360/138‑1 (MQ=255)                                                      ataatGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCACCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTg                                                                                                       >  2:233529/1‑138 (MQ=255)                                                      ataatGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCACCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTg                                                                                                       >  2:229360/1‑138 (MQ=255)                                                          tGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATc                                                                                        <  2:182601/149‑1 (MQ=255)                                                             ttCGCCCTCTGACTGTTTCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCa                                                                                     <  4:144747/149‑1 (MQ=255)                                                               cGCCCTCTGACTGTTTCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCAc                                                                                    <  6:222991/148‑1 (MQ=255)                                                                gCCCTCTGACTGTTTCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCAccc                                                                                  >  2:109371/1‑149 (MQ=255)                                                                       gACTGTTTCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCACCCCACGCag                                                                           >  7:139021/1‑149 (MQ=255)                                                                        aCTGTTTCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTCAATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACATTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGATTCATGATCGAACCCCACg                                                                              >  1:82294/1‑146 (MQ=255)                                                                            tttCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCAGCCCACGCAGAGtca                                                                      >  3:203838/1‑149 (MQ=255)                                                                            tttCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCACCCCACGCAGAGtc                                                                       <  5:273354/148‑1 (MQ=255)                                                                               cGATGTCCAGTTGATTGGCCAGCTTGGTGATTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCACCCCACGCAGAGtcatc                                                                    <  2:82294/148‑1 (MQ=255)                                                                                gATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCACCCCACGCAGAGTCATCAc                                                                  <  8:82023/149‑1 (MQ=255)                                                                                gATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCACCCCACGCAGAGTCATCAc                                                                  >  3:266210/1‑149 (MQ=255)                                                                                  tGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCACCCCACGCAGAGTCATCACGt                                                                <  7:258147/149‑1 (MQ=255)                                                                                  tGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCACCCCACGCAGAGTCATCACGt                                                                <  3:15993/149‑1 (MQ=255)                                                                                  tGTCCAGGTGATTGGCCAGCTTGGTGAGTGCAATCCGGCCTTC‑ATACGCGTCTCCACGATTGCCGCCAGTGAGGAGCACAGTCGGGCCTATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCACCCCACGCAGAGTCATCACGt                                                                <  6:87644/149‑1 (MQ=255)                                                                                   gTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCGGTTTTTCATGACGGTGATAGGAATCATGATCGCACCCCACGCAGAGTCATCAcat                                                                <  6:304844/148‑3 (MQ=255)                                                                                   gTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCACCCCACGCAGAGTCATCACGt                                                                >  6:285932/1‑148 (MQ=255)                                                                                    tCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCACCCCACGCAGAGTCATCACGTGa                                                              <  8:139021/149‑1 (MQ=255)                                                                                      cAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAg                                                                                                   <  4:40098/110‑1 (MQ=255)                                                                                      cAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAg                                                                                                   >  3:40098/1‑110 (MQ=255)                                                                                        gTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCACCCCACGCAGAGTCATCACGTGAGt                                                            >  4:123391/1‑147 (MQ=255)                                                                                          tGATTGGCGAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCACCCCACGCAGAGTCATCACGTGAGTAAg                                                         <  1:109371/148‑1 (MQ=255)                                                                                                tccAGCTTGGTGAGTGCAATCGTGCCTTC‑ATACTCGTCTCCATGATTGCCGCTAGTGAGGATCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGAt                                                                                         <  1:245899/109‑1 (MQ=255)                                                                                                gacAGCTTGGTGAGTGCAATCGGGCCTTC‑AAACTCGTCTCCATGATTGCCGCCAGTGAGGAGCAAAGTCGGGCCGATACCTTTTTTCATGACGGTGATAGGAATCATGAt                                                                                         >  2:245899/3‑110 (MQ=255)                                                                                                   aGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCACCCCACGCAGAGTCATCACGTGAGTAAGGCAGCTTCaa                                               <  5:285932/149‑1 (MQ=255)                                                                                                               gCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCACCCCACGCAGAGTCATCACGt                                                                >  2:161753/1‑120 (MQ=255)                                                                                                               gCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCACCCCACGCAGAGTCATCACGt                                                                <  1:161753/120‑1 (MQ=255)                                                                                                                  aTCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCACCCCACGCAGAGTCATCACGTGAGTAAGGCAGCTTCAAGAACCCGTGCGacac                                >  6:176010/1‑149 (MQ=255)                                                                                                                     gggCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCACCCCACGCAGAGTCATCACGTGAGTAAGGCAGCTTCAAGAACCCGTGCGACAcccc                             <  5:80584/149‑1 (MQ=255)                                                                                                                              aTACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCACCCCACGCAGAGTCATCACGTGAGTAAGGCAGCTTCAAGAACCCGTGCGACACCCCTTGTTGGt                     >  4:262716/1‑149 (MQ=255)                                                                                                                               tACTCGTCTCCATGATTGCCGCCAGAGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCACCCCACGCAGAGTCATCACGTGAGTAAGGCAGCTTCAAGAACCCGTGCGACACCCCTTGTTGGTc                    <  2:137823/149‑1 (MQ=255)                                                                                                                                 cTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTc                                                                                                                >  6:313936/1‑55 (MQ=255)                                                                                                                                 cTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTc                                                                                                                <  5:313936/55‑1 (MQ=255)                                                                                                                                  tggTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCCATCCCGTTTTTCATGACGGTGATAGGGATCATGATCGCACCCCACGCAGCGTCAGCACCTGAGTAAGGCAGCTTCAAGCACCCGTGCGACACCCCTGTTTGGTCGa                  >  6:73386/3‑148 (MQ=255)                                                                                                                                    gTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCACCCCACGCAGAGTAATCACGTGAGTAAGGCAGCTTCAAGAACCCGTGCGACACCCCTTGTTGGTCGAAAt               >  6:202526/1‑149 (MQ=255)                                                                                                                                             aTTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCACCCCACGCAGAGTCATCACGTGAGTAAGGCAGCTTCAAGAACCCGTGCGACACCCCTTGTTGGTCGAAATCGATAGTTg      >  2:296638/1‑149 (MQ=255)                                                                                                                                                  cgccAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCACCCCACGCAGAGTCATCACGTGAGTAAGGCAGCTTCAAGAACCCGTGCGACACCCCTTGTTGGTCGAAATCGATAGTTGGCtt  >  3:82940/1‑148 (MQ=255)                                                                                                                                                  cgccAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCACCCCACGCAGAGTCATCACGTGAGTAAGGCAGCTTCAAGAACCCGTGCGACACCCCTTGTTGGTCGAAATCGATAGTTGGCtt  <  4:82940/148‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   TCGGAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTC‑ATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCACCCCACGCAGAGTCATCACGTGAGTAAGGCAGCTTCAAGAACCCGTGCGACACCCCTTGTTGGTCGAAATCGATAGTTGGCTT  >  NC_002947/5019802‑5020093

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