Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,041,200	T→A	9.0%	intergenic (+178/‑484)	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,200	0	T	A	9.0%	219.8 / 10.7	78	intergenic (+178/‑484)	PP_5680/PP_4447	Cro/CI family transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (40/31);  new base A (5/2);  total (45/33)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.64e-01

GTACTCCACACCCAGTCTCACATCTGACCCT‑GTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGT  >  NC_002947/5041054‑5041346                                                                                                                                                    |                                                                                                                                                    gTACTCCACACCCAGTCTCACATCTGACGCT‑GTCAGCGACTCTAATGACGGTTTTGATCCAAAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCtt                                                                                                                                                     <  3:159726/147‑1 (MQ=255)    cTCCACACCCAGTCTCACATCTGACCCT‑GTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGcc                                                                                                                                                >  5:51526/1‑149 (MQ=255)      ccACACCCAGTCTCACATCTGACCCTGGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGccc                                                                                                                                               >  1:17665/1‑149 (MQ=255)      ccACACCCAGTCTCACATCTGACCCT‑GTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCccg                                                                                                                                              >  3:191619/1‑149 (MQ=255)        acacCCAGTCTCACATCTGACCCT‑GTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCccgcc                                                                                                                                            >  6:8605/1‑149 (MQ=255)          acCCAGTCTCACATCTGACCCT‑GTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCccgcc                                                                                                                                            >  2:231315/1‑147 (MQ=255)           cccAGTCTCACATCTGACCCT‑GTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTC‑TCGGCGGGCTTAAGGCGCGCCgatg                                                                                                                                        <  2:332621/149‑3 (MQ=255)            ccAGTCTCACATCTGACCCT‑GTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGGGCTTAAGGCGCGCCgat                                                                                                                                         >  7:171640/1‑147 (MQ=255)            ccAGTCTCACATCTGACCCT‑GTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGGGCTTAAGGCGCGCCgat                                                                                                                                         <  8:171640/148‑2 (MQ=255)            ccAGTCTCACATCTGACCCT‑GTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCga‑g                                                                                                                                        >  7:68216/1‑148 (MQ=255)            ccAGTCTCACATCTGACCCT‑GTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGAGGCCAGCACGTCTGCGATGAGGACATGACCTCTTCAGCGGCGCGCCGTAGGGCCGCCga‑g                                                                                                                                        >  6:338546/1‑148 (MQ=255)               gTCTCACATCTGACCCT‑GTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGa‑gaaga                                                                                                                                    <  5:8605/149‑1 (MQ=255)               gTCTCACATCTGACCCT‑GTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGa‑gaaga                                                                                                                                    <  1:346161/149‑1 (MQ=255)                  tcACATCTGACCCT‑GTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGc                                                                                                                                  >  8:76869/1‑148 (MQ=255)                   cacaTCTGACCCT‑GTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGc                                                                                                                                  >  7:159637/1‑147 (MQ=255)                    acttCTGACCCT‑GTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGc                                                                                                                                  <  7:66040/143‑1 (MQ=255)                    acaTCTGACCCT‑GTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGc                                                                                                                                  >  1:9823/1‑146 (MQ=255)                           aCCCT‑GTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGAcccc                                                                                                                        <  6:338770/149‑1 (MQ=255)                           aCCCT‑GTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGAc                                                                                                                           <  8:159637/146‑1 (MQ=255)                                       gACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGAccc                                                                                                             <  1:43672/149‑1 (MQ=255)                                        aCGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGAcccc                                                                                                            >  6:185502/1‑149 (MQ=255)                                          gCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCt                                                                                                           >  2:295320/1‑148 (MQ=255)                                            tAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGa‑gaaga                                                                                                                                    >  7:101091/1‑121 (MQ=255)                                            tAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGa‑gaaga                                                                                                                                    <  8:101091/121‑1 (MQ=255)                                                       ttGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATgg                                                                                              <  2:17665/148‑1 (MQ=255)                                                         gATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCg                                                                                            <  1:231315/148‑1 (MQ=255)                                                             cacaGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTtc                                                                                       >  3:148854/1‑149 (MQ=255)                                                              acaGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTtc                                                                                       >  4:381026/1‑148 (MQ=255)                                                              acaGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTtc                                                                                       >  8:279582/1‑148 (MQ=255)                                                              acaGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCATAAGGCGGttagg                                                                                      <  6:51526/149‑4 (MQ=255)                                                                agagGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTctct                                                                                    <  5:14412/149‑1 (MQ=255)                                                                    ggTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTtctc                                                                                     <  4:57661/144‑1 (MQ=255)                                                                    ggTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTtctc                                                                                     >  3:57661/1‑144 (MQ=255)                                                                    ggTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCt                                                                                <  8:8261/149‑1 (MQ=255)                                                                         gggATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCa                                                                            <  6:15971/148‑1 (MQ=255)                                                                          ggATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGt                                                                          >  4:104252/1‑149 (MQ=255)                                                                            aTAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTct                                                                        <  7:76869/149‑1 (MQ=255)                                                                             tAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTctg                                                                       <  8:295038/149‑1 (MQ=255)                                                                                ccTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCg                                                                                            >  2:174416/1‑125 (MQ=255)                                                                                  ttAAGGCGTAATTGTCTGATGCCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCtgctgc                                                                   >  3:375534/1‑148 (MQ=255)                                                                                     aGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAg                                                               <  2:9823/149‑1 (MQ=255)                                                                                       gCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCa                                                                            >  7:307894/1‑134 (MQ=255)                                                                                       gCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCa                                                                            <  8:307894/134‑1 (MQ=255)                                                                                       gCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGc                                                             >  2:388435/1‑149 (MQ=255)                                                                                            aTTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATa                                                        <  7:369898/149‑1 (MQ=255)                                                                                             ttGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAg                                                       >  1:392942/1‑149 (MQ=255)                                                                                                  tGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGCAGGCCGATAGCCGAg                                                  >  6:271489/1‑149 (MQ=255)                                                                                                    aTACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGCCCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTc                                                >  4:68197/1‑149 (MQ=255)                                                                                                        cAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCTCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGcc                                                     <  4:218338/140‑1 (MQ=255)                                                                                                        cAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCTCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGcc                                                     >  3:218338/1‑140 (MQ=255)                                                                                                        cAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTccaac                                            <  7:201803/149‑1 (MQ=255)                                                                                                          acacCTCTGCGATGAGGACATGACCTCTTCACCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTccaacc                                           >  3:19284/1‑148 (MQ=255)                                                                                                              ctctGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGCTAGCCGAGGCCCACCAATTc                                      >  1:182225/1‑149 (MQ=255)                                                                                                                   cGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCTGATAGCCGAGTCCAACCAATTCTAt                                   <  5:338546/147‑1 (MQ=255)                                                                                                                   cGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTAtt                                  <  3:104252/148‑1 (MQ=255)                                                                                                                     aTGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTAttttt                               >  1:217594/1‑149 (MQ=255)                                                                                                                     aTGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTAttttt                               >  3:72780/1‑149 (MQ=255)                                                                                                                      tGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTc                              >  3:243418/1‑149 (MQ=255)                                                                                                                      tGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTc                              <  2:217594/149‑1 (MQ=255)                                                                                                                        aGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCg                             <  1:295320/148‑1 (MQ=255)                                                                                                                        aGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCg                             <  5:185502/148‑1 (MQ=255)                                                                                                                        aGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGa                            <  4:375534/149‑1 (MQ=255)                                                                                                                              tGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGa                                                   >  4:335551/1‑120 (MQ=255)                                                                                                                              tGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGc                      <  2:392942/149‑1 (MQ=255)                                                                                                                               tggctctTC‑TCGGCGGGCTTAAGGCGCGCCGATGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGcc                     >  8:249695/4‑149 (MQ=38)                                                                                                                                aCCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTg                        <  5:236529/145‑1 (MQ=255)                                                                                                                                aCCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTg                        >  6:236529/1‑145 (MQ=255)                                                                                                                                 ccTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCt                    >  4:103081/1‑148 (MQ=255)                                                                                                                                  ctctTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCACCCAATAATATTTTTCGACCTGGCCTGa                  >  6:30807/1‑149 (MQ=255)                                                                                                                                   tctTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGAc                 >  3:361985/1‑149 (MQ=255)                                                                                                                                      ttctCGGCGGGCTTAAGGCGCGCCGATGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCg               >  7:360367/4‑149 (MQ=255)                                                                                                                                       tctCGGCGGGCTTAAGGCGCGCCGATGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCg               >  3:112802/3‑148 (MQ=255)                                                                                                                                       tctCGGCGGGCTTAAGGCGCGCCGATGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCg               >  3:118277/3‑148 (MQ=255)                                                                                                                                         tCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGt           <  8:68216/149‑1 (MQ=255)                                                                                                                                            gcgcgcCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTaga        <  3:381026/149‑1 (MQ=255)                                                                                                                                              gcgcCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGc      >  3:152855/1‑149 (MQ=255)                                                                                                                                                  cTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCgtgt  >  2:391848/1‑149 (MQ=255)                                                                                                                                                  cTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCgtgt  >  5:397350/1‑149 (MQ=255)                                                                                                                                                  cTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTATCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCgagt  >  4:171660/1‑149 (MQ=255)                                                                                                                                                   ttAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTaga        >  5:393347/1‑142 (MQ=255)                                                                                                                                                   ttAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTaga        <  6:393347/142‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GTACTCCACACCCAGTCTCACATCTGACCCT‑GTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGA‑GAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGT  >  NC_002947/5041054‑5041346

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