Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,182,072	A→C	7.0%	intergenic (‑130/+24)	PP_4560 ← / ← csbD	ribonuclease BN/stress response protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,182,072	0	A	C	7.0%	166.4 / 3.6	57	intergenic (‑130/+24)	PP_4560/csbD	ribonuclease BN/stress response protein
Reads supporting (aligned to +/- strand):  ref base A (29/24);  new base C (2/2);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

AGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAACACCCT  >  NC_002947/5181926‑5182214                                                                                                                                                    |                                                                                                                                                aGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCTC‑GGGTGAACCCgcga                                                                                                                                               <  6:131094/149‑1 (MQ=255)   gTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCTC‑GGGTGAACCCGcgaag                                                                                                                                             <  8:41437/149‑1 (MQ=255)    tCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGC‑GGGTTCACCCGagaaga                                                                                                                                            <  5:348421/149‑1 (MQ=255)       gggAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGC‑GGGTTCACCCGAGAAGAGGc                                                                                                                                         >  8:278319/1‑149 (MQ=255)                          cTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCtcaag                                                                                                                          >  8:9541/1‑145 (MQ=255)                          cTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCtcaag                                                                                                                          <  7:9541/145‑1 (MQ=255)                                ccgcTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCCAGGcttc                                                                                                                >  2:52218/1‑149 (MQ=255)                                  gcTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCTC‑GGGTGAACCCGcgaa                                                                                                                                              <  8:131937/117‑1 (MQ=39)                                  gcTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCtcaagtcaag                                                                                                                     <  5:84660/142‑1 (MQ=255)                                  gcTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCtcaagtcaag                                                                                                                     >  6:84660/1‑142 (MQ=255)                                        tCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCg                                                                                                         >  6:120298/1‑148 (MQ=255)                                          gTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGGGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCa                                                                                                       <  5:257254/148‑1 (MQ=255)                                              cacaGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCt                                                                                                   >  4:214911/1‑148 (MQ=255)                                                 aGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCTC‑GGGTGAACCCGCGAAGAGGCCGGc                                                                                                                                     >  6:271093/1‑111 (MQ=38)                                                 aGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCTC‑GGGTGAACCCGCGAAGAGGCCGGc                                                                                                                                     <  5:271093/111‑1 (MQ=38)                                                 aGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGC‑GGGTTCACCCGagaa                                                                                                                                              >  3:73449/1‑102 (MQ=255)                                                 aGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGC‑GGGTTCACCCGagaa                                                                                                                                              <  4:73449/102‑1 (MQ=255)                                                 aGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTcac                                                                                               <  7:278319/149‑1 (MQ=255)                                                          tAATCGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCttcttc                                                                                                             <  3:51847/126‑1 (MQ=255)                                                            aTAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTgc                                                                                   >  8:262669/1‑148 (MQ=255)                                                                 ccGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTGACctgc                                                                              <  3:214911/149‑1 (MQ=255)                                                                     gCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCt                                                                                      >  2:301039/1‑137 (MQ=255)                                                                     gCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCt                                                                                      <  1:301039/137‑1 (MQ=255)                                                                      cccGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTGACCTGCTgcgc                                                                         >  6:109546/1‑149 (MQ=255)                                                                        cGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTGACCTGCTGcgcct                                                                       >  1:149879/1‑149 (MQ=255)                                                                                      cccGGCCGTGGTTGAAGAATTTGCTTACGAATAGGTGCCGGCCTCTTCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCtcaagtcaa                                                                                                                      >  6:331941/1‑89 (MQ=255)                                                                                      cccGGCCGTGGTTGAAGAATTTGCTTACGAATACGTGCCGGCCTCTTCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCtcaagtcaa                                                                                                                      <  5:331941/89‑1 (MQ=255)                                                                                      cccGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTGACCTGCTGCGCCTCGCCTTTGAGCTcc                                                         <  5:120298/149‑1 (MQ=255)                                                                                       ccGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCt                                                        >  5:180791/1‑149 (MQ=255)                                                                                       ccGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCt                                                        <  6:180791/149‑1 (MQ=255)                                                                                       ccGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGC‑GGGTTCACCCGAAAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCt                                                        >  1:101214/1‑149 (MQ=255)                                                                                         ggCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTGACCTGCTgcgc                                                                         >  7:119763/1‑130 (MQ=255)                                                                                          gCCGTGGTTGAATAATTTGCTTACGAATAACTGCCGGCCTGTTCGC‑GGGCTCACCCGATAAGCGGCCGGCACCGCCTGAAGTCAAGGCTTCTTCACCCCATCCATCACATTTCTCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTgc                                                      <  1:52218/149‑1 (MQ=255)                                                                                             gTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCgcct                                                  <  3:270049/149‑1 (MQ=255)                                                                                             gTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTGACCTGCTGCGCCTCACCTTTGAGCTCCTGCgcct                                                  >  1:157180/1‑149 (MQ=255)                                                                                                gTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCgccttg                                                >  1:243201/1‑149 (MQ=255)                                                                                                   gaagaaTTTGCTTACGAATAAGTGCCGGCCTCTTCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCgccttgcct                                             <  2:243201/148‑1 (MQ=255)                                                                                                        aTTTGCTTACGAATAAGTGCCGGCCTCTTCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGc                                        <  1:8506/148‑1 (MQ=255)                                                                                                            gCTTACCAATAAGTGCCGGCCTCTTCGC‑GGGTTCACCCGAGAAGTGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGc                                    <  2:101214/148‑1 (MQ=255)                                                                                                              ttACGAATAAGTGCCGGCCTCTTCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGc                                 >  7:181562/1‑149 (MQ=255)                                                                                                              ttACGAATAAGTGCCGGCCTCTTCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGc                                 >  3:176883/1‑149 (MQ=255)                                                                                                                aCGAATAAGTGCCGGCCTCTTCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCAACGCATCCTTCACATTTC‑CCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCtt                               >  2:22902/1‑149 (MQ=255)                                                                                                                 cGAATAAGTGCCGGCCTCTTCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCttt                              >  2:53713/1‑149 (MQ=255)                                                                                                                        gTGCCGGCCTCTTCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCgccttgc                                               <  7:273202/125‑1 (MQ=255)                                                                                                                        gTGCCGGCCTCTTCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCgccttgc                                               >  8:273202/1‑125 (MQ=255)                                                                                                                          gCCGGCCTCTTCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCg                     >  4:182536/1‑149 (MQ=255)                                                                                                                          gCCGGCCTCTTCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCg                     <  4:176883/149‑1 (MQ=255)                                                                                                                           ccGGCCTCTTCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCgccttgc                                               <  8:234186/122‑1 (MQ=255)                                                                                                                             ggCCTCTTCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGt                   >  7:20794/1‑148 (MQ=255)                                                                                                                              gCCTCTTCTC‑GGGTGAACCCGCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTGa                                                                                   <  2:64185/83‑1 (MQ=25)                                                                                                                              gCCTCTTCTC‑GGGTGAACCCGCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTGa                                                                                   >  1:64185/1‑83 (MQ=25)                                                                                                                                   ttCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTc                            >  2:144823/1‑133 (MQ=255)                                                                                                                                   ttCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTc                            <  1:144823/133‑1 (MQ=255)                                                                                                                                   ttCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCg                           <  6:23590/134‑1 (MQ=255)                                                                                                                                   ttCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCg                           >  5:23590/1‑134 (MQ=255)                                                                                                                                   ttCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTg            >  4:146166/1‑149 (MQ=255)                                                                                                                                      gc‑gGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGcc          >  4:176920/1‑148 (MQ=255)                                                                                                                                          ggTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAaca      <  5:109546/149‑1 (MQ=255)                                                                                                                                            ttCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAACAcc    <  2:149879/149‑1 (MQ=255)                                                                                                                                            ttCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAACAcc    >  5:78653/1‑149 (MQ=255)                                                                                                                                               aCCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAACACCCt  <  1:53713/148‑1 (MQ=255)                                                                                                                                                 ccGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCg                     <  6:280396/127‑1 (MQ=255)                                                                                                                                                 ccGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCg                     >  5:280396/1‑127 (MQ=255)                                                                                                                                                    |                                                                                                                                                AGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGC‑GGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTC‑CCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAACACCCT  >  NC_002947/5181926‑5182214

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