Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	807,881	G→A	5.5%	F51F (TTC→TTT) ‡	PP_0694 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	807,881	0	G	A	5.5%	279.6 / 3.8	90	F51F (TTC→TTT) ‡	PP_0694	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (43/42);  new base A (3/2);  total (46/44)
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.20e-01

GCACGGCCATCGGGTCGCGGCCCTGGTCAGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGAT  >  NC_002947/807733‑808029                                                                                                                                                     |                                                                                                                                                       gCACGGCCATCGGGTCGCGGCCCTGGTCAGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTACTAGAAACCACGGGTATCGCACACCAGTTCg                                                                                                                                                        <  1:33892/149‑1 (MQ=255)    cGGCCATCGGGTCGCGGCCCTGGTCAGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAg                                                                                                                                                    <  4:135930/149‑1 (MQ=255)       ccATCGGGTCGCGGCCCTGGTCAGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGcc                                                                                                                                                  >  6:180588/1‑148 (MQ=255)       ccATCGGGTCGCGGCCCTGGTCAGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGcc                                                                                                                                                  >  7:76081/1‑148 (MQ=255)         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                                                                                                                 <  7:252142/146‑1 (MQ=255)                                         gACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCt                                                                                                                 >  8:252142/1‑146 (MQ=255)                                          aCCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTgg                                                                                                               >  3:32476/1‑147 (MQ=255)                                          aCCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCAC‑CCTGTTCA‑AGGCCCTGGGCTTCGAACTGGTGTGCGATAccc                                                                                                                      >  4:223967/1‑138 (MQ=25)                                          aCCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCAC‑CCTGTTCA‑AGGCCCTGGGCTTCGAACTGGTGTGCGATAccc                                                                                                                      <  3:223967/140‑3 (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‑a                                                                                                                                                      <  5:384448/93‑1 (MQ=255)                                                          aCAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑a                                                                                                                                                      >  6:384448/1‑93 (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                                                       >  6:87628/1‑140 (MQ=255)                                                                                                         tCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCttt                                                       <  5:87628/140‑1 (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)                                                                                                                                                     |                                                                                                                                                       GCACGGCCATCGGGTCGCGGCCCTGGTCAGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGAT  >  NC_002947/807733‑808029

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