Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	807,888	A→T	5.6%	L49Q (CTG→CAG)	PP_0694 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	807,888	0	A	T	5.6%	268.8 / 4.3	89	L49Q (CTG→CAG)	PP_0694	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (43/40);  new base T (3/2);  total (46/43)
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.16e-01

ATCGGGTCGCGGCCCTGGTCAGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATG  >  NC_002947/807741‑808030                                                                                                                                                     |                                                                                                                                                aTCGGGTCGCGGCCCTGGTCAGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAg                                                                                                                                                >  6:66872/1‑149 (MQ=255)      gTCGCGGCCCTGGTCAGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCt                                                                                                                                           <  2:184041/149‑1 (MQ=255)       tCGCGGCCCTGGTCAGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCtt                                                                                                                                          <  6:107287/149‑1 (MQ=255)        cgcgGCCCTGGTCAGCAAGGTGCTCGACCAGGATGATGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGAAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTg                                                                                                                                         <  7:90370/149‑1 (MQ=255)             cccTGGTCAGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAg                                                                                                                                    >  6:170569/1‑149 (MQ=255)                tGGTCAGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCCACAGGCGCTGCGCGGTCTTGATGACCTGCGCGGCGGTCATGTCCCGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGATCAGG‑Gt                                                                                                                                >  7:165315/1‑149 (MQ=255)                     aGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCgg                                                                                                                            >  6:4389/1‑148 (MQ=255)                      gCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTa                                                                                                                          >  2:310399/1‑149 (MQ=255)                        aaGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGTACAGG‑GTGCGGTACt                                                                                                                        <  7:349652/149‑1 (MQ=255)                        aaGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGCTCATGTCCGGCACGATGTAGTAGAAACCACGGGTATCGCACACAAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACt                                                                                                                        >  6:366928/1‑149 (MQ=255)                          ggTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTgg                                                                                                                      >  3:319049/1‑149 (MQ=255)                             gCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTcc                                                                                                                   >  1:151453/1‑149 (MQ=255)                              cTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCt                                                                                                                  <  4:61754/149‑1 (MQ=255)                               tCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACAACAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTg                                                                                                                 >  5:183270/1‑149 (MQ=255)                                 gACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCt                                                                                                                  >  8:252142/1‑146 (MQ=255)                                 gACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCt                                                                                                                  <  7:252142/146‑1 (MQ=255)                                  aCCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTgg                                                                                                                >  3:32476/1‑147 (MQ=255)                                  aCCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCAC‑CCTGTTCA‑AGGCCCTGGGCTTCGAACTGGTGTGCGATAccc                                                                                                                       <  3:223967/140‑3 (MQ=25)                                  aCCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCAC‑CCTGTTCA‑AGGCCCTGGGCTTCGAACTGGTGTGCGATAccc                                                                                                                       >  4:223967/1‑138 (MQ=25)                                     aGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGtt                                                                                                           <  4:32476/149‑1 (MQ=255)                                        aTGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGa                                                                                                        <  2:116847/149‑1 (MQ=255)                                         tGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGa                                                                                                        >  2:336574/1‑148 (MQ=255)                                         tGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACAAGTTCG‑AAGCCCAGGCCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGa                                                                                                        >  5:194476/1‑148 (MQ=255)                                                 aaCAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTggg                                                                                                <  4:319049/148‑1 (MQ=255)                                                 aaCAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGc                                                                                               >  1:355153/1‑149 (MQ=255)                                                 aaCAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGc                                                                                               <  5:170569/149‑1 (MQ=255)                                                  aCAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGc                                                                                               <  1:171383/148‑1 (MQ=255)                                                     gCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTAc                                                                                           <  6:265899/149‑1 (MQ=255)                                                      cGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTa                                                                                            >  3:142164/1‑147 (MQ=255)                                                            gcgcTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGTCAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCa                                                                                     <  1:310399/149‑1 (MQ=255)                                                             cgcTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAgg                                                                                   <  2:151453/149‑1 (MQ=255)                                                             cgcTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAgg                                                                                   <  8:76081/149‑1 (MQ=255)                                                                tGCGCGGTCTCGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCGGGGTc                                                                                >  3:253965/1‑149 (MQ=255)                                                                  cgcgGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTcgc                                                                              <  2:84123/149‑1 (MQ=255)                                                                   ggggTGTTGTTGACCTGCGCGGGGGTCATGTCCGGCACGAAGTATTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTcgcg                                                                             <  6:194476/147‑1 (MQ=255)                                                                     ggTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCAGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGcc                                                                           >  2:143636/1‑149 (MQ=255)                                                                         ttgttgACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCt                                                                       >  7:10279/1‑149 (MQ=255)                                                                         gacatgACCGCCGCGCAGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCt                                                                       >  4:93474/5‑149 (MQ=255)                                                                         gacatgACCGCCGCGCAGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCt                                                                       >  2:221444/5‑149 (MQ=255)                                                                               aCCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTg                                                                 <  1:336574/149‑1 (MQ=255)                                                                               aCCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTAG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGCAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTg                                                                 >  6:289396/1‑149 (MQ=255)                                                                                 cTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTgg                                                                >  8:67872/1‑148 (MQ=255)                                                                                 cTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTgg                                                                >  3:35609/1‑148 (MQ=255)                                                                                    cgcgGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAg                                                             <  5:66872/148‑1 (MQ=255)                                                                                      cggcggTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCt                                                          >  8:270216/1‑149 (MQ=255)                                                                                        gcggTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCttt                                                        <  2:138920/149‑1 (MQ=255)                                                                                        gcggTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCttt                                                        <  7:67872/149‑1 (MQ=255)                                                                                         cggTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCGGGGTCGCGCCGGCTGACGTGGAAGCCTTTg                                                       <  4:253965/149‑1 (MQ=255)                                                                                         cggTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GGGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTg                                                       <  4:51058/149‑1 (MQ=255)                                                                                            tCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGaa                                                     >  2:60983/1‑148 (MQ=255)                                                                                               tGTCCGGAACTAAGTAGTAGAAACCCCGGGTATCGCACACCCGTTCG‑AAGCCCCGGGCCTTGAACAGG‑GGGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGc                                                 <  8:28121/149‑1 (MQ=255)                                                                                                gTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGc                                                 >  4:97630/1‑148 (MQ=255)                                                                                                 tCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCttt                                                        <  5:87628/140‑1 (MQ=255)                                                                                                 tCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCttt                                                        >  6:87628/1‑140 (MQ=255)                                                                                                   cGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAgg                                                                                                                                   >  5:26893/1‑64 (MQ=255)                                                                                                   cGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTcgc                                             <  1:60983/149‑1 (MQ=255)                                                                                                    ggCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTTACGTGGAAGCCTTTGAACAGCTcgcg                                            <  6:183270/149‑1 (MQ=255)                                                                                                      cACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCgg                                           >  7:284163/1‑148 (MQ=255)                                                                                                      cACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCgg                                           >  2:286983/1‑148 (MQ=255)                                                                                                      cACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGa                                          <  5:4389/149‑1 (MQ=255)                                                                                                         gAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCgg                                           >  2:141098/1‑145 (MQ=255)                                                                                                              agtagAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCgg                                   <  4:35609/148‑1 (MQ=255)                                                                                                               gtagAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGgcg                                 >  7:220357/1‑149 (MQ=255)                                                                                                               gtagAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGgc                                  <  7:270216/148‑1 (MQ=255)                                                                                                                tagAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGgcgc                                >  1:149557/1‑149 (MQ=255)                                                                                                                  gAAACCACGGGTATCGCACACCGGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTGGCGGAAGATCGgcgc                                >  7:59424/1‑147 (MQ=255)                                                                                                                  gAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGgcgc                                >  6:33613/1‑147 (MQ=255)                                                                                                                     aCCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCt                           >  2:342658/1‑149 (MQ=255)                                                                                                                      ccACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCTGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTg                          <  1:286983/149‑1 (MQ=255)                                                                                                                       cACGGGCATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGGCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTgg                         <  8:158044/149‑1 (MQ=255)                                                                                                                          cagTACCGCAC‑CCTGTTCA‑AGGCCCTGGGCTTCGAACTGGTGTGCGATACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGa                       >  1:56034/3‑148 (MQ=25)                                                                                                                           agTACCGCAC‑CCTGTTCA‑AGGCCCTGGGCTTCGAACTGGTGTGCGATACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGa                       >  3:95580/2‑147 (MQ=25)                                                                                                                             tATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGt                   >  5:21947/1‑149 (MQ=255)                                                                                                                               tCGCACACCATAGCG‑TACCCCAGGGCCTTGCACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACACCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGtt                  <  8:165315/148‑1 (MQ=255)                                                                                                                               tCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGTCAGtt                  <  8:284163/148‑1 (MQ=255)                                                                                                                                 gCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTc                 <  1:82327/147‑1 (MQ=255)                                                                                                                                 gCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTc                 >  2:82327/1‑147 (MQ=255)                                                                                                                                 gCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCgg               <  1:141098/149‑1 (MQ=255)                                                                                                                                  cacaCCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGa              <  3:97630/149‑1 (MQ=255)                                                                                                                                     ccctgTTCA‑AGGCCCTGGGCTTCGAACTGGTGTGCGATACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGaaa            <  4:95580/145‑1 (MQ=37)                                                                                                                                     acCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGTAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAg           <  8:220357/149‑1 (MQ=255)                                                                                                                                     acCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGa              >  3:37701/1‑146 (MQ=255)                                                                                                                                      ccAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGa          <  2:355153/149‑1 (MQ=255)                                                                                                                                       cAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTCAACAGCTCGCGGAAGATCGGCGCGAGCTGGTACAGTTCGGAAAGa          >  4:49877/1‑148 (MQ=255)                                                                                                                                          ttCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCaa      <  5:33613/149‑1 (MQ=255)                                                                                                                                             g‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCAGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGac   >  3:128662/1‑148 (MQ=255)                                                                                                                                               aaGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATg  >  8:50674/1‑149 (MQ=255)                                                                                                                                                   ccAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAg     >  8:190377/1‑142 (MQ=255)                                                                                                                                                   ccAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAg     <  7:190377/142‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                ATCGGGTCGCGGCCCTGGTCAGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATG  >  NC_002947/807741‑808030

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