Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,075,382	A→G	10.0%	intergenic (+146/+174)	PP_4468 → / ← PP_t58	Cro/CI family transcriptional regulator/tRNA‑Arg

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,075,382	0	A	G	10.0%	155.9 / 8.1	60	intergenic (+146/+174)	PP_4468/PP_t58	Cro/CI family transcriptional regulator/tRNA‑Arg
Reads supporting (aligned to +/- strand):  ref base A (25/29);  new base G (2/4);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.87e-01

TCAGGCACAAAGGCGGACAAGCAAACGTAGATCAATTCCACTCTGCCCATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCG‑AC‑ACCCATAGGAGAAGCGCCGAGAGCAGGAACTTCCACCCCACTCAACTCAAGCAAAGCGGATATGAAGAACACAGAGAAGGCAGCAAACCTGGTGCTTC  >  NC_002947/5075237‑5075526                                                                                                                                                  |                                                                                                                                                   tCAGGCACAAAGGCGGACAAGCAAACGTAGATCAATTCCACTCTGCCCATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGcc                                                                                                                                                 >  8:286029/1‑149 (MQ=255)     gCACAAAGGCGGACAAGCAAACGTAGATCAATTCCACTCTGCCCATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAg                                                                                                                                             >  1:105001/1‑149 (MQ=255)      cacaAAGGCGGACAAGCAAACGTAGATCAATTCCACTCTGCCCATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCTGGTCGGCCGaact                                                                                                                                            <  7:94233/149‑3 (MQ=255)       acaAAGGCGGACAAGCAAACGTAGATCAATTCCACTCTGCCCATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCTGATCGGCCGaactc                                                                                                                                           <  6:282117/149‑4 (MQ=255)       acaAAGGCGGACAAGCAAACGTAGATCAATTCCACTCTGCCCATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCTGATCGGCCGaactc                                                                                                                                           <  7:199308/149‑4 (MQ=255)       acaAAGGCGGACAAGCAAACGTAGATCAATTCCACTCTGCCCATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCa                                                                                                                                           >  3:17310/1‑149 (MQ=255)       acaAAGGCGGACAAGCAAACGTAGAACAAATCCACTCTGCCCATCAATCGAACGACATCCATTACCCCCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGATGCTGGAAAGCGAGTACGGCCGATCACCCGAAGCa                                                                                                                                           >  3:15556/1‑149 (MQ=255)        caAAGGCGGACAAGCAAACGTAGATCAATTCCACTCTGCCCATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCAt                                                                                                                                          <  7:38116/149‑1 (MQ=255)         aaaGGCGGACAAGCAAACGTAGATCAATTCCACTCTGCCCATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCTGATCGGCCGaactcgc                                                                                                                                         <  1:236692/149‑6 (MQ=255)         aaaGGCGGACAAGCAAACGTAGATCAATTCCACTCTGCCCATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATg                                                                                                                                         >  6:351657/1‑149 (MQ=255)                gACAAGCAAACGTAGATCAATTCCACTCTGCCCATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAagag                                                                                                                                  <  7:71234/149‑1 (MQ=255)                gACAAGCAAACGTAGATCAATTCCACTCTGCCCATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAagag                                                                                                                                  <  1:321694/149‑1 (MQ=255)                     gCAAACGTAGATCAATTCCACTCTGCCCATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCg                                                                                                                               >  5:200954/1‑147 (MQ=255)                      cAAACGTAGATCAATTCCACTCTGCCCATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGATATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCg                                                                                                                               >  6:286215/1‑146 (MQ=255)                       aaaCGTAGATCAATTCCACTCTGCCCATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTaga                                                                                                                           <  6:58121/149‑1 (MQ=255)                           gTAGATCAATTCCACTCTTCCCATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTagaagag                                                                                                                       <  8:179238/149‑1 (MQ=255)                           gTAGATCAATTCCACTCTGCCCATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTagaagag                                                                                                                       <  8:210921/149‑1 (MQ=255)                              gATCAATTCCACTATGCCCATCAATCGAACGACATCCATTTCCCGCAGCATCTCATATATTTTCCAGGAATTTGCCAGGGTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTagaagaga                                                                                                                      <  4:15556/147‑1 (MQ=255)                                tCAATTCCACTCTGCCCATCAATCGAACGCCATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGCTCAGCCGAAGCATGTTACGAGGCGTAGAAGAGAAGct                                                                                                                  >  6:352245/1‑148 (MQ=255)                                   aTTCCACTCTGCCCATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGa                                                                                                                >  6:130611/1‑148 (MQ=255)                                    ttCCACTCTGCCCATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGaa                                                                                                               >  7:143648/1‑148 (MQ=255)                                          tctGCCCATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGATAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATAAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGccc                                                                                                         >  1:43914/1‑148 (MQ=255)                                              cccATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGc                                                                                                                                <  2:335983/121‑1 (MQ=255)                                              cccATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGc                                                                                                                                >  1:335983/1‑121 (MQ=255)                                              cccATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCG‑ac‑ac                                                                                                  <  5:130611/149‑1 (MQ=255)                                              cccATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCG‑ac‑ac                                                                                                  <  5:352245/149‑1 (MQ=255)                                               ccATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGccc                                                                                                         <  2:13424/143‑1 (MQ=255)                                               ccATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGccc                                                                                                         >  1:13424/1‑143 (MQ=255)                                                  tCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCG‑AC‑ACCCAt                                                                                              <  5:351657/149‑1 (MQ=255)                                                    aaTCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCG‑AC‑ACCCATAg                                                                                            <  7:286029/149‑1 (MQ=255)                                                       cGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCG‑AC‑ACCCATAg                                                                                            >  5:130131/1‑146 (MQ=255)                                                         aaCGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCG‑AC‑ACCCATAGGAGaa                                                                                       <  6:200954/149‑1 (MQ=255)                                                         aaCGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCG‑AC‑ACCCATAGGAGaa                                                                                       <  8:38106/149‑1 (MQ=255)                                                         aaCGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGCAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGACGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGATAAGCAGAATAGCCCG‑AC‑ACACATAGgag                                                                                         >  5:351421/1‑147 (MQ=255)                                                             aCATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGcc                                                                                                          >  6:240669/1‑128 (MQ=255)                                                             aCATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGcc                                                                                                          <  5:240669/128‑1 (MQ=255)                                                              cATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCG‑AC‑ACCCATAGGAGAAGCGcc                                                                                  <  4:17310/149‑1 (MQ=255)                                                                tCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCG‑AC‑ACCCATAGGAGAAGCGcc                                                                                  >  6:161483/1‑147 (MQ=255)                                                                tCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCG‑AC‑ACCCATAGGAGAAGCGcc                                                                                  >  2:107888/1‑147 (MQ=255)                                                                  cATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCG‑AC‑ACCCATAGGAGAAGCGCCgaga                                                                              <  2:43914/149‑1 (MQ=255)                                                                                  tcAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCG‑AC‑ACCCATAGGAGAAGCGCCGAGAGCAGGAACTTCCAccc                                                              >  2:182913/1‑149 (MQ=255)                                                                                  tcAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCG‑AC‑ACCCATAGGAGAAGCGCCGAGAGCAGGAACTTCCAcc                                                               >  5:289916/1‑148 (MQ=255)                                                                                    aGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCG‑AC‑ACCCATAGGAGAAGCGCCGAGAGCAGGAACTTCCACCCCa                                                            <  6:130131/149‑1 (MQ=255)                                                                                       aaTTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCG‑AC‑ACCCATAGGAGAAGCGCCGAGATCAGGAACTTCCACCCCactc                                                         <  6:289916/149‑1 (MQ=255)                                                                                            tCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCG‑AC‑ACCCATAGGAGAAgcgc                                                                                   >  1:91694/1‑118 (MQ=255)                                                                                            tCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCG‑AC‑ACCCATAGGAGAAgcgc                                                                                   <  2:91694/118‑1 (MQ=255)                                                                                              cAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTACGCGAGTGCGGCCGATCAGCCGAAGCATGTTAAGAGGGGTAGAAGAGAAGCAGAATAGCCCG‑AC‑ACCCATAGGAGAAGCGCCGAGAGCAGGAACTTCCACCCCActcaactcaa                                                  <  8:143648/149‑1 (MQ=255)                                                                                               aGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCG‑AC‑ACCCATAGGAGAAGCGCCGAGAGCAGGAACTTCCACCCCACTCAACTCAAg                                                 >  8:122177/1‑149 (MQ=255)                                                                                                  aaaTTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCG‑AC‑ACCCATAGGAGAAGCGCCGAGAGCa                                                                           <  6:120356/120‑1 (MQ=255)                                                                                                  aaaTTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCG‑AC‑ACCCATAGGAGAAGCGCCGAGAGCa                                                                           >  5:120356/1‑120 (MQ=255)                                                                                                    aTTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCG‑AC‑ACCCATAGGAGaa                                                                                       <  5:333885/106‑1 (MQ=255)                                                                                                    aTTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCG‑AC‑ACCCATAGGAGaa                                                                                       >  6:333885/1‑106 (MQ=255)                                                                                                                ggAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCG‑AC‑ACCCATAGGAGAAGCGCCGAGAGCAGGAACTTCCACCCCACTCAACTCAAGCAAAGCGGATATgaag                                 >  5:230127/1‑148 (MQ=255)                                                                                                                    ttatctgcttGTAATCGAGTTCCGCCGATCATCCGAATGATTTTAATATGCGTATAATTTAATCAAAATATCCCGTACAATCCATAGGAGAAGCGCCGAGAGCAGGAACTTCCACCCCACTCAACTCAAGCAAAGCGGATATgaagaa                               <  1:107888/144‑1 (MQ=255)                                                                                                                     gAGCTCCTGGAAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCG‑AC‑ACCCATAGGAGAAGCGCCGAGAGCAGGAACTTCCACCCCACTCAACTCAAGCAAAGCGGATATGAAGAAcaca                           <  5:257267/149‑1 (MQ=255)                                                                                                                      aGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCG‑AC‑ACCCATAGGAGAAGCGCCGAGAGCAGGAACTTCCACCCCACTCAACTCAAGCAAAGCGGATATGAAGAAcaca                           <  1:182913/148‑1 (MQ=255)                                                                                                                                  catgcttCGGCTGATCGGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCG‑AC‑ACCCATAGGAGAAGCGCCGAGAGCAGGAACTTCCACCCCACTCAACTCAAGCAAAGCGGATATGAAGAACACAGAGAAGgcagca               >  7:185628/6‑148 (MQ=255)                                                                                                                                      cttCGGCTGATCGGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCG‑AC‑ACCCATAGGAGAAGCGCCGAGAGCAGGAACTTCCACCCCACTCAACTCAAGCAAAGCGGATATGAAGAACACAGAGAAGGCAGCAAACCt          >  1:44850/2‑149 (MQ=255)                                                                                                                                           gCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCG‑AC‑ACCCATAGGAGAAGCGCCGAGAGCAGGAACTTCCACCCCACTCAACTCAAGCAAAGCGGATATGAAGAACACAGAGAAGGCAGCAAACCTGGTGc     <  2:105001/149‑1 (MQ=255)                                                                                                                                              gATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCG‑AC‑ACCCATAGGAGAAGCGCCGAGAGCAGGAACTTCCACCCCACTCAACTCAAGCAAAGCGGATATGAAGAACACAGAGAAGGCAGCAAACCTGGTGCTTc  <  5:286215/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                   TCAGGCACAAAGGCGGACAAGCAAACGTAGATCAATTCCACTCTGCCCATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCG‑AC‑ACCCATAGGAGAAGCGCCGAGAGCAGGAACTTCCACCCCACTCAACTCAAGCAAAGCGGATATGAAGAACACAGAGAAGGCAGCAAACCTGGTGCTTC  >  NC_002947/5075237‑5075526

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