Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,143,883	T→A	7.1%	intergenic (+36/‑87)	PP_3647 → / → PP_3648	oxidoreductase/carboxymuconolactone decarboxylase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,143,883	0	T	A	7.1%	162.0 / 2.7	61	intergenic (+36/‑87)	PP_3647/PP_3648	oxidoreductase/carboxymuconolactone decarboxylase family protein
Reads supporting (aligned to +/- strand):  ref base T (29/26);  new base A (2/2);  total (31/30)
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

CACGCGATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAG‑CAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGAG  >  NC_002947/4143740‑4144026                                                                                                                                                  |                                                                                                                                                cACGCGATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAG‑CAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCcg                                                                                                                                            >  3:155565/1‑149 (MQ=255) cACGCGATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAG‑CAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCAAGCCcg                                                                                                                                            >  1:234791/1‑149 (MQ=255)    gcgATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAG‑CAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCcgcg                                                                                                                                          >  5:369902/1‑148 (MQ=255)     cgATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAG‑CAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCCGCGaa                                                                                                                                        <  3:98755/149‑1 (MQ=255)     cgATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAG‑CAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCCGCGaa                                                                                                                                        <  5:214408/149‑1 (MQ=255)      gATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAG‑CAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCcgcg                                                                                                                                          >  7:28029/1‑146 (MQ=255)          aTTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAG‑CAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCAAGCCCGCGAAGATg                                                                                                                                    <  2:234791/148‑1 (MQ=255)           ttGTCGGCGCGGTTAAGAACCTGGGGTTGCCGGAG‑CAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGcc                                                                                                                                  <  6:5132/149‑1 (MQ=255)           ttGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGTCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGc                                                                                                                                   >  1:244955/1‑149 (MQ=255)                ggCGCGGTGAAGAACCTGGTGTTGCCTGAG‑CAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCAAGCCCGCGAAGATGccaacat                                                                                                                             <  4:94640/149‑6 (MQ=255)                     ggTGAAGAATCTGGTGTTGCCTGAG‑CAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCa                                                                                                                        <  3:168411/149‑1 (MQ=255)                     ggTGAAGAACCTGGTGTTGCCTGAG‑CAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCa                                                                                                                        <  7:100268/149‑1 (MQ=255)                             aCCTGGTGTTGCCTGAG‑CAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGTGTcgc                                                                                                                >  1:272416/1‑149 (MQ=255)                                     ttGCCTGAG‑CAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGCCAGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATc                                                                                                         <  1:298753/149‑1 (MQ=255)                                         cTGAG‑CAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGtct                                                                                                     >  6:89129/1‑148 (MQ=255)                                         cTGAG‑CAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGtct                                                                                                     >  2:192718/1‑148 (MQ=255)                                           gAG‑CAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCAc                                                                                                  >  2:240444/1‑149 (MQ=255)                                                          cccGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGt                                                                                                       <  8:335077/130‑1 (MQ=255)                                                          cccGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGt                                                                                                       >  7:335077/1‑130 (MQ=255)                                                          cccGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGcc                                                                                    >  4:320767/1‑149 (MQ=255)                                                                ggTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCtttt                                                                              >  5:274634/1‑149 (MQ=255)                                                                ggTGTACCACAAGGGCAAGCTCGGCGCCCTGCAGCCAATTGGCTGATACATCTTGAGTC‑GATGTTGGCATCTTAGCGGGCTTGCCCGCGAAGATGCACCTGTTTTCAGGGGTAGAGCTAATCGTCTCACCCTCGTTTGCGGGCCGtttt                                                                              >  7:152569/1‑149 (MQ=255)                                                                 gtgtACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCttttt                                                                             <  2:244955/149‑1 (MQ=255)                                                                 gtgtACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCttttt                                                                             >  8:346540/1‑149 (MQ=255)                                                                 gtgtACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCGGATACATCTTGAGGC‑GATGTGGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCttttt                                                                             <  7:38985/149‑1 (MQ=255)                                                                  tgtACCACAAGTGCAAGCTCTGAGCCCTGCCGCCACTTGGCTGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCGTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCttttt                                                                             <  3:44170/148‑1 (MQ=255)                                                                      ccACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCttttt                                                                             >  6:281213/1‑144 (MQ=255)                                                                      ccACAAGGGCAAGCGCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCttttt                                                                             <  5:281213/144‑1 (MQ=255)                                                                       cacaAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCATCGGTTGCGGCCCCTTTTTTATGcc                                                                       >  8:320352/1‑149 (MQ=255)                                                                       cacaAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTCAGTC‑GATGTTGGCATCTTCCCCGGCTTGCCCCCGAAGATGCCCCTGTTTTCAGGGGTCGCGCAAATAGACTCACCCACCGTTGCGGCCCCAGTTTTaatcc                                                                       >  6:52855/1‑145 (MQ=255)                                                                        acaAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCg                                                                      <  6:274634/149‑1 (MQ=255)                                                                        acaAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCg                                                                      <  1:240444/149‑1 (MQ=255)                                                                          aaGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGa                                                                    <  5:335291/149‑1 (MQ=255)                                                                          aaGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGa                                                                    <  4:155565/149‑1 (MQ=255)                                                                          aaGGGCAAGCGCGGCGCCCTGCCGCCACTTGGCTGAAACATCTTCAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGa                                                                    <  4:255177/149‑1 (MQ=255)                                                                           aGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAg                                                                   <  6:369902/149‑1 (MQ=255)                                                                             ggCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTa                                                                 >  3:184147/1‑149 (MQ=255)                                                                                aaGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACt                                                               >  2:350186/1‑148 (MQ=255)                                                                                         gcCCTGCCGCCACTTGGCTGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAg                                                      >  5:161187/1‑148 (MQ=255)                                                                                         gcCCCGCCGCCACTTGGCAGCGACCTCTTGAGAC‑CCGGTTGGCATCTTCGCGGGCGTGCCCGCGAAGATGCCCCGGTTTTCAGGGGTCGCGCCAAGCGTCTCCCCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTCCTGCACATGAGc                                                     <  4:243863/149‑1 (MQ=255)                                                                                           ccTGCCGCCACTTGGCTGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCaa                                                   <  5:89129/149‑1 (MQ=255)                                                                                                 ggaaCTTGGCTTATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCCGCTATGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCGCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGaaaa                                             <  1:155819/146‑1 (MQ=255)                                                                                                    aCTTGGCTGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTTCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTAACCCTCGGTTGCGGACCGGTTTTTATTCCGGAGTACTGCAGATGAGGAATGAAAAGt                                           >  8:296623/1‑148 (MQ=255)                                                                                                           tGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAg                                   >  7:354078/1‑149 (MQ=255)                                                                                                           tGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAg                                   >  1:177783/1‑149 (MQ=255)                                                                                                           tGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAg                                   >  3:196899/1‑149 (MQ=255)                                                                                                            gATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAgg                                  >  2:295173/1‑149 (MQ=255)                                                                                                                   cTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGaa                                      >  4:22854/1‑138 (MQ=255)                                                                                                                   cTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGaa                                      <  3:22854/138‑1 (MQ=255)                                                                                                                        gTC‑GATGTTGTCATCTTCGAGGTCTTGCCCGAGAAGATGCCCCTTTTTTAATGTGTCGCGCTAATCGTCTAACACTCGTTTGCGGCCCATTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATc                       <  7:320352/148‑1 (MQ=255)                                                                                                                        gTC‑GATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATcc                      >  4:180035/1‑149 (MQ=255)                                                                                                                             aTGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGGTTTCAGGGGTCGCGCTAATCGTCTCACCCGCGGTTGCGGCCCCTGTTTTATGCCGGAGTACTGGACATGAGGAATGAAAAGGACGATACGGGCCGGCAAATGCGCAc                  >  8:136392/1‑149 (MQ=255)                                                                                                                                   gCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAgg             <  1:350186/148‑1 (MQ=255)                                                                                                                                       agtCGCGGGCACGCCCGCGAAGATGACCCTGTTTTCATGGGTCGCGCTAATCCTCTCACCCTCGGTTGCGGCCACTTTTTTATTCCTGAGTACTTCACATGAGCAATGAAAAGTACGAAAATGGCCTGAAAATCCGCACCCATGTACTg        <  8:152569/147‑1 (MQ=255)                                                                                                                                        ttCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTgg       >  8:168168/1‑149 (MQ=255)                                                                                                                                        ttCGCGGGCAAGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTgg       >  8:293169/1‑149 (MQ=255)                                                                                                                                          cgcgGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTggg      <  3:320767/148‑1 (MQ=255)                                                                                                                                            cgGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGTCTGAAAATCAGCACCCAGGTACTGGGc     >  8:69630/1‑147 (MQ=255)                                                                                                                                            cgGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGa   <  4:196899/149‑1 (MQ=255)                                                                                                                                             gggCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGAg  <  8:28029/149‑1 (MQ=255)                                                                                                                                             gggCTGCCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGAg  <  7:296623/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                CACGCGATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAG‑CAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTC‑GATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGAG  >  NC_002947/4143740‑4144026

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑