Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,653,199	T→G	6.0%	L362* (TTA→TGA)	PP_3218 →	NtaA/SnaA/SoxA family monooxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,653,199	0	T	G	6.0%	208.7 / 3.6	67	L362* (TTA→TGA)	PP_3218	NtaA/SnaA/SoxA family monooxygenase
Reads supporting (aligned to +/- strand):  ref base T (27/36);  new base G (2/2);  total (29/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

CCAGTACGACCTGGATGCCCCGTTCCCGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTC‑ACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTACCGAAGAAGTGGTACCGTTGCT  >  NC_002947/3653056‑3653346                                                                                                                                                 |                                                                                                                                                    ccAGTACGACCTGGATGCCCCGTTCCCGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTC‑ACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGc                                                                                                                                                >  8:35728/1‑149 (MQ=255)       cGACCTGGATGCCCCGTTCCCGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTC‑ACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTAACTGACGGCCaaact                                                                                                                                          <  3:69482/149‑5 (MQ=255)        gACCTGGATGCCCCGTTCCCGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTC‑ACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGc                                                                                                                                         >  7:150294/1‑149 (MQ=255)         aCCTGGATGCCCCGTTCCCGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTC‑ACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGcc                                                                                                                                        >  7:195114/1‑149 (MQ=255)          ccTGGATGCCCCGTTCCCGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTC‑ACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCt                                                                                                                                       >  1:378016/1‑149 (MQ=255)           cTGGATGCCCCGTTCCCGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTC‑ACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAgg                                                                                                                                            >  3:76500/1‑143 (MQ=255)           cTGGATGCCCCGTTCCCGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTC‑ACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAgg                                                                                                                                            <  4:76500/143‑1 (MQ=255)               aTGCCCCGTTCCCGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTC‑ACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGtt                                                                                                                                  <  4:316763/149‑1 (MQ=255)                tGCCCCGTTCCCGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTC‑ACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGTCCGTCAGTTACGGCCAGGGGCCTTCGTTg                                                                                                                                 >  7:64527/1‑149 (MQ=255)                  gcccGTTCACTCAAGAGAGCCTGTAGTTTGCACGGAGCAGCTTTTTC‑ACGAATGCAAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTgg                                                                                                                                <  8:367876/147‑1 (MQ=255)                   cccGTTCCCGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTC‑ACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGtt                                                                                                                              <  5:103484/149‑1 (MQ=255)                         cccGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTCTCCCAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGc                                                                                                                          <  1:363977/148‑1 (MQ=255)                          ccGCAAGAGACCCTGGAGTTTGCACGGAGCAGCTTTTTC‑AGGAATGCCAAGCGTATTGCGGAGCAGGCGGGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCg                                                                                                                        <  3:270869/148‑1 (MQ=255)                           cGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTC‑ACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGa                                                                                                                       >  2:432977/1‑148 (MQ=255)                           cGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTC‑ACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGa                                                                                                                       >  5:429825/1‑148 (MQ=255)                           cGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTC‑ACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGa                                                                                                                       >  6:292475/1‑148 (MQ=255)                            gCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTC‑ACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGaaa                                                                                                                     <  7:35728/149‑1 (MQ=255)                            gCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTC‑ACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGaaa                                                                                                                     >  8:26107/1‑149 (MQ=255)                             cAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTC‑ACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAAc                                                                                                                    >  3:354561/1‑149 (MQ=255)                                 agagCCTGGAGTTTGCACGGAGCAGCTTTTTC‑ACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTg                                                                                                                <  8:64527/149‑1 (MQ=255)                                      cTGGAGTTTGCACGGAGCAGCTTTTTC‑ACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTg                                                                                                                                 <  6:420600/127‑1 (MQ=255)                                      cTGGAGTTTGCACGGAGCAGCTTTTTC‑ACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTg                                                                                                                                 >  5:420600/1‑127 (MQ=255)                                         gAGTTTGCACGGAGCAGCTTTTTC‑ACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGa                                                                                                         <  1:377007/148‑1 (MQ=255)                                         gAGTTTGCACGGAGCAGCTTTTTC‑ACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGa                                                                                                         <  5:180637/148‑1 (MQ=255)                                            tttGCACGGAGCAGCTTTTTC‑ACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGa                                                                                                      <  6:321237/148‑1 (MQ=255)                                                cACGGAGCAGCTTTTTC‑ACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGcc                                                                                                 <  4:206630/149‑1 (MQ=255)                                                cACGGAGCAGCTTTTTC‑ACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGcc                                                                                                 <  4:199736/149‑1 (MQ=255)                                                cACGGAGCAGCTTTTTC‑ACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGcc                                                                                                 <  2:378016/149‑1 (MQ=255)                                                               tC‑ACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAgggg                                                                                                                                          <  2:424675/93‑1 (MQ=255)                                                               tC‑ACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAgggg                                                                                                                                          >  1:424675/1‑93 (MQ=255)                                                                  aCGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAgcgc                                                                                 <  7:343447/148‑1 (MQ=255)                                                                  aCGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGcgcg                                                                                >  6:371423/1‑149 (MQ=255)                                                                         ccAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTg                                                                         >  5:25531/1‑149 (MQ=255)                                                                          cgaGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTgg                                                                        <  5:407140/147‑1 (MQ=255)                                                                           aaGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGa                                                                       >  8:234133/1‑149 (MQ=255)                                                                           aaGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGa                                                                       >  5:185733/1‑149 (MQ=255)                                                                             gCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATg                                                                     <  6:182605/149‑1 (MQ=255)                                                                             gCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATg                                                                     >  5:182605/1‑149 (MQ=255)                                                                               gTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATggg                                                                   <  8:195114/149‑1 (MQ=255)                                                                                tATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGt                                                                  <  5:371423/149‑1 (MQ=255)                                                                                 aTTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGt                                                                  <  8:16574/148‑1 (MQ=255)                                                                                    gCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCa                                                               <  5:292475/148‑1 (MQ=255)                                                                                        aCCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATc                                                          >  3:368712/1‑149 (MQ=255)                                                                                         ccAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCt                                                         <  1:432977/149‑1 (MQ=255)                                                                                             gcgcGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACAt                                                     >  6:399516/1‑149 (MQ=255)                                                                                             gcgcGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTGGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACAt                                                     <  3:365869/149‑1 (MQ=255)                                                                                                        ggCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTg                                          <  5:399516/149‑1 (MQ=255)                                                                                                          ccTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGc                                        >  5:294824/1‑149 (MQ=255)                                                                                                          ccTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGc                                        <  2:105809/149‑1 (MQ=255)                                                                                                           cTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGc                                        >  2:66735/1‑148 (MQ=255)                                                                                                           cTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGc                                        >  4:180229/1‑148 (MQ=255)                                                                                                              cAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCa                                      <  7:26107/147‑1 (MQ=255)                                                                                                                 tCTGCGTCAGGCTGTTGAGTTTGGCAGTCAGTTACGGCCAGGGGCCTTCGATGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCAAGATcc                                 >  6:191808/1‑149 (MQ=255)                                                                                                                    gCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGa                                                                                               >  6:402773/1‑84 (MQ=255)                                                                                                                    gCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGa                                                                                               <  5:402773/84‑1 (MQ=255)                                                                                                                         aGGCTGTTGAGTTTGGCCGTCAGTTACTGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAatat                                                           <  8:31249/115‑1 (MQ=255)                                                                                                                         aGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGtt                          >  2:35299/1‑148 (MQ=255)                                                                                                                           gCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTa                        <  4:354561/148‑1 (MQ=255)                                                                                                                           gCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTAc                       <  6:185733/149‑1 (MQ=255)                                                                                                                               ttGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCa                                                               <  1:248744/105‑1 (MQ=255)                                                                                                                               ttGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCa                                                               >  2:248744/1‑105 (MQ=255)                                                                                                                                  cccctGGCCGTAACTGACGGCCAGCGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTACCgaag                  >  2:268135/5‑147 (MQ=255)                                                                                                                                   gTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGGTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTACCGaagaag               <  6:25531/149‑1 (MQ=255)                                                                                                                                   ccctGGCCGTAACTGACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTACCgaagaa                <  5:249050/145‑1 (MQ=255)                                                                                                                                   ccctGGCCGTAACTGACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTACCgaag                  >  6:249050/4‑146 (MQ=255)                                                                                                                                         ccGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTACCGAAGAAGTGGTAc         <  6:429825/149‑1 (MQ=255)                                                                                                                                              aGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTACCGAAGAAGTGGTACCGTtg    <  7:234133/149‑1 (MQ=255)                                                                                                                                                ttACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTACCGAAGAAGTGGTACCGTtgct  >  3:56668/1‑149 (MQ=255)                                                                                                                                                ttACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTACCGAAGAAGTGGTACCGTtgct  >  1:181410/1‑149 (MQ=255)                                                                                                                                                ttACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTACCGAAGAAGTGGTACCGTtgct  >  8:371421/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                    CCAGTACGACCTGGATGCCCCGTTCCCGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTC‑ACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTACCGAAGAAGTGGTACCGTTGCT  >  NC_002947/3653056‑3653346

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