Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,469,853	A→T	13.9%	intergenic (+22/+13)	vacJ → / ← PP_2164	lipoprotein VacJ/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,469,853	0	A	T	13.9%	228.0 / 20.7	90	intergenic (+22/+13)	vacJ/PP_2164	lipoprotein VacJ/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (37/40);  new base T (4/9);  total (41/49)
Fisher's exact test for biased strand distribution p-value = 3.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCCGCGCCAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTACCTCGATCTGATC  >  NC_002947/2469706‑2470000                                                                                                                                                      |                                                                                                                                                        acCCGCGCCAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAg                                                                                                                                                        >  4:130927/1‑149 (MQ=255)          aGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGGTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCCGAGGCCGCTTTGCGGc‑‑cccag                                                                                                                                             <  5:301748/149‑5 (MQ=255)          aGCCTGCTCTCGGCAGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCCGAGGCCGCTTTGCGGc‑‑cccag                                                                                                                                             <  7:351981/149‑5 (MQ=255)           gCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCggg                                                                                                                                            <  2:221649/149‑1 (MQ=255)            ccTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGt                                                                                                                                           <  6:35070/149‑1 (MQ=255)             cTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGtt                                                                                                                                          >  6:187031/1‑149 (MQ=255)             cTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGtt                                                                                                                                          <  3:112919/149‑1 (MQ=255)                ctctCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCCGAGGCCGCTTTGCGGC‑‑CCCAGGTTCac                                                                                                                                       <  5:175709/149‑2 (MQ=255)                   tcGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAAt                                                                                                                                      <  5:264365/147‑1 (MQ=255)                   tcGCCCGAAAAGCTGATTCAGGGTGACAAGTACATCCTCATTCGCAATGCCTACCTGCAAAACCTCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCTAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAAt                                                                                                                                      <  3:272812/147‑1 (MQ=255)                     ggCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGCCACGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CGCGGGTTCaa                                                                                                                                       <  1:115962/146‑1 (MQ=255)                     ggCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCCGAGGCCGCTTTGCGGC‑‑CCCAGGTTCacagat                                                                                                                                   <  7:92505/148‑6 (MQ=38)                     ggCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCCGAGGCCGCTTTGCGGC‑‑CCCAGGTTCacagat                                                                                                                                   <  4:60701/148‑6 (MQ=38)                         gAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGc                                                                                                                               >  1:227113/1‑148 (MQ=255)                                 gATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCAt                                                                                                                                  <  1:50645/137‑1 (MQ=255)                                 gATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCAt                                                                                                                                  >  2:50645/1‑137 (MQ=255)                                 gATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCg                                                                                                                      >  4:82351/1‑149 (MQ=255)                                  aTTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATcc                                                                                                                       <  2:279730/147‑1 (MQ=255)                                  aTTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATcc                                                                                                                       <  2:285323/147‑1 (MQ=255)                                      aGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGcc                                                                                                                 >  7:63784/1‑149 (MQ=255)                                        ggTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCg                                                                                                               <  5:296505/149‑1 (MQ=255)                                           gACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCaaa                                                                                                                                                         <  5:395587/106‑1 (MQ=255)                                           gACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCaaa                                                                                                                                                         >  6:395587/1‑106 (MQ=255)                                              aaGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTc                                                                                                         >  2:325529/1‑149 (MQ=255)                                                  aCAACTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAACGTGAACGATG‑‑GCGACGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTctgct                                                                                                     >  3:326756/1‑149 (MQ=255)                                                   cATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTctgct                                                                                                     <  6:143287/148‑1 (MQ=255)                                                    atcatcCTTCGCAATGCCTCCCTGCAAACCCGCGAGTTCAAGGGGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCtgctgc                                                                                                   <  4:212939/149‑1 (MQ=255)                                                    aTCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTctgct                                                                                                     <  6:319281/147‑1 (MQ=255)                                                    aTCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTctgct                                                                                                     >  5:319281/1‑147 (MQ=255)                                                     tCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTgcc                                                                                                  >  4:260736/1‑149 (MQ=255)                                                          aTTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCtgctgc                                                                                                   >  2:275743/1‑143 (MQ=255)                                                          aTTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCtgctgc                                                                                                   <  1:275743/143‑1 (MQ=255)                                                          aTTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTgcccgccc                                                                                             >  4:276937/1‑149 (MQ=255)                                                          aTTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTgcccgccc                                                                                             <  3:82351/149‑1 (MQ=255)                                                           ttCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGcccgcccg                                                                                            >  5:130776/1‑149 (MQ=255)                                                             cGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCgg                                                                                           <  2:365550/148‑1 (MQ=255)                                                              gCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATT‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCATTGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑t                                                                                       <  1:294198/149‑1 (MQ=255)                                                              gCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑t                                                                                       >  1:393510/1‑149 (MQ=255)                                                               cAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑Tg                                                                                      >  5:241006/1‑149 (MQ=255)                                                               cAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑Tg                                                                                      >  3:341531/1‑149 (MQ=255)                                                                 aTGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGt                                                                                    <  4:136010/149‑1 (MQ=255)                                                                     cTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCgg                                                                                 >  8:43160/1‑148 (MQ=255)                                                                        ccTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCa                                                                               <  3:276937/147‑1 (MQ=255)                                                                         cTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATc                                                                            <  4:331411/149‑1 (MQ=255)                                                                           gCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCa                                                                          >  2:236003/1‑149 (MQ=255)                                                                            cAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCAt                                                                         <  2:393510/149‑1 (MQ=255)                                                                              aaaCCGCGAGTTCTTGGTGTAGGAGG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACGTGGGGGGGCATTGCGGC‑‑CGCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGGCCGCCCGGC‑‑TGGTCGGCAATCCATAc                                                                       <  4:326756/149‑1 (MQ=255)                                                                              aaaCCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATAc                                                                       <  8:63784/149‑1 (MQ=255)                                                                              aaaCCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGACAGAATCCGCAGCCCGAACTTCTGCTGCCCGCACGGC‑‑TGGTCGGCAATCCATAc                                                                       >  7:221239/1‑149 (MQ=255)                                                                                      aGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCgg                                                               <  5:187031/149‑1 (MQ=255)                                                                                        ttCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTg                                                             >  4:312133/1‑149 (MQ=255)                                                                                           aaGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCGGCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGc                                                            >  1:120726/1‑149 (MQ=255)                                                                                           aaGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTg                                                          <  1:311536/149‑1 (MQ=255)                                                                                           aaGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTg                                                          <  6:241006/149‑1 (MQ=255)                                                                                           aaGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTg                                                          <  8:312011/149‑1 (MQ=255)                                                                                               tGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCAt                                                      >  4:223068/1‑149 (MQ=255)                                                                                               tGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCAt                                                      <  3:260736/149‑1 (MQ=255)                                                                                                 aaGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATg                                                     <  1:394289/148‑1 (MQ=255)                                                                                                   ggATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAg                                                  >  5:281427/1‑149 (MQ=255)                                                                                                      tG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGccc                                               <  1:236003/149‑1 (MQ=255)                                                                                                               gTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAgc                                        >  6:87414/1‑149 (MQ=255)                                                                                                                tCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAgcg                                       >  3:13500/1‑149 (MQ=255)                                                                                                                tCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAgcg                                       >  8:51534/1‑149 (MQ=255)                                                                                                                  gAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCa                                          >  8:223436/1‑144 (MQ=255)                                                                                                                  gAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCa                                          <  7:223436/144‑1 (MQ=255)                                                                                                                          cTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGAc                             >  2:297177/1‑149 (MQ=255)                                                                                                                          cTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGAc                             >  6:277423/1‑149 (MQ=255)                                                                                                                           ttCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATg                                                                                                                                 >  8:10650/1‑50 (MQ=255)                                                                                                                           ttCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATg                                                                                                                                 <  7:10650/50‑1 (MQ=255)                                                                                                                           ttCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACg                            >  7:992/1‑149 (MQ=255)                                                                                                                           ttCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACg                            >  8:54089/1‑149 (MQ=255)                                                                                                                            tCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGcccc                                              >  1:306178/1‑130 (MQ=255)                                                                                                                            tCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGcccc                                              <  2:306178/130‑1 (MQ=255)                                                                                                                             cTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACg                            >  4:80620/1‑147 (MQ=255)                                                                                                                              catgcattgAACCCGAGGCCGCTTTGCGGA‑‑CCCAGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAg                         >  7:371498/8‑149 (MQ=38)                                                                                                                                aTCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGa                          <  6:130776/146‑1 (MQ=255)                                                                                                                                 tctggtcaCCTGGGGGCCGAAACCGGG‑‑CGCGGGGTCAATGGATTGCCATTAACCGCATCCCGAACTTCTGCCGCCCGCCCGGCCTTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAg                         <  1:325529/141‑1 (MQ=255)                                                                                                                                   attgAACCCGAGGCCGCTTTGCGGC‑‑CCCAGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCg                                    >  1:79502/3‑133 (MQ=38)                                                                                                                                   attgAACCCGAGGCCGCTTTGCGGC‑‑CCCAGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCg                                    <  2:79502/131‑1 (MQ=38)                                                                                                                                   attgAACCCGAGGCCGCTTTGCGGC‑‑CCCAGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGAt                    >  3:24633/3‑149 (MQ=38)                                                                                                                                   attgAACCAGAGGCCGCTTTGCGGA‑‑CCCACGTTCAATCCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCCGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTCGCATGCAGCCCCTTCAGCGCAGGATTGACGGAGTCGAt                    >  8:221269/3‑149 (MQ=255)                                                                                                                                     tgAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTc                       <  4:13500/144‑1 (MQ=255)                                                                                                                                     tgAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATcc                  <  7:43160/149‑1 (MQ=255)                                                                                                                                      gAACCCGAGGCCGCTTTGCGGC‑‑CCCAGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATcc                  <  4:24633/148‑1 (MQ=255)                                                                                                                                        acatGGGGCCGGAATGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCTCAGCGCGAACTTCTGCTGCACGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGt                <  8:221239/145‑1 (MQ=255)                                                                                                                                         ccTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCGATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTAc              <  8:138070/149‑1 (MQ=255)                                                                                                                                          cTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTAcc             <  5:277423/149‑1 (MQ=255)                                                                                                                                            ggggCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCACCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTACCTc           >  4:148741/1‑149 (MQ=255)                                                                                                                                               gCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTACCTCGAt        <  3:130927/149‑1 (MQ=255)                                                                                                                                                 cGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCTGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCAGGATGGACTGAGTCGATCCGGAACTCGAtat      >  2:291673/1‑147 (MQ=255)                                                                                                                                                 cGCAAAGCGGC‑‑CTCGGGTTCAATGCAAGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTACCTCGATCt      <  2:346946/149‑1 (MQ=255)                                                                                                                                                   cAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTACCTCGATCTGa    >  7:151078/1‑149 (MQ=255)                                                                                                                                                     aaGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTACCTCGATCTGATc  <  3:80620/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                        ACCCGCGCCAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTACCTCGATCTGATC  >  NC_002947/2469706‑2470000

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