Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,605,418	A→C	7.0%	Q263P (CAG→CCG) ‡	PP_4077 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,605,418	0	A	C	7.0%	284.1 / 10.8	99	Q263P (CAG→CCG) ‡	PP_4077	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (49/43);  new base C (4/3);  total (53/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.12e-01

TTCAGCCCTGAGGAGGCCTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGC  >  NC_002947/4605270‑4605559                                                                                                                                                       |                                                                                                                                              ttCAGCCCTGAGGAGGCCTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAg                                                                                                                                              <  1:15695/150‑1 (MQ=255)   cgatgggctcggcactgtgcttgaacaggttgacagccggggtgcagaacaggctgaaatcttctgtgccaacgcgaagggtgtctggcagttgacgcgtgaaatggaactcgatgcccacCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGt                                                                                                                                            >  1:1899237‑M1/122‑150 (MQ=255)   cAGCCCTGAGGAGGCCTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGt                                                                                                                                            >  1:1782937/1‑150 (MQ=255)    aGCCCTGAGGAGGCCTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTg                                                                                                                                           <  2:1694010/150‑1 (MQ=255)     gCCCTGAGGAGGCCTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTgg                                                                                                                                          >  2:2182426/1‑150 (MQ=255)      cccTGAGGAGGCCTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTggg                                                                                                                                         >  2:846949/1‑150 (MQ=255)      cccTGAGGAGGCCTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTggg                                                                                                                                         >  2:919022/1‑150 (MQ=255)         tGAGGAGGCCTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACCTGTGGGCag                                                                                                                                      <  1:368741/150‑2 (MQ=255)            ggaggCCTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGa                                                                                                                                   <  1:133131/150‑1 (MQ=255)              gcactgtgcttgaacaggttgacagccggggtgcagaacaggctgaaatcttctgtgccaacgcgaagggtgtctggcagttgacgcgtgaaatggaactcgatgcccacCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTgg                                                                                                                                          >  2:221472‑M1/111‑141 (MQ=255)              gcactgtgcttgaacaggttgacagccggggtgcagaacaggctgaaatcttctgtgccaacgcgaagggtgtctggcagttgacgcgtgaaatggaactcgatgcccacCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTgg                                                                                                                                          <  1:221472‑M1/31‑1 (MQ=255)              aggCCTTGTTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGt                                                                                                                                 <  2:1160897/150‑1 (MQ=255)               ggCCTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGtt                                                                                                                                <  2:1100539/150‑1 (MQ=255)               ggCCTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGtt                                                                                                                                >  1:1990063/1‑150 (MQ=255)                gCCTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATc                                                                                                                                     <  2:65146/144‑1 (MQ=255)                gCCTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATc                                                                                                                                     >  1:65146/1‑144 (MQ=255)                 ccTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTcc                                                                                                                              >  2:1295021/1‑150 (MQ=255)                  cTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCa                                                                                                                             >  1:274501/1‑150 (MQ=255)                     gctgcCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CT‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCAttt                                                                                                                          >  2:1837559/1‑150 (MQ=255)                     gctgcCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCAttt                                                                                                                          <  1:846949/150‑1 (MQ=255)                            tACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGt                                                                                                                   >  2:1184053/1‑150 (MQ=255)                             ccccGCAGAACGCGTTCCATGGCTACCGTATTCTGCAGGAATCCTTCGTGTTCCCGAAGCGCTTCCACTTCTTCCGCATTACCGGCCAGGAAAAGCTGTGG‑CC‑GGCACCAGCCTGCCCACAGTTGGGCATCCAGTTCCATTCCACGCGTc                                                                                                                  <  2:91632/149‑1 (MQ=255)                                cGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAAc                                                                                                               <  2:1990063/150‑1 (MQ=255)                                   aGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGc                                                                                                            <  2:2184142/150‑1 (MQ=255)                                     aaCGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCGTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCa                                                                                                          >  1:804786/1‑150 (MQ=255)                                      aCGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAg                                                                                                         <  1:1477866/150‑1 (MQ=255)                                       cGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGa                                                                                                        <  2:1875867/150‑1 (MQ=255)                                       cGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGa                                                                                                        <  2:2301779/150‑1 (MQ=255)                                           ttCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACAcc                                                                                                    >  2:898937/1‑150 (MQ=255)                                              gATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGGAAAGCGGTGG‑CC‑GGCACAAGCCTGCCCACAGGTgg                                                                                                                                          >  1:1498603/1‑109 (MQ=255)                                              gATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTgg                                                                                                                                          <  2:1498603/109‑1 (MQ=255)                                               aTGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTc                                                                                                <  1:1630670/150‑1 (MQ=255)                                                 ggCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTcgc                                                                                              <  2:994936/150‑1 (MQ=255)                                                    tACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACCTGTGGGCagg                                                                                                                                     >  2:625428/1‑106 (MQ=255)                                                    tACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACCTGTGGGCagg                                                                                                                                     <  1:625428/108‑3 (MQ=255)                                                     aCCGTATTCTGCAGGAATACTTCGTGTGCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTg                                                                                          >  2:417836/1‑150 (MQ=255)                                                      ccGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGt                                                                                                                                            >  2:937469/1‑99 (MQ=255)                                                      ccGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGt                                                                                                                                            <  1:937469/99‑1 (MQ=255)                                                      ccGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTgg                                                                                         <  1:1184053/150‑1 (MQ=255)                                                      ccGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTgg                                                                                         <  1:2182426/150‑1 (MQ=255)                                                         tATTCTGCAGGAATACTTCGTGTTCCCTAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGcac                                                                                      >  2:2376280/1‑150 (MQ=255)                                                         tATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑TGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGcac                                                                                      <  1:898937/150‑1 (MQ=255)                                                         tATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGcac                                                                                      <  2:658151/150‑1 (MQ=255)                                                           ttCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACag                                                                                    >  1:2253902/1‑150 (MQ=255)                                                            tCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACaga                                                                                   >  1:1427294/1‑150 (MQ=255)                                                             cTGCAGGAATACTTCGTGTTCCCTAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACagaa                                                                                  <  1:2376280/150‑1 (MQ=255)                                                              tGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACagaag                                                                                 >  2:1760705/1‑150 (MQ=255)                                                               gCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACagaaga                                                                                >  1:275742/1‑150 (MQ=255)                                                                ttctgtgccaacgcgaagggtgtctggcagttgacgcgtgaaatggaactcgatgcccacCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTgg                                                                                                                                          >  2:283223‑M1/61‑91 (MQ=255)                                                                ttctgtgccaacgcgaagggtgtctggcagttgacgcgtgaaatggaactcgatgcccacCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTgg                                                                                                                                          <  1:283223‑M1/31‑1 (MQ=255)                                                                  ggAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAACGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGAttt                                                                             >  2:936026/1‑150 (MQ=255)                                                                  ggAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACagaa                                                                                  <  2:2220213/145‑1 (MQ=255)                                                                  ggAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACagaa                                                                                  >  1:2220213/1‑145 (MQ=255)                                                                         ttCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACCTGTGGGCag                                                                                                                                      <  2:2407220/86‑2 (MQ=255)                                                                         ttCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACCTGTGGGCag                                                                                                                                      >  1:2407220/1‑85 (MQ=255)                                                                         ttCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCa                                                                                                          <  2:2271507/114‑1 (MQ=255)                                                                         ttCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCa                                                                                                          >  1:2271507/1‑114 (MQ=255)                                                                            gtgtTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTc                                                                   <  1:919022/150‑1 (MQ=255)                                                                             tgtTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGGACAAGCCTGCCCACACGTGGGCATCGAGTTCCATTTCACGCGTAAACTGCCCGACACCCTTCGCGTTGGCACACAAGATTTCAGCCTGTTCt                                                                  >  2:1666837/1‑150 (MQ=255)                                                                               ttCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGc                                                                                                            >  2:2207306/1‑106 (MQ=255)                                                                               ttCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGc                                                                                                            <  1:2207306/106‑1 (MQ=255)                                                                                 ctcGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAACGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCAc                                                              >  2:653188/3‑150 (MQ=255)                                                                                 cccGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCAc                                                              >  2:998568/1‑150 (MQ=255)                                                                                  tcGAAGCGCGTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCAcc                                                             >  2:1705272/2‑150 (MQ=255)                                                                                  ccGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATc                                                                                                                                     <  2:1817144/78‑1 (MQ=255)                                                                                  ccGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATc                                                                                                                                     >  1:1817144/1‑78 (MQ=255)                                                                                  ccGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCa                                                                           >  2:1577129/1‑136 (MQ=255)                                                                                  ccGAAGCGCTTCCACTGCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCa                                                                           <  1:1577129/136‑1 (MQ=255)                                                                                    gAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑AC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACCCCCTTCGCGTTGGCACAGAAGATTCCAGCCTGTTCTCCGcccc                                                           >  2:516691/1‑150 (MQ=255)                                                                                       gcgcTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGc                                                        <  2:1913152/150‑1 (MQ=255)                                                                                       gcgcTTCAACTTCTTCATCATTACGGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCGCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACGCCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGc                                                        <  2:1701869/150‑1 (MQ=255)                                                                                          cTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGt                                                     >  1:772813/1‑150 (MQ=255)                                                                                          cTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGt                                                     >  2:1183395/1‑150 (MQ=255)                                                                                             ccACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAAGTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTTCACCCCGGCTGTCaa                                                  >  2:1522288/1‑150 (MQ=255)                                                                                             ccACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTg                                                                                          >  2:2137236/1‑110 (MQ=255)                                                                                             ccACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTg                                                                                          <  1:2137236/110‑1 (MQ=255)                                                                                              cACTTCTTCAGCATTACCGGTCTGGAAAAGCTTTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAAc                                                 <  1:653188/150‑1 (MQ=255)                                                                                              cACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAAc                                                 <  2:1674658/150‑1 (MQ=255)                                                                                              cACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAAc                                                 <  2:1389902/150‑1 (MQ=255)                                                                                               aCTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAAcc                                                >  1:473570/1‑150 (MQ=255)                                                                                                cttcttCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCt                                               >  1:372325/1‑150 (MQ=255)                                                                                                  tcttcAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGt                                             <  1:1295021/150‑1 (MQ=255)                                                                                                      cAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCaa                                         <  2:1427294/150‑1 (MQ=255)                                                                                                          atggaactcgatgcccacCTGTGGGCA‑GGCTTGTGCCGGCC‑ACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAg                                        <  2:626055‑M1/129‑1 (MQ=255)                                                                                                          atggaactcgatgcccacCTGTGG‑GCAGGCTTGTGCCGGCC‑ACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAg                                        >  1:626055‑M1/19‑147 (MQ=255)                                                                                                                  ccTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGa                             <  1:998568/150‑1 (MQ=255)                                                                                                                   cTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAg                            >  1:2007973/1‑150 (MQ=255)                                                                                                                    tGGAATAGCGGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTGCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGc                           >  2:1299602/1‑150 (MQ=255)                                                                                                                         aaGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATc                      <  1:1561394/150‑1 (MQ=255)                                                                                                                            CTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGAc                   <  2:2398381‑M1/150‑1 (MQ=255)                                                                                                                               TGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTg                >  2:1494607‑M1/1‑150 (MQ=255)                                                                                                                                 g‑CC‑AGCACAAGCCTGCCCACAGGGGGACATCGAGTTCTATTTCACGCGTCAACTGCCAGACACCCTTCCCGATCGCACAGACGACTTCTGACGAATCGCCACCCCCGCTCCCCACTCCGTCAAGCACAGTCCCGACCCCCTCAACCTgac              >  2:572620/1‑148 (MQ=255)                                                                                                                                         acaAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGa        >  1:1553996/1‑150 (MQ=255)                                                                                                                                          caAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAg       >  1:2085712/1‑150 (MQ=255)                                                                                                                                           ctcgatGCCCACCTGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAgg      >  2:1111666/6‑150 (MQ=255)                                                                                                                                           aaGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAgg      <  2:1553996/150‑1 (MQ=255)                                                                                                                                           aaGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAgg      >  2:843781/1‑150 (MQ=255)                                                                                                                                               cTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGc  >  1:377924/1‑150 (MQ=255)                                                                                                                                               atGCCCACCTGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGc  >  1:1859498/2‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                              TTCAGCCCTGAGGAGGCCTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGG‑CC‑GGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGC  >  NC_002947/4605270‑4605559

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑