Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,533,852	C→T	6.4%	S621S (AGC→AGT)	secA →	protein translocase subunit SecA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,533,852	0	C	T	6.4%	187.1 / 2.3	62	S621S (AGC→AGT)	secA	protein translocase subunit SecA
Reads supporting (aligned to +/- strand):  ref base C (30/28);  new base T (2/2);  total (32/30)
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.35e-01

CAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCGGGCGAGGCCATCGAGCACCGCATGGTCACCAACGCTATCGAGAAAGCCCAGCGCAAGGTCGAAGGCCGCAACTTCGACATTCG  >  NC_002947/1533709‑1533998                                                                                                                                                |                                                                                                                                                   cAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGAt                                                                                                                                               >  8:1331/1‑149 (MQ=255)   gCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTTTACCTGTCGCTGGAAGACAGcc                                                                                                                                                   >  8:199066/1‑143 (MQ=255)   gCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTTTACCTGTCGCTGGAAGACAGcc                                                                                                                                                   <  7:199066/143‑1 (MQ=255)   gCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATgc                                                                                                                                             >  1:137233/1‑149 (MQ=255)      aGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCa                                                                                                                                          >  2:155606/1‑149 (MQ=255)             gagaCGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGc                                                                                                                                   <  1:155606/149‑1 (MQ=255)                       gCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGTCTGATGCGCATCTTCGCCTCTGACCgg                                                                                                                         >  3:125542/1‑149 (MQ=255)                       gCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGTCTGATGCGCATCTTCGCCTCTGACCgg                                                                                                                         >  3:120854/1‑149 (MQ=255)                       gCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCgg                                                                                                                         <  2:137233/149‑1 (MQ=255)                         cTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGt                                                                                                                       >  3:52958/1‑149 (MQ=255)                           gCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTg                                                                                                                      >  2:28776/1‑148 (MQ=255)                            cATGTAATCACTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATgc                                                                                                                                             <  1:274069/124‑1 (MQ=255)                            cATGTAATCACTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATgc                                                                                                                                             >  2:274069/1‑124 (MQ=255)                              tGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTgaag                                                                                                                   >  8:249676/1‑148 (MQ=255)                                tAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAAc                                                                                                                >  7:1339/1‑149 (MQ=255)                                  aTCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACtt                                                                                                              >  8:91390/1‑149 (MQ=255)                                     gCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCa                                                                                                            <  7:249676/148‑1 (MQ=255)                                     gCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCAt                                                                                                           >  1:273902/1‑149 (MQ=255)                                          cGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAg                                                                                                       <  8:1339/148‑1 (MQ=255)                                          cGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAg                                                                                                       <  5:82663/148‑1 (MQ=255)                                                 cACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGAt                                                                                                                                               <  7:88617/101‑1 (MQ=255)                                                 cACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGAt                                                                                                                                               >  8:88617/1‑101 (MQ=255)                                                 cACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTc                                                                                                                                     <  2:222682/111‑1 (MQ=255)                                                 cACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTc                                                                                                                                     >  1:222682/1‑111 (MQ=255)                                                 cACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTggg                                                                                               >  7:136874/1‑149 (MQ=255)                                                     aaTCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATg                                                                                           <  8:69517/149‑1 (MQ=255)                                                         gcgcCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGt                                                                                       >  7:17623/1‑149 (MQ=255)                                                            ccGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGTCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCgg                                                                                    <  4:120854/149‑1 (MQ=255)                                                            ccGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGTCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCgg                                                                                    <  4:125542/149‑1 (MQ=255)                                                               tATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCGGGc                                                                                  >  4:307308/1‑148 (MQ=255)                                                                  cGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACtt                                                                                                              <  3:144534/117‑1 (MQ=255)                                                                  cGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACtt                                                                                                              >  4:144534/1‑117 (MQ=255)                                                                  cGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCGGGCGAg                                                                               <  4:52958/148‑1 (MQ=255)                                                                   gACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCGGGCGAgg                                                                              <  1:28776/148‑1 (MQ=255)                                                                    aCAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCGGGCGAGGc                                                                             >  8:22344/1‑148 (MQ=255)                                                                    aCAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCGGGCGAGGc                                                                             >  7:85101/1‑148 (MQ=255)                                                                      aaCCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCGGGCGAGGCCAt                                                                          >  2:234054/1‑149 (MQ=255)                                                                      aaCCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCGGGCGAGGCCAt                                                                          <  7:22344/149‑1 (MQ=255)                                                                       aCCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCGGGCGAGGCCATc                                                                         <  8:136874/149‑1 (MQ=255)                                                                       aCCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCGGGCGAGGCCATc                                                                         <  8:85101/149‑1 (MQ=255)                                                                                     cGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCGGGCGAGGCCATCGAGCACCGCATGGt                                                           <  2:273902/149‑1 (MQ=255)                                                                                      gCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCGGGCGAGGCCATCGAGCACCGCATGGt                                                           >  2:188823/1‑148 (MQ=255)                                                                                      gCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCGGGCGAGGCCATCGAGCACCGCATGGt                                                           >  8:292473/1‑148 (MQ=255)                                                                                      gCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCGGGCGAGGCCATCGAGCACCGCATGGTc                                                          <  3:270102/149‑1 (MQ=255)                                                                                          cGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCGGGCGAGGCCATCGAGCACCGCATGGTCAc                                                        <  8:17623/147‑1 (MQ=255)                                                                                          cGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCGGGCGAGGCCATCGAGCACCGCATGGTCACCa                                                      <  1:234054/149‑1 (MQ=255)                                                                                            gtcgtcAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCGGGCGAGGCCATCGAGCACCGCATGGTCACCAAc                                                    >  5:106413/1‑149 (MQ=255)                                                                                             tcgtcAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCGGGCGAGGCCATCGAGCACCGCATGGTCACCAACg                                                   >  2:295922/1‑149 (MQ=255)                                                                                                  ctgacGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCGGGCGAGGCCATCGAGCACCGCATg                                                             >  3:293022/5‑134 (MQ=255)                                                                                                  ctgacGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCGGGCGAGGCCATCGAGCACCGCATg                                                             <  4:293022/130‑1 (MQ=255)                                                                                                  aGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCGGGCGAGGCCATCGAGCACCGCATGGTCACCAACGCTATc                                              <  1:295922/149‑1 (MQ=255)                                                                                                         cccGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCGGGCGAGGCCATCGAGCACCGCATGGTCACCAACGCTATCGAGAAAg                                       <  3:307308/149‑1 (MQ=255)                                                                                                           cGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCGGGCGAGGCCATCGAGCACCGCATGGTCACCAACGCTATCGAGAAAGcc                                     <  7:292473/149‑1 (MQ=255)                                                                                                               ttCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCGGGCGAGGCCATCGAGCACCGCATGGTCACCAACGCTATCGAGAAAGCCCAgc                                 <  7:91390/149‑1 (MQ=255)                                                                                                                            tACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACAGGGTGAAGAACTTCATGAAGTCACTGTGCATGCAGTCGGGCGATGACATCGAGCACAGCATGGTCACCAACGCGATCGACAAAGCACAGCGCAAGGTCGAAgg                    >  4:49866/1‑149 (MQ=255)                                                                                                                              ccTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCGGGCGAGGCCATCGAGCACCGCATGGTCACCAACGCTATCGAGAAAGCCCAg                                  >  2:84666/1‑133 (MQ=255)                                                                                                                              ccTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCGGGCGAGGCCATCGAGCACCGCATGGTCACCAACGCTATCGAGAAAGCCCAg                                  <  1:84666/133‑1 (MQ=255)                                                                                                                              ccTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCGGGCGAGGCCATCGAGCACCGCATGGTCACCAACGCTATCGAGAAAGCCCAGCGCAAGGTCGAAGGcc                  >  3:169191/1‑149 (MQ=255)                                                                                                                                      tGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCGGGCGAGGCCATCGAGCACCGCATGGTCACCAACGCTATCGAGAAAGCCCAGCGCAAGGTCGAAGGCCGCAACTTc          >  7:236036/1‑149 (MQ=255)                                                                                                                                        gAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCGGGCGAGGCCATCGAGCACCGCATGGTCACCAACGCTATCGAGAAAGCCCAgc                                 <  4:260438/124‑1 (MQ=255)                                                                                                                                        gAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCGGGCGAGGCCATCGAGCACCGCATGGTCACCAACGCTATCGAGAAAGCCCAgc                                 >  3:260438/1‑124 (MQ=255)                                                                                                                                             cAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCGGGCGAGGCCATCGAGCACCGCATGGTCACCAACGCTATCGAGAAAGCCCAGCGCAAGGTCGAAGGCCGCAACTTCGACATTc   <  4:169191/149‑1 (MQ=255)                                                                                                                                             cAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCGGGCGAGGCCATCGAGCACCGCATGGTCACCAACGCTATCGAGAAAGCCCAGCGCAAGGTCGAAGGCCGCAACTTCGACATTc   <  7:1331/149‑1 (MQ=255)                                                                                                                                              aGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCGGGCGAGGCCATCGAGCACCGCATGGTCACCAACGCTATCGAGAAAGCCCAGCGCAAGGTCGAAGGCCGCAACTTCGACATTCg  >  8:228549/1‑149 (MQ=255)                                                                                                                                                ccTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCGGGCGAGGCCATCGAGCACCGCATGGTCACCAACGCTATc                                              <  4:4035/103‑1 (MQ=255)                                                                                                                                                ccTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCGGGCGAGGCCATCGAGCACCGCATGGTCACCAACGCTATc                                              >  3:4035/1‑103 (MQ=255)                                                                                                                                                |                                                                                                                                                   CAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCGGGCGAGGCCATCGAGCACCGCATGGTCACCAACGCTATCGAGAAAGCCCAGCGCAAGGTCGAAGGCCGCAACTTCGACATTCG  >  NC_002947/1533709‑1533998

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