Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	951,217	A→T	5.0%	I16I (ATA→ATT)	cyoB →	cytochrome bo terminal oxidase subunit I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	951,217	0	A	T	5.0%	251.0 / 2.0	79	I16I (ATA→ATT)	cyoB	cytochrome bo terminal oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base A (36/39);  new base T (2/2);  total (38/41)
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.99e-01

CGCGGTCGTAGCCACGAAGAAGCAGGCAGCAAAGATCTGGCCACAACCAAGGGTGTGGAATCGAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCG  >  NC_002947/951071‑951364                                                                                                                                                    |                                                                                                                                                     cgcgGTCGTAGCCACGAAGAAGCAGGCAGCAAAGATCTGGCCACAACCAAGGGTGTGGAATCGAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGt                                                                                                                                                    >  5:57782/1‑149 (MQ=255)    gggcGTAGCCACGAAGAAGCAGGCAGCAAAGATCTGGCCACAACCAAGGGTGTGGAATCGAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCAt                                                                                                                                                 <  8:195228/146‑1 (MQ=255)    ggTCGTAGCCACGAAGAAGCAGGCAGCAAAGATCTGGCCACAACCAAGGGTGTGGAATCGAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCAt                                                                                                                                                 <  3:230140/149‑1 (MQ=255)     gTCGTAGCCACGAAGAAGCAGGCAGCAAAGATCTGGCCACAACCAAGGGTGTGGAATCTAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATg                                                                                                                                                <  5:54239/149‑1 (MQ=255)      tCGTAGCCACGAAGAAGCAGGCAGCAAAGATGTGGCCACAACCAAGGGTGTGGAATCGAGTATGCAACCTGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATgg                                                                                                                                               <  5:121740/149‑1 (MQ=255)        gTAGCCACGAAGAAGCAGGCAGCAAAGATCTGGCCACAACCAAGGGTGTGGAATCGAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATTGTCATGGt                                                                                                                                              >  4:49172/1‑148 (MQ=255)           gCCACGAAGAAGCAGGCAGCAAAGATCTGGCCACAACAAAGTGTGTGGAATCTAGTATGCAACCAGCTGCCGGCGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACg                                                                                                                                          <  3:401012/149‑1 (MQ=255)            ccACGAAGAAGCAGGCAGCAAAGATCTGGCCACAACCAAGGGTGTGGAATCGAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGc                                                                                                                                         <  5:411814/149‑1 (MQ=255)                 aagaagCAGGCAGCAAAGATCTGGCCACAACCAAGGGTGTGGAATCGAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGcc                                                                                                                                    <  5:152584/149‑1 (MQ=255)                 aagaagCAGGCAGCAAAGATCTGGCCACAACCAAGGGTGTGGAATCGAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGc                                                                                                                                     >  1:73547/1‑148 (MQ=255)                      gcaggcagCAAAGATCTGGCCACAACCAAGGGTGTGGAATCGAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGAt                                                                                                                               >  8:282482/1‑149 (MQ=255)                       caggcagCAAAGATCTGGCCACAACCAAGGGTGTGGAATCGAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGTTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATc                                                                                                                              <  6:240694/149‑1 (MQ=255)                          gcagcaAAGATCTGGCCACAACCAAGGGTGTGGAATCGAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATcgcg                                                                                                                           >  2:119966/1‑149 (MQ=255)                                     cTGGCCACAACCAAGGGTGTGGAATCGAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTAt                                                                                                                >  3:113750/1‑149 (MQ=255)                                      tGGCCACAACCAAGGGTGTGGAATCGAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATc                                                                                                               <  6:441843/149‑1 (MQ=255)                                        gCCACAACCAAGGGTGTGGAATCGAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCg                                                                                                              <  6:57782/148‑1 (MQ=255)                                          cacaACCAAGGGTGTGGAATCGAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTg                                                                                                           >  7:387054/1‑149 (MQ=255)                                          cacaACCAAGGGTGTGGAATCGAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTg                                                                                                           <  4:113750/149‑1 (MQ=255)                                           acaACCAAGGGTGTGGAATCGAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTATCCACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTg                                                                                                           <  1:119966/149‑1 (MQ=255)                                           acaACCAAGGGTGTGGAATCGAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTgt                                                                                                          <  1:31773/149‑1 (MQ=255)                                             aaCCAAGGGTGTGGAATCGAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTcggc                                                                                                                     >  2:238987/1‑136 (MQ=255)                                             aaCCAAGGGTGTGGAATCGAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTcggc                                                                                                                     <  1:238987/136‑1 (MQ=255)                                                cAAGGGTGTGGAATCGAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGtcgt‑cg                                                                                                    <  4:283077/149‑1 (MQ=255)                                                cAAGGGTGTGGAATCGAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGtcgt‑cg                                                                                                    >  4:52819/1‑149 (MQ=255)                                                cAAGGGTGTGGAATCGAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGtcgt‑cg                                                                                                    >  1:360555/1‑149 (MQ=255)                                                   gggTGTGGAATCGAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTc                                                                                                 <  3:70264/149‑1 (MQ=255)                                                        tgGAATCGAGTATGCAACCAGCTGCCGGTGAAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTAt                                                                                             >  8:266647/1‑148 (MQ=255)                                                         ggAATCGAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCa                                                                                           >  3:167329/1‑149 (MQ=255)                                                          gAATCGAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCAc                                                                                          <  4:167329/149‑1 (MQ=255)                                                            aTCGAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTg                                                                                                                                      <  3:350069/104‑1 (MQ=255)                                                            aTCGAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTg                                                                                                                                      >  4:350069/1‑104 (MQ=255)                                                             tCGAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTa                                                                                       >  2:33112/1‑149 (MQ=255)                                                             tCGAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTa                                                                                       >  6:439015/1‑149 (MQ=255)                                                               gAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATc                                                                                                                              >  4:303973/1‑109 (MQ=255)                                                               gAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTAtt                                                                                     >  1:386914/1‑149 (MQ=255)                                                                  tATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTcc                                                                                  >  6:419794/1‑149 (MQ=255)                                                                    tGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCg                                                                                 <  8:387054/148‑1 (MQ=255)                                                                        aCCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGTTACCCTCT‑CACGACCGGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGt                                                                            <  1:142761/149‑1 (MQ=255)                                                                        aCCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGt                                                                            <  2:74137/149‑1 (MQ=255)                                                                        aCCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGCCGT‑CGGTCTTATCACCTATTTCCGCAAGt                                                                            <  1:33112/149‑1 (MQ=255)                                                                            gCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGTGGAc                                                                        >  8:243050/1‑149 (MQ=255)                                                                             cTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGTGGAcc                                                                       <  5:439015/149‑1 (MQ=255)                                                                             cTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGTGGAcc                                                                       >  4:341964/1‑149 (MQ=255)                                                                                ccGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGTGGACCTAc                                                                    >  3:422152/1‑149 (MQ=255)                                                                                ccGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGTGGACCTAc                                                                    >  1:343722/1‑149 (MQ=255)                                                                                      cAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGTGGACCTACTtgtg                                                               <  2:303071/148‑1 (MQ=255)                                                                                      cAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGTGGACCTACTtgtg                                                               >  2:161189/1‑148 (MQ=255)                                                                                       agagGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGa                                                             <  3:341964/149‑1 (MQ=255)                                                                                           gagTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACTCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGa                                                         <  7:282482/149‑1 (MQ=255)                                                                                              tAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTg                                                      >  5:7109/1‑149 (MQ=255)                                                                                              tAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTg                                                      >  3:333697/1‑149 (MQ=255)                                                                                              tAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTg                                                      >  1:389306/1‑149 (MQ=255)                                                                                              tAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGAGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCAGCAAGTGGACCTACTTGTGGAGCGAGtt                                                      >  3:37470/1‑148 (MQ=255)                                                                                                agagATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATTGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGc                                                    >  3:427236/1‑149 (MQ=255)                                                                                                   gaTGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGa                                                 <  2:360555/149‑1 (MQ=255)                                                                                                    aTGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGTgg                                                                          >  5:34885/1‑123 (MQ=255)                                                                                                    aTGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGTgg                                                                          <  6:34885/123‑1 (MQ=255)                                                                                                       ttCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATTGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGacta                                              <  3:49172/148‑1 (MQ=255)                                                                                                          ggTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGt                                          <  8:146423/149‑1 (MQ=255)                                                                                                           gTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATTGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTc                                         >  8:205927/1‑149 (MQ=255)                                                                                                           gTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTc                                         <  3:52819/149‑1 (MQ=255)                                                                                                               aTTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGGCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGGCGTAAGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGAc                                      <  1:347917/149‑1 (MQ=255)                                                                                                                    gCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACaa                                 >  3:424485/1‑148 (MQ=255)                                                                                                                      cTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGATCGAGAGGCTGACTACGGTCGCCCACaaa                                >  2:356819/1‑147 (MQ=255)                                                                                                                         gAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGAt                           <  7:243050/149‑1 (MQ=255)                                                                                                                              gATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTc                                         <  4:325253/130‑1 (MQ=255)                                                                                                                              gATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTc                                         >  3:325253/1‑130 (MQ=255)                                                                                                                                tACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTg                    <  4:333697/149‑1 (MQ=255)                                                                                                                                  cccTAT‑CACGAGCCGATTGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGa                   <  7:205927/148‑1 (MQ=255)                                                                                                                                   ccTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACGTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTTATg                 >  6:42482/1‑149 (MQ=255)                                                                                                                                    cTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTGGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATg                 >  1:307735/1‑148 (MQ=255)                                                                                                                                       t‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTAc              <  4:422152/148‑1 (MQ=255)                                                                                                                                       t‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTAc              <  6:7109/148‑1 (MQ=255)                                                                                                                                       t‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTAc              <  2:343722/148‑1 (MQ=255)                                                                                                                                            gAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTAcatcat         >  8:240610/1‑149 (MQ=255)                                                                                                                                                 gATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCAtcgtcg    >  8:141780/1‑149 (MQ=255)                                                                                                                                                  aTAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGa                                                             <  4:324082/91‑1 (MQ=255)                                                                                                                                                  aTAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGa                                                             >  3:324082/1‑91 (MQ=255)                                                                                                                                                   tAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTcgcg  <  1:161189/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     CGCGGTCGTAGCCACGAAGAAGCAGGCAGCAAAGATCTGGCCACAACCAAGGGTGTGGAATCGAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTAT‑CACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGT‑CGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCG  >  NC_002947/951071‑951364

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