Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,303,894	A→T	6.3%	D177V (GAC→GTC)	cysS →	cysteine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,303,894	0	A	T	6.3%	189.0 / 2.2	63	D177V (GAC→GTC)	cysS	cysteine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (30/29);  new base T (2/2);  total (32/31)
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.66e-01

GATCGACAAAGGGTACGCCTACGCGCCAGGCAATGGCGACGTGTACTACCGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCAC  >  NC_002947/3303748‑3304036                                                                                                                                                   |                                                                                                                                               gATCGACAAAGGGTACGCCTACGCGCCAGGCAATGGCGACGTGTACTACCGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGAc                                                                                                                                               <  6:167235/148‑1 (MQ=255)             gTACGCCTACGCGCCAGGCAATGGCGACGTGTACTACCGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCg                                                                                                                                  >  7:239505/1‑149 (MQ=255)               aCGCCTACGCGCCAGGCAATGGCGACGTGTACTACCGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTc                                                                                                                                >  2:215785/1‑149 (MQ=255)                 gCCTACGCGCCAGGCAATGGCGACGTGTACTACCGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCt                                                                                                                              >  5:145249/1‑149 (MQ=255)                   cTACGCGCCAGGCAATGGCGACGTGAACTACCGGGTTGTCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCtgt                                                                                                                            <  1:237044/149‑1 (MQ=255)                      agcgcCAGGCAATTGCGACGTGTACTACAGGGGTGTCAAGTGCGTCGGCTACGGCAAGCTGTCGAGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTAGACTTCGTCCTGTGGa                                                                                                                         <  6:234294/148‑1 (MQ=255)                         gcCAGGCAATGGCGACGTGTACTACCGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAgg                                                                                                                      >  4:15817/1‑149 (MQ=255)                            aGGCAATAGCGACGTGTACTACCGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCg                                                                                                                   <  1:19375/149‑1 (MQ=255)                                   ggCGACGTGTACTACCGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGcc                                                                                                            >  8:47009/1‑149 (MQ=255)                                            tactacCGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGcc                                                                                                   <  7:47009/149‑1 (MQ=255)                                                 ccGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGc                                                                                                    <  4:10035/143‑1 (MQ=255)                                                 ccGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGc                                                                                                    >  3:10035/1‑143 (MQ=255)                                                   gggTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTgg                                                                                            <  5:369431/149‑1 (MQ=255)                                                            aGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCa                                                                                   <  6:145249/149‑1 (MQ=255)                                                                  tcgGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGc                                                                             >  6:357894/1‑149 (MQ=255)                                                                    ggCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCt                                                                                              <  3:159842/130‑1 (MQ=255)                                                                    ggCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCt                                                                                              >  4:159842/1‑130 (MQ=255)                                                                     gCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCg                                                                          >  3:174642/1‑149 (MQ=255)                                                                         cGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTc                                                                      <  6:35924/149‑1 (MQ=255)                                                                          ggCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCgggcga                                                                     >  7:64275/1‑146 (MQ=255)                                                                          ggCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCgg                                                                         <  6:210205/145‑1 (MQ=255)                                                                          ggCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCgg                                                                         >  5:210205/1‑145 (MQ=255)                                                                               gCTGTCGCGCAAGAGGAGCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCgg                                                                <  1:60248/149‑1 (MQ=255)                                                                                   tCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCt                                                             <  3:259588/148‑1 (MQ=255)                                                                                   tCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCt                                                             >  4:259588/1‑148 (MQ=255)                                                                                     gcgcAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTgg                                                           <  8:239505/148‑1 (MQ=255)                                                                                      cgcAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGca                                                         <  3:15817/149‑1 (MQ=255)                                                                                      cgcAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGCGACCCGGGccccat                                                         <  8:64275/149‑6 (MQ=255)                                                                                         aaGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCt                                                                                              <  7:33272/109‑1 (MQ=255)                                                                                         aaGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCt                                                                                              >  8:33272/1‑109 (MQ=255)                                                                                          agagGATCGAAGACCTGCGTATCGGGGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATc                                                     <  7:116480/149‑1 (MQ=255)                                                                                            agGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATc                                                     >  2:275200/1‑147 (MQ=255)                                                                                             ggATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGa                                                   <  3:24043/148‑1 (MQ=255)                                                                                              gATCGAAGACCTGCGTATCGGAGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGt                                                 >  2:257015/1‑149 (MQ=255)                                                                                               aTCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTg                                                >  7:238953/1‑149 (MQ=255)                                                                                                tCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGc                                               >  2:72346/1‑149 (MQ=255)                                                                                                      aCCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTc                                         <  1:215785/149‑1 (MQ=255)                                                                                                        cTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCAt                                       <  3:264726/149‑1 (MQ=255)                                                                                                          gCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGt                                     >  5:205770/1‑149 (MQ=255)                                                                                                              aTCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCAc                                 >  3:137049/1‑149 (MQ=255)                                                                                                              aTCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCAc                                 >  4:187190/1‑149 (MQ=255)                                                                                                                  gTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCCCATCGAGTGCTCGGTCATGTcc                                   >  7:242402/1‑143 (MQ=255)                                                                                                                  gTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTcc                                   <  8:242402/143‑1 (MQ=255)                                                                                                                     cACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCtgctgc                          <  1:72346/149‑1 (MQ=255)                                                                                                                        gCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCg                                                                 >  5:23787/1‑107 (MQ=255)                                                                                                                        gCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCg                                                                 <  6:23787/107‑1 (MQ=255)                                                                                                                             gAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGAGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGc                    >  8:193909/1‑147 (MQ=255)                                                                                                                                caggacgaAGTCGAGCGGGTCCTGCTTGGCTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCAc                                 >  7:42791/4‑131 (MQ=38)                                                                                                                                caggacgaAGTCGAGCGGGTCCTGCTTGGCTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCAc                                 <  8:42791/128‑1 (MQ=38)                                                                                                                                 tCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCgag                 <  7:119594/146‑1 (MQ=255)                                                                                                                                 tCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCgag                 >  8:119594/1‑146 (MQ=255)                                                                                                                                 aggacgaAGTCGAGCGGGTCCTGCTTGGCTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGc              >  7:118388/3‑149 (MQ=38)                                                                                                                                   gacgaAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCtt            >  7:31188/1‑149 (MQ=255)                                                                                                                                   gacgaAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCtt            >  8:306875/1‑149 (MQ=255)                                                                                                                                   gacgaAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCtt            <  3:187190/149‑1 (MQ=255)                                                                                                                                   gacgaAGCCAAGCAGGACCCACTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTgg                                                                                            <  7:242814/69‑1 (MQ=255)                                                                                                                                   gacgaAGCCAAGCAGGACCCACTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTgg                                                                                            >  8:242814/1‑69 (MQ=255)                                                                                                                                    acgaAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTc           >  3:130308/1‑149 (MQ=255)                                                                                                                                            aaGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGc                    <  1:100368/132‑1 (MQ=255)                                                                                                                                            aaGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGc                    >  2:100368/1‑132 (MQ=255)                                                                                                                                             aGCGGGTCCTGCTTGGCTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCAc  <  8:118388/149‑1 (MQ=255)                                                                                                                                              gCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACctgct                            >  2:155557/1‑122 (MQ=255)                                                                                                                                              gCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACctgct                            <  1:155557/122‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                               GATCGACAAAGGGTACGCCTACGCGCCAGGCAATGGCGACGTGTACTACCGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCAC  >  NC_002947/3303748‑3304036

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