Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,414,830	T→G	6.2%	C10G (TGC→GGC)	PP_2115 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,414,830	0	T	G	6.2%	159.2 / 2.2	62	C10G (TGC→GGC)	PP_2115	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (25/32);  new base G (2/2);  total (28/34)
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.71e-01

GCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGC‑‑‑‑GCTTGA‑CGGGC‑ATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGAAAGCCGTGCTTGGACCG  >  NC_002947/2414684‑2414977                                                                                                                                                         |                                                                                                                                                    gCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGC‑‑‑‑GCTTGA‑CGGGC‑ATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGc                                                                                                                                                   <  5:334915/149‑1 (MQ=32)      gcTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGC‑‑‑‑GCTTGA‑CGGGC‑ATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCg                                                                                                                                              <  5:6823/149‑1 (MQ=32)      gcTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGC‑‑‑‑GCTTGA‑CGGGC‑ATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCg                                                                                                                                              <  6:70455/149‑1 (MQ=32)       cTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCCCGCCCAAGGCCCCATCCAGGGGCCACAGC‑‑‑‑GCATGA‑CGGGC‑ATCCATGCCCGTCACCTGCCATCGCAAGGCCTGCGTTAgg                                                                                                                                             >  2:307230/1‑149 (MQ=11)        tCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTACCGGCAGGACCGCGCCCAAGGCCCCATCCTGGGGCCTCAGC‑‑‑‑GCTTGA‑CGGGC‑ATCCATGCCCGTCACCTGCCATCCCAAGGCCAGCGTTCGGc                                                                                                                                            >  8:345666/1‑149 (MQ=11)          ccATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGC‑‑‑‑GCTTGA‑CGGGC‑ATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGc                                                                                                                                            <  4:279396/147‑1 (MQ=32)           catcatCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGC‑‑‑‑GCTTGA‑CGGGC‑ATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGcctc                                                                                                                                         >  6:118798/1‑149 (MQ=32)            atcatcCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGC‑‑‑‑GCTTGA‑CGGGC‑ATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGcctcc                                                                                                                                        >  7:287328/1‑149 (MQ=32)             tcatcCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGC‑‑‑‑GCTTGA‑CGGGC‑ATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGcctcct                                                                                                                                       <  8:221057/149‑1 (MQ=32)                        ctgcgctgcgCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGC‑‑‑‑GCTTGA‑CGGGC‑ATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaa                                                                                                                            <  5:118798/149‑1 (MQ=32)                        ctgcgctgcgCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGC‑‑‑‑GCTTGA‑CGGGC‑ATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaa                                                                                                                            <  3:162363/149‑1 (MQ=32)                                cgctCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGC‑‑‑‑GCTTGA‑CGGGC‑ATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaa                                                                                                                      <  2:14105/147‑1 (MQ=32)                                  ctCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGTGGCCTCAGC‑‑‑‑GCTTGA‑CGGGC‑ATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatc                                                                                                                  <  5:244506/149‑1 (MQ=17)                                   tCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGC‑‑‑‑GCTTGA‑CGGGC‑ATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatc                                                                                                                  >  6:304576/1‑148 (MQ=32)                                     cGGGGTCCCCTCACTCCGGCCTTTCTCCCGGCAGGACCGCGGCGAAGGCCCCATCCTGGGGCCTCAGC‑‑‑‑GCTTGA‑CGGGC‑ATCCATGCCAGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatc                                                                                                                  <  7:40935/146‑1 (MQ=12)                                           tcaggaggccgaacgcaggccttgcggaggcaggtgacgggcatggatgcccgtcaagcgctgaggccccaggatggggccttcggcgcggtcctgccgggaGCGAGGCCGGAGTTCGGCCTCCTGAAGTCGCGaagatcaagatc                                                                                                                  <  8:292436‑M1/44‑1 (MQ=255)                                            cccTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCCGCCTGGGGCCTCAGC‑‑‑‑GCTTGA‑CGGGC‑CTCCATGCCCGGCACCTGCCTTCGCAAGGCCGGCGTTCGGCCTCCTGCAGTCGCGcagagcaaggtcaagttcaag                                                                                                         >  3:212676/1‑148 (MQ=11)                                            cccTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGC‑‑‑‑GCTTGA‑CGGGC‑ATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaaga                                                                                                        >  7:203238/1‑149 (MQ=32)                                              cTCAATCCGGCCTTGCTGCCGGCATGACCGCGCCGAAGGCCCCATTCTGGGGGCTCACCCACGGTTTGAGCCACCTATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatc                                                                                                            <  1:307230/149‑1 (MQ=2)                                                   tCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGC‑‑‑‑GCTTGA‑CGGGC‑ATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatcaaga                                                                                                  >  8:315153/1‑148 (MQ=32)                                                   tCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGC‑‑‑‑GCTTGA‑CGGGC‑ATCCATGCCCGTCACCTGCCTTCACAAGGCCTGCGTTCGGcctcct                                                                                                                                       <  1:133198/111‑1 (MQ=17)                                                   tCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGC‑‑‑‑GCTTGA‑CGGGC‑ATCCATGCCCGTCACCTGCCTTCACAAGGCCTGCGTTCGGcctcct                                                                                                                                       >  2:133198/1‑111 (MQ=17)                                                     cGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGC‑‑‑‑GCTTGA‑CGGGC‑ATCCATTCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatcaagatc                                                                                                <  7:345666/148‑1 (MQ=17)                                                      aacgcaggccttgcggaggcaggtgacgggcatggatgcccgtcaagcgctgaggccccaggatggggccttcggcgcggtcctgccgggaGCGAGGCCGGAGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatcaa                                                                                                    >  7:292436‑M1/92‑149 (MQ=255)                                                         gcaggccttgcggaggcaggtgacgggcatggatgcccgtcaagcgctgaggccccaggatggggccttcggcgcggtcctgccgggaGCGAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaag                                                                                                               >  4:216508‑M1/89‑135 (MQ=255)                                                         gcaggccttgcggaggcaggtgacgggcatggatgcccgtcaagcgctgaggccccaggatggggccttcggcgcggtcctgccgggaGCGAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaag                                                                                                               <  3:216508‑M1/47‑1 (MQ=255)                                                            gCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGC‑‑‑‑GCTTGA‑CGGGC‑ATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatcaagatcaagatca                                                                                         <  4:102622/148‑1 (MQ=32)                                                               cccGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGC‑‑‑‑GCTTGA‑CGGGC‑ATCCCTGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatcaagatcaagatc                                                                                          >  1:350513/1‑144 (MQ=17)                                                                ccGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGC‑‑‑‑GCTTGA‑CGGGC‑ATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatcaagatcaagatcaagatc                                                                                    <  3:218339/149‑1 (MQ=32)                                                                ccGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGC‑‑‑‑GCTTGA‑CGGGC‑ATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatcaagatcaagatcaagatc                                                                                    <  5:304576/149‑1 (MQ=32)                                                                ccGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGC‑‑‑‑GCTTGA‑CGGGC‑ATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatcaagatcaagatcaagat                                                                                     >  8:31099/1‑148 (MQ=32)                                                                 cGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGC‑‑‑‑GCTTGA‑CGGGC‑ATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatcaagatcaagatcaagatca                                                                                   <  8:287328/149‑1 (MQ=32)                                                                         ccGCGCCGAAGGCCCCATCCTGGGGCCTCAGC‑‑‑‑CCTTGA‑CGGGC‑ATACATGCCCGTCACCTGCCTTCGCAACGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaacatcaagaacaagatcaagaacaacatc                                                                              >  2:109977/1‑146 (MQ=11)                                                                          aggtgacgggcatggatgcccgtcaagcgctgaggccccaggatggggccttcggcgcggtcctgccgggaGCAAGGCCGGAGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatcaagatcaagatcaagatcaa                                                                                  <  1:172036‑M1/76‑1 (MQ=255)                                                                            gtgacgggcatggatgcccgtcaagcgctgaggccccaggatggggccttcggcgcggtcctgccgggaGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatcaagatcaagatcaagatcaagat                                                                               <  6:46213‑M1/79‑1 (MQ=255)                                                                                gAAGGCCCCATCCTGGGGCCTCAGC‑‑‑‑GCTTGA‑CGGGC‑ATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatcaagatcaagatcaagatcaagatcaagatcaaga                                                                    <  7:31099/149‑1 (MQ=33)                                                                                      cccATCCTGGGGCCTCAGC‑‑‑‑GCTTGA‑CGGGC‑ATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatcaagatcaagatcaagatcaagatcaagatcaagatc                                                                  >  7:152402/1‑145 (MQ=35)                                                                                        cATCCTGGGGCCTCAGC‑‑‑‑GCTTGA‑CGGGC‑ATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatcaagatcaagatcaagatcaagatcaagatcaagatcaagatc                                                            <  8:203238/149‑1 (MQ=40)                                                                                               gtgCCTCAGC‑‑‑‑GGTTGA‑CGGGC‑ATCCATGCGAGTCACCTGCCTTCGCAAGGCCTGCGTTAGGCCTCCTGAAGTCGCGAAGTTCAAGATCTAGATCAAGATCAAGATCAAGATCTAGATCTAGATCAagatcaagatcaagatcgagatc                                                      <  1:109977/146‑1 (MQ=14)                                                                                                     cAGC‑‑‑‑GCTTGA‑CGGGC‑ATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGGCTCCTGAAGTCCCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAagatcaagatcaagatcgagatcaagata                                                >  1:170364/1‑147 (MQ=255)                                                                                                      aGC‑‑‑‑GCTTGA‑CGGGC‑ATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAagatcaagatcaagatcgagatcaagatc                                                >  2:183142/1‑147 (MQ=255)                                                                                                        c‑‑‑‑gcTTGA‑CGCGC‑ATCCATGCCCGTCACCAGCCGTCGCACCGCCTGCCTTCCGCCTCCAGCAGTCCCCAAGATCAAGATCAAGATCAAGATCAAGAGCAAGATCAACATCAAGATCAagatcaagatcaagatcgagatcaagatcaaga                                            <  6:267523/149‑1 (MQ=17)                                                                                                             gcTTGA‑CGGGC‑ATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatcaagatcaagatcaagatcaagatcaag                                                                           >  4:307740/1‑117 (MQ=32)                                                                                                             gcTTGA‑CGGGC‑ATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatcaagatcaagatcaagatcaagatcaag                                                                           <  3:307740/117‑1 (MQ=32)                                                                                                                 gA‑CGGGC‑ATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAagatcaagatcaagatcgagatcaagatcaagatcg                                         <  8:152402/147‑1 (MQ=255)                                                                                                                  a‑CGGGC‑ATGGATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAggatcaagatcaagatcgagatcaagatcaagatc                                          >  2:42965/1‑145 (MQ=255)                                                                                                                    cGGGC‑ATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGAtcaagatcaagatcgagatcaagatcaagatcgaga                                      <  1:183142/148‑1 (MQ=255)                                                                                                                     gggC‑ATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCaagatcaagatcgagatcaagatcaagatcgagatc                                    <  4:212676/149‑1 (MQ=255)                                                                                                                              aTGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAagatcaagatcaagatcgagatcaagatcaaga                                            >  8:255201/1‑133 (MQ=255)                                                                                                                              aTGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAagatcaagatcaagatcgagatcaagatcaaga                                            <  7:255201/133‑1 (MQ=255)                                                                                                                              aTGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAgcgc                            <  7:315153/149‑1 (MQ=255)                                                                                                                              aTGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAgcgc                            >  8:228664/1‑149 (MQ=255)                                                                                                                              aTGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGAGCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGTTCAAGATCGAGATCAATATCAAGATCGAGATCAAGAgcgc                            >  3:37573/1‑149 (MQ=255)                                                                                                                               tGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCt                           >  1:124542/1‑149 (MQ=255)                                                                                                                                   cGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGt                       <  2:124542/149‑1 (MQ=255)                                                                                                                                       aCCTGCCTTCGCAAGGACTGCTTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGTTCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGa                   <  4:37573/149‑1 (MQ=255)                                                                                                                                         cTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGaaa                 >  2:379218/1‑149 (MQ=255)                                                                                                                                           gCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAagatcaagatcaagatcgaga                                                        >  8:92901/1‑108 (MQ=255)                                                                                                                                           ccgggaGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGaaa                 >  5:46213‑M1/7‑147 (MQ=255)                                                                                                                                            ccTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGAGCAAGATCAAGATCAAGATCACGATGAAGGTCGAGGTGAAGATGAAGATCGAGATCAAGAGCGCGGGGTTTGAAAGcg              >  8:79025/1‑148 (MQ=18)                                                                                                                                                      GCCTGCGTTGCGCCTCCTGAAGTCGCGAAGATTAAGTTCATCATCAAGATCAAGATGAAGCTCAAGTTCAATATCAAGATCAAGTTCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGAAAGCCGTGCTTTGAc    >  2:170364‑M1/1‑149 (MQ=255)                                                                                                                                                        CTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGAAAGCCGTGCTTGGACCg  <  3:123183‑M1/149‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                                    GCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGC‑‑‑‑GCTTGA‑CGGGC‑ATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGAAAGCCGTGCTTGGACCG  >  NC_002947/2414684‑2414977

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