Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,041,195	C→G	6.0%	intergenic (+173/‑489)	PP_5680 → / → PP_4447	Cro/CI family transcriptional regulator/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,041,195	0	C	G	6.0%	215.6 / 2.6	66	intergenic (+173/‑489)	PP_5680/PP_4447	Cro/CI family transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (34/28);  new base G (2/2);  total (36/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 = 8.06e-01

TTGTACTCCACACCCAGTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGA  >  NC_002947/5041052‑5041340                                                                                                                                                |                                                                                                                                                   ttGTACTCCACACCCAGTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCt                                                                                                                                                >  6:331435/1‑148 (MQ=255)   gTACTCCACACCCAGTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCtt                                                                                                                                               >  8:155964/1‑147 (MQ=255)      cTCCACACCCAGTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGcc                                                                                                                                          >  3:43386/1‑149 (MQ=255)         cacaCCCAGTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCccgc                                                                                                                                       >  3:2715/1‑149 (MQ=255)         cacaCCCAGTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTC‑TCGGCGGGCTTAAGGCgcgcc                                                                                                                                      >  1:69798/1‑149 (MQ=255)         cacaCCCAGTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCAACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCccgc                                                                                                                                       >  2:2033/1‑149 (MQ=255)              ccAGTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCga‑g                                                                                                                                  >  8:85744/1‑148 (MQ=255)               cAGTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGa‑gaa                                                                                                                                >  2:60922/1‑149 (MQ=255)                aGTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTC‑TCGGCGGGCTTAAGGCGCGCCgatgaag                                                                                                                               <  2:69798/149‑6 (MQ=255)                 gTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGa‑gaaga                                                                                                                              >  5:284575/1‑149 (MQ=255)                    tcACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCAACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGc                                                                                                                            >  6:161050/1‑148 (MQ=255)                     cacaTCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGTGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGc                                                                                                                            >  6:87869/1‑147 (MQ=255)                      acaTCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGTGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGc                                                                                                                            <  1:60922/146‑1 (MQ=255)                       caTCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTC‑TCGGCGGGCTTAAGGCGCGCCGATGAagaggtcat                                                                                                                        <  5:63289/149‑6 (MQ=37)                                   tCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGc                                                                                                              >  2:226282/1‑147 (MQ=255)                                      gCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGAc                                                                                                         <  7:85744/149‑1 (MQ=255)                                       cGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGAcc                                                                                                        <  7:267342/149‑1 (MQ=255)                                           gCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCt                                                                                                     >  2:224520/1‑148 (MQ=255)                                                     gtgtTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTAt                                                                                          <  5:331435/149‑1 (MQ=255)                                                       gtTGAGCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCTCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATgg                                                                                        >  7:311345/1‑149 (MQ=255)                                                           aTCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGACCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGCCCCCTAACCGCCTTATGGCGGt                                                                                    >  7:5695/1‑149 (MQ=255)                                                              cacaGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGcc                                                                                              <  2:14624/136‑1 (MQ=255)                                                              cacaGAGGGTTAGGGATAGCCTTAAGGCATAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGcc                                                                                              >  1:14624/1‑136 (MQ=255)                                                               acaGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTtct                                                                                >  2:56375/1‑149 (MQ=255)                                                                caGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTtctc                                                                               <  8:311345/149‑1 (MQ=255)                                                                 agagGGTTAGTGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTctct                                                                              <  5:87869/149‑1 (MQ=255)                                                                 agagGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGTACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTctct                                                                              <  4:43386/149‑1 (MQ=255)                                                                   agGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGc                                                                            <  6:284575/149‑1 (MQ=255)                                                                       ttAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTcc                                                                        >  1:116760/1‑149 (MQ=255)                                                                          gggATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCa                                                                      <  7:155964/148‑1 (MQ=255)                                                                          gggATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCa                                                                      <  1:56375/148‑1 (MQ=255)                                                                           ggATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAg                                                                     >  4:221575/1‑148 (MQ=255)                                                                           ggATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGt                                                                    >  3:50634/1‑149 (MQ=255)                                                                             aTAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTct                                                                  <  3:173872/149‑1 (MQ=255)                                                                                     aaGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGa                                                           <  5:161050/148‑1 (MQ=255)                                                                                      aGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAg                                                         <  3:141408/149‑1 (MQ=255)                                                                                      aGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAg                                                         <  1:2033/149‑1 (MQ=255)                                                                                          gTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGc                                                       >  7:103911/1‑147 (MQ=255)                                                                                          gTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGc                                                       <  8:103911/147‑1 (MQ=255)                                                                                             aTTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTTCGAAGGCCGATa                                                  <  8:5695/149‑1 (MQ=255)                                                                                                 tCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGcc                                               <  4:2715/148‑1 (MQ=255)                                                                                                   tGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTctct                                                                              >  7:133696/1‑115 (MQ=255)                                                                                                   tGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTctct                                                                              <  8:133696/115‑1 (MQ=255)                                                                                                         cAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTccaac                                      >  1:213175/1‑149 (MQ=255)                                                                                                           acacCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTctct                                                                              <  4:268839/107‑1 (MQ=255)                                                                                                           acacCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTctct                                                                              >  3:268839/1‑107 (MQ=255)                                                                                                             acCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCCCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAAt                                  >  5:238683/1‑149 (MQ=255)                                                                                                              ccTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAAtt                                 <  4:50634/149‑1 (MQ=255)                                                                                                                tctGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCCA‑GAAGAGCCAGAGACCCCGAGCAAGCCCCCTAACCGCCTTATCGCGGTTTATCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCACCCAATTc                                >  3:95081/1‑148 (MQ=255)                                                                                                                         aGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGa                      <  1:224520/149‑1 (MQ=255)                                                                                                                          ggACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCCCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGAc                     <  6:238683/149‑1 (MQ=255)                                                                                                                                gACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGcc                                                      >  4:32490/1‑110 (MQ=255)                                                                                                                                gACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGcc                                                      <  3:32490/110‑1 (MQ=255)                                                                                                                                gACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATa                                                  >  7:61799/1‑114 (MQ=255)                                                                                                                                gACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATa                                                  <  8:61799/114‑1 (MQ=255)                                                                                                                                gACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGcc               >  5:256183/1‑149 (MQ=255)                                                                                                                                gACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGcc               >  7:317740/1‑149 (MQ=255)                                                                                                                                    tctTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGAc           >  7:170402/1‑149 (MQ=255)                                                                                                                                     ctTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCg                                                     <  2:72803/106‑1 (MQ=255)                                                                                                                                     ctTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCg                                                     >  1:72803/1‑106 (MQ=255)                                                                                                                                      ttCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCg         >  6:185134/1‑149 (MQ=255)                                                                                                                                         ctCGGCGGGCTTAAGGCGCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAg      >  1:11385/2‑149 (MQ=255)                                                                                                                                          tCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGa            <  8:88470/142‑1 (MQ=255)                                                                                                                                          tCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGa            >  7:88470/1‑142 (MQ=255)                                                                                                                                            ggCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGa                                             <  1:79350/107‑1 (MQ=255)                                                                                                                                            ggCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGa                                             >  2:79350/1‑107 (MQ=255)                                                                                                                                             gcgcgcCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTaga  <  8:317740/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                   TTGTACTCCACACCCAGTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGA  >  NC_002947/5041052‑5041340

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