Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,199,963:1	+C	7.7%	coding (829/1278 nt)	PP_1944 →	aminomethyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,199,963	1	.	C	7.7%	248.6 / 14.0	65	coding (829/1278 nt)	PP_1944	aminomethyltransferase
Reads supporting (aligned to +/- strand):  ref base . (34/26);  new base C (3/2);  total (37/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.01e-01

ACTACATGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGCCTCGCTGTT  >  NC_002947/2199823‑2200111                                                                                                                                              |                                                                                                                                                     aCTACATGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAgctg                                                                                                                                              >  6:380032/1‑149 (MQ=255)      aTGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGCTCGa‑ctgggctg                                                                                                                                         <  5:200018/149‑1 (MQ=255)                  ttCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCa                                                                                                                              <  2:187566/148‑1 (MQ=255)                      aCCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCAGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAg                                                                                                                         >  1:66412/1‑149 (MQ=255)                          agagATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTc                                                                                                                      >  6:139710/1‑148 (MQ=255)                           gagaTGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGa                                                                                                                    <  5:380032/149‑1 (MQ=255)                                   gggTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACt                                                                                                            >  5:150439/1‑149 (MQ=255)                                           tcAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGt                                                                                                    >  1:166364/1‑149 (MQ=255)                                           tcAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGt                                                                                                    >  4:316505/1‑149 (MQ=255)                                            cAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTc                                                                                                   >  1:179950/1‑149 (MQ=255)                                            cAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGCTATCGCAGCCCAG‑TCGAGCAGGGCAGGGCCAAAACCAACAAGTTCGATCACGACTTAATCGGTc                                                                                                   >  1:174242/1‑149 (MQ=255)                                             aaGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAgctggg                                                                                                                                            <  1:199990/107‑1 (MQ=255)                                             aaGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAgctggg                                                                                                                                            >  2:199990/1‑107 (MQ=255)                                                 aaTACCGCGCGGGTATGCCATGTATTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGAATGGTATCGCAGCCCAG‑ACGACCTGGGCTGGGCCAAAACCAACTACTTCGATCACGACTATATCGGTCGttgc                                                                                              >  2:339630/1‑146 (MQ=255)                                                      cgcgcgGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGt                                                                                          <  2:66412/148‑1 (MQ=255)                                                      cgcgcgGGTATGCCAAGTATTGCCTCGACTTTCAATGTTGCGGCGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCAACAAGTTCGATCACGACATTATCGGTCGTGCCGCGtt                                                                                         >  8:173413/1‑149 (MQ=255)                                                                 gCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGc                                                                                               >  6:28053/1‑132 (MQ=255)                                                                 gCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGc                                                                                               <  5:28053/132‑1 (MQ=255)                                                                 gCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAgg                                                                              >  2:74498/1‑149 (MQ=255)                                                                   cAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGtt                                                                            >  6:334912/1‑149 (MQ=255)                                                                    aatttttGCCTCGTCTTTCAATGTTGGGGGGATCTTTGTGTCGAATAATGTATGCGATTGGTATCGCAGCCCAG‑TGGAGGTGGGGTGGGGCAAAACCATCTAGTTCGATCACGACTTTATCGGTCGTGGCGCGTTGGAAAAAGAGGTTg                                                                           <  2:174242/146‑1 (MQ=255)                                                                    aaGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTg                                                                           <  2:179950/149‑1 (MQ=255)                                                                         ttGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCaga                                                                       <  5:184963/148‑1 (MQ=255)                                                                             cTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGt                                                                             <  3:325673/138‑1 (MQ=255)                                                                             cTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGt                                                                             >  4:325673/1‑138 (MQ=255)                                                                             cTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTc                                                                  >  6:101564/1‑149 (MQ=255)                                                                                   tttCAATGTTGCGGGGAGCTTTGAGTCTATTAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCTATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACggg                                                            <  1:23460/149‑1 (MQ=255)                                                                                   tttCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACggg                                                            >  4:51197/1‑149 (MQ=255)                                                                                     taaaTGTTGCCGGGACCTTAGACTCGAATAACGTAAGCGATTGGTATCGCAGCCCAC‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTACGATCACGACTTTATAGGTCGAGCCGCGATGGAAAAAGAGGTTGGAGAGACTAGAcagct                                                           >  4:359659/3‑144 (MQ=255)                                                                                     tCAATGTTGCGGGGAGCTTTGAGTCGAATAAGTTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCGGGTCTTGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGGCGTTCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACggg                                                            <  8:297503/147‑1 (MQ=255)                                                                                          gTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATgcg                                                     >  7:234391/1‑149 (MQ=255)                                                                                                gggAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCtt                                               >  4:368249/1‑149 (MQ=255)                                                                                                gggAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCtt                                               >  3:123121/1‑149 (MQ=255)                                                                                                 ggAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTg                                              >  2:199564/1‑149 (MQ=255)                                                                                                 ggAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTg                                              <  1:199564/149‑1 (MQ=255)                                                                                                   aGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGa                                                                                                                    <  5:113230/77‑1 (MQ=255)                                                                                                   aGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGa                                                                                                                    >  6:113230/1‑77 (MQ=255)                                                                                                     cTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCCCGACGGGTGATGCGCACCTTGGttt                                          >  2:111908/1‑149 (MQ=255)                                                                                                      tttGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTg                                         <  3:368249/149‑1 (MQ=255)                                                                                                        tGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTg                                         <  1:161267/147‑1 (MQ=255)                                                                                                        tGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTg                                         <  3:51197/147‑1 (MQ=255)                                                                                                        ggAGTCGAATTCGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCTCGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTg                                         <  1:339630/146‑1 (MQ=255)                                                                                                              gAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTg                                         <  5:9107/141‑1 (MQ=255)                                                                                                              gAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTg                                         >  6:9107/1‑141 (MQ=255)                                                                                                              gAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGtg                                 <  5:119331/149‑1 (MQ=255)                                                                                                              gAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGt                                  <  5:334912/148‑1 (MQ=255)                                                                                                               aaTAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACggg                                                            >  3:223774/1‑121 (MQ=255)                                                                                                               aaTAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACggg                                                            <  4:223774/121‑1 (MQ=255)                                                                                                               aaTAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGtga                                >  2:261159/1‑149 (MQ=255)                                                                                                                aTAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGtgat                               >  8:235283/1‑149 (MQ=255)                                                                                                                aTAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGtga                                <  3:67556/148‑1 (MQ=255)                                                                                                                  aaGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGtgatga                             >  3:366183/1‑149 (MQ=255)                                                                                                                    ggTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGaaaaa                                                                                  <  8:200660/94‑1 (MQ=255)                                                                                                                    ggTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGaaaaa                                                                                  >  7:200660/1‑94 (MQ=255)                                                                                                                       aaGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTg                        >  5:222099/1‑149 (MQ=255)                                                                                                                          cGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTg                        >  6:376491/1‑146 (MQ=255)                                                                                                                                gTATGGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGtt                <  3:123356/148‑1 (MQ=255)                                                                                                                                gTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTc               <  6:222099/149‑1 (MQ=255)                                                                                                                                  ggcccAGCCCAGCTCGA‑CTGGGCTGCGCCAAAACCATCAAGTTCGATCACGACTTTTTCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTGGGTTTGGAATGGTGATGACGTTGCAGATGTTCAt             >  1:228249/5‑149 (MQ=255)                                                                                                                                   gcccAGCCCAGCTCGA‑CTGGGCTGCGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCa                      <  7:159369/136‑1 (MQ=255)                                                                                                                                   gcccAGCCCAGCTCGA‑CTGGGCTGCGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCa                      >  8:159369/4‑139 (MQ=255)                                                                                                                                    cGCAGCCCAG‑TCGAGCCGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGc           >  8:281095/1‑149 (MQ=255)                                                                                                                                     gCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGt                                                           >  4:130021/1‑100 (MQ=255)                                                                                                                                     gCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGt                                                           <  3:130021/100‑1 (MQ=255)                                                                                                                                     ccAGCCCAGCTCGA‑CTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGcc          >  1:333666/2‑149 (MQ=255)                                                                                                                                             g‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGCCTCGCTGtt  <  3:316505/149‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                     ACTACATGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGCCTCGCTGTT  >  NC_002947/2199823‑2200111

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