Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,449,986	T→A	6.0%	P20P (CCT→CCA)	nicA →	nicotinate dehydrogenase subunit A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,449,986	0	T	A	6.0%	205.1 / 2.4	67	P20P (CCT→CCA)	nicA	nicotinate dehydrogenase subunit A
Reads supporting (aligned to +/- strand):  ref base T (32/31);  new base A (2/2);  total (34/33)
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.70e-01

ACTCATTTCATGTAGTGAGTACTTACCGAAACGCTCGCGTCATTTGAGTCCTTTGGTAGCGCAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCC  >  NC_002947/4449839‑4450133                                                                                                                                                    |                                                                                                                                                    actcaTTTCATGTAGTGAGTACTTACCGAAACGCTCGCGTCATTTGAGTCCTTTGGTAGCGCAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACCGCCAGCCCGTCGAAGTCAGCGCCATGCCt                                                                                                                                                     >  3:194827/1‑148 (MQ=255)  ctcaTTTCATGTAGTGAGTACTTACCGAAACGCTCGCGTCATTTGAGTCCTTTGGTAGCGCAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGa                                                                                                                                                   <  2:77532/149‑1 (MQ=255)     aTTTCATGTAGTGAGTACTTACCGAAACGCTCGCGTCATTTGAGTCCTTTGGTAGCGCAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGac                                                                                                                                                  <  1:189995/147‑1 (MQ=255)      tttCATGTAGTGAGTACTTACCGAAACGCTCGCGTCATTTGAGTCCTTTGGTAGCGCAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACAcc                                                                                                                                               >  7:23150/1‑149 (MQ=255)        tCATGTAGTGAGTACTTACCGAAACGCTCGCGTCATTTGAGTCCTTTGGTAGCGCAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACAcccc                                                                                                                                             >  3:69983/1‑149 (MQ=255)          aTGTAGTGAGTACTTACCGAAACGCTCGCGTCATTTGAGTCCTTTGGTAGCGCAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCg                                                                                                                                            <  6:85142/148‑1 (MQ=255)           tGTAGTGAGTACTTACCGAAACGCTCGCGTCATTTGAGTCCTTTGGTAGCGCAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCgct                                                                                                                                          >  4:57207/1‑149 (MQ=255)             tAGTGAGTACTTACCGAAACGCTCGCGTCATTTGAGTCCTTTGGTAGCGCAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCgctgc                                                                                                                                        <  7:251940/149‑1 (MQ=255)                   gTACTTACCGAAACGCTCGCGTCATTTGAGTCCTTTGGTAGCGCAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGa                                                                                                                                  <  4:132458/149‑1 (MQ=255)                      cTTACCGAAACGCTCGCGTCATTTGAGTCCTTTGGTAGCGCAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATc                                                                                                                                <  3:117315/148‑1 (MQ=255)                          ccGAAACGCTCGCGTCATTTGAGTCCTTTGGTAGCGCAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTgc                                                                                                                            <  1:123336/148‑1 (MQ=255)                          ccGAAACGCTCGCGTCATTTGAGTCCTTTGGTAGCGCAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTgc                                                                                                                            <  2:1453/148‑1 (MQ=255)                           cGAAACGCTCGCGTCATTTGAGTCCTTTGGTAGCGCAAGCAACCTCCAGAGATGAGCCAATAATGCCAACACCCATCTCCCAGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGAGCCTgcg                                                                                                                           >  5:62663/1‑148 (MQ=255)                           cGAAACGCTCGCGTCATTTGAGTCCTTTGGTAGCGCAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTgcg                                                                                                                           >  5:289815/1‑148 (MQ=255)                           cGAAACGCTCGCGTCATTTGAGTCCTTTGGTAGCGCAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTgcg                                                                                                                           >  1:212785/1‑148 (MQ=255)                           cGAAACGCTCGCGTCATTTGAGTCCTTTGGTAGCGCAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGAGCCATGCCAGACACCCCGCTGCTGCTGATCCTgcg                                                                                                                           >  3:56355/1‑148 (MQ=255)                                     gcgTCATTTGAGTCCTTTGGTAGCGCAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCtgt                                                                                                                <  6:289815/149‑1 (MQ=255)                                      cgTCATTTGAGTCCTTTGGTAGCGCAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCtgtg                                                                                                               >  6:68591/1‑149 (MQ=255)                                       gTCATTTGAGTCCTTTGGTAGCGCAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGc                                                                                                              <  2:163263/149‑1 (MQ=255)                                             tGAGTCCTTTGGTAGCGCAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGaa                                                                                                         <  2:212785/148‑1 (MQ=255)                                                   cTTTGGTAGCGCAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAAGGACCTGTGCCTGAACGGccc                                                                                                   >  7:254051/1‑148 (MQ=255)                                                    tttGGTAGCGCAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAAc                                                                                                                       <  2:22704/127‑1 (MQ=255)                                                           gcgcAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGaa                                                                                                         >  2:174303/1‑134 (MQ=255)                                                             gcaagcaaCCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCt                                                                                        >  4:63475/1‑149 (MQ=255)                                                              caagcaaCCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTg                                                                                       >  6:92915/1‑149 (MQ=255)                                                              caagcaaCCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTg                                                                                       >  4:260741/1‑149 (MQ=255)                                                                     ccTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTg                                                                                >  7:179554/1‑149 (MQ=255)                                                                        ccAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAAc                                                                                                                       <  4:291816/107‑1 (MQ=255)                                                                        ccAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAAc                                                                                                                       >  3:291816/1‑107 (MQ=255)                                                                         cAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGc                                                                             <  8:179554/148‑1 (MQ=255)                                                                          agagATGAACCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCg                                                                            >  2:87311/1‑148 (MQ=255)                                                                                 aGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCggg                                                                    >  5:34179/1‑149 (MQ=255)                                                                                 aGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCgg                                                                     <  8:23150/148‑1 (MQ=255)                                                                                  gCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCgggg                                                                   <  4:69983/149‑1 (MQ=255)                                                                                   ccAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATTTTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGc                                                                  <  1:87311/149‑1 (MQ=255)                                                                                   ccAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGc                                                                  <  3:63475/149‑1 (MQ=255)                                                                                     aataatGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCgt                                                                >  6:80466/1‑149 (MQ=255)                                                                                      ataatGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCAGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGc                                                                  <  6:2471/146‑1 (MQ=255)                                                                                      ataatGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCAGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGc                                                                  >  5:2471/1‑146 (MQ=255)                                                                                         atGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTa                                                             <  3:260741/148‑1 (MQ=255)                                                                                          tGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGTCCCCAAGTACGGCTGCGGCCTGCGCGAATGAGGGGCGTGTACg                                                           >  8:266018/1‑149 (MQ=255)                                                                                          tGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACg                                                           >  5:154036/1‑149 (MQ=255)                                                                                          tGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACg                                                           >  5:149688/1‑149 (MQ=255)                                                                                           gCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACgg                                                          >  5:19072/1‑149 (MQ=255)                                                                                             aaaCAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGt                                                         >  6:286532/1‑148 (MQ=255)                                                                                             aaaCAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTg                                                        <  6:62663/149‑1 (MQ=255)                                                                                              aacaacCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGa                                                       <  5:80466/149‑1 (MQ=255)                                                                                                caacCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGatc                                                     <  1:188578/149‑1 (MQ=255)                                                                                                 aacCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCTCCATGCCAGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACTGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGatca                                                    <  4:56355/149‑1 (MQ=255)                                                                                                 aacCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGatca                                                    >  5:135156/1‑149 (MQ=255)                                                                                                   ccATCTCCCT‑CAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGc                                                                                                              >  3:3863/1‑88 (MQ=255)                                                                                                   ccATCTCCCT‑CAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGc                                                                                                              <  4:3863/88‑1 (MQ=255)                                                                                                        tcCCTGCAGGTCAACGGCCAGCACGTCGAAGTCAGCGCCATGACTGACCCCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATTCGGGGCGAGTACGGTGAACATCGACgg                                             >  5:307021/1‑149 (MQ=255)                                                                                                           cTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGt                                          <  5:68591/149‑1 (MQ=255)                                                                                                              cAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTg                                         <  6:38904/147‑1 (MQ=255)                                                                                                                  tCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGgccgcc                                   <  6:34179/149‑1 (MQ=255)                                                                                                                   cAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGgccg                                     >  8:288144/1‑146 (MQ=255)                                                                                                                   cAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGgccg                                     <  7:288144/146‑1 (MQ=255)                                                                                                                       ggCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTc                              <  6:135156/149‑1 (MQ=255)                                                                                                                         ccAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCt                            >  7:174545/1‑149 (MQ=255)                                                                                                                         ccAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCt                            >  6:87966/1‑149 (MQ=255)                                                                                                                             cccGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGt                        >  8:91922/1‑149 (MQ=255)                                                                                                                             cccGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGt                        >  1:127913/1‑149 (MQ=255)                                                                                                                                 tCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATc                    <  5:286532/149‑1 (MQ=255)                                                                                                                                    aaGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATcccc                 <  5:92915/149‑1 (MQ=255)                                                                                                                                             ccATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCt                            <  4:170306/129‑1 (MQ=255)                                                                                                                                             ccATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCt                            >  3:170306/1‑129 (MQ=255)                                                                                                                                                   cTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCgccgcc  >  4:90037/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                    ACTCATTTCATGTAGTGAGTACTTACCGAAACGCTCGCGTCATTTGAGTCCTTTGGTAGCGCAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCC  >  NC_002947/4449839‑4450133

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