Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,935,557	C→G	7.8%	S193W (TCG→TGG)	fau →	5‑formyltetrahydrofolate cyclo‑ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,935,557	0	C	G	7.8%	197.6 / 4.4	62	S193W (TCG→TGG)	fau	5‑formyltetrahydrofolate cyclo‑ligase
Reads supporting (aligned to +/- strand):  ref base C (30/27);  new base G (3/2);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.47e-01

GGCGGTGGCTTCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGACCA  >  NC_002947/5935409‑5935699                                                                                                                                                     |                                                                                                                                                ggCGGTGGCTTCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGtctc                                                                                                                                                 <  3:286716/149‑1 (MQ=255) ggCGGTGGCTTCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGtctc                                                                                                                                                 <  3:122987/149‑1 (MQ=255) ggCGGTGGCTTCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGtctc                                                                                                                                                 <  1:51858/149‑1 (MQ=255)    ggTGGCTTCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGa                                                                                                                                              <  4:297369/149‑1 (MQ=255)          ttCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGGATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTgg                                                                                                                                        >  4:107167/1‑149 (MQ=255)             tATGAGCGCGGCCTGGCCTACCCGGCGCGTCGGCACACCTGGCAAAGACGCTTGCTGCTGGGGCTGGCGCATGAGGGCGGGCAGGTCGGGGGGCGGGCGCAGGGCGGGGGGGAGGTTGCGCGGCAGGGCACGGTGTGGGACGGTGGcgg                                                                                                                                     >  8:63632/1‑147 (MQ=255)              aTGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGTCGCAGTCCAGTTGGTATGTTCCGCTGCAGTCCACGGTCCGAGACCGTgccct                                                                                                                                     <  4:217993/148‑5 (MQ=255)                 aTCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTAc                                                                                                                                 <  8:24488/149‑1 (MQ=255)                        ccTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGgcgc                                                                                                                          >  7:16054/1‑149 (MQ=255)                              cTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGtt                                                                                                                    >  6:189076/1‑149 (MQ=255)                              cTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGtt                                                                                                                    >  4:91897/1‑149 (MQ=255)                                aCCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGc                                                                                                                                       >  7:116095/1‑128 (MQ=255)                                aCCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGc                                                                                                                                       <  8:116095/128‑1 (MQ=255)                                    ggCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAgcgg                                                                                                              >  6:12144/1‑149 (MQ=255)                                       gcgTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGgcgc                                                                                                           <  4:200493/149‑1 (MQ=255)                                            ggCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCAc                                                                                                                       <  6:294422/132‑1 (MQ=255)                                            ggCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCAc                                                                                                                       >  5:294422/1‑132 (MQ=255)                                              cAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCggg                                                                                                    <  4:48953/149‑1 (MQ=255)                                               aaaCCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCgggg                                                                                                   <  6:14511/149‑1 (MQ=255)                                                aaCCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCggggg                                                                                                  <  5:189076/149‑1 (MQ=255)                                                     ggAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑Tc                                                                                            <  5:12144/149‑1 (MQ=255)                                                      gAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑Tc                                                                                            >  7:132528/1‑148 (MQ=255)                                                          aaaCCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGt                                                                                       >  3:70863/1‑149 (MQ=255)                                                            aCCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTtg                                                                                     >  1:16583/1‑149 (MQ=255)                                                              ccTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTtgtg                                                                                   >  3:198066/1‑149 (MQ=255)                                                                   ctgctgGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGtggctg                                                                               <  2:16583/148‑1 (MQ=255)                                                                   cggctgGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGtggctg                                                                               <  4:59276/146‑1 (MQ=255)                                                                           gCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGtggctggctggc                                                                         <  8:33821/146‑1 (MQ=255)                                                                           gCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGtggctggctggc                                                                         >  7:33821/1‑146 (MQ=255)                                                                                 gcATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGTGGCTGGCTGGCGTCGACCgg                                                                >  7:76332/1‑149 (MQ=255)                                                                                 gcATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGTGGCTGGCTGGCGTCGACCgg                                                                >  7:202152/1‑149 (MQ=255)                                                                                          ccAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGa                                                         >  6:66583/1‑147 (MQ=255)                                                                                            aGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAgg                                                     <  8:16054/149‑1 (MQ=255)                                                                                             gAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGc                                                    >  4:128206/1‑149 (MQ=255)                                                                                                    gAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTc                                              >  8:223148/1‑148 (MQ=255)                                                                                                         gCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGtt                                                 <  7:166406/140‑1 (MQ=255)                                                                                                         gCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGtt                                                 <  3:90222/140‑1 (MQ=255)                                                                                                         gCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGtt                                                 >  8:166406/1‑140 (MQ=255)                                                                                                         gCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGtt                                                 >  4:90222/1‑140 (MQ=255)                                                                                                         gCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCtt                                        >  8:59206/1‑149 (MQ=255)                                                                                                             gcgcAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTgcgc                                    >  1:269587/1‑149 (MQ=255)                                                                                                               gcAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTc                                  >  2:186146/1‑149 (MQ=255)                                                                                                               gcAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGTGGCTGGCTGGCGTCGACCGGCATATGATAGGCGTTGTCCAGCTTGCGCTc                                  <  3:91897/149‑1 (MQ=255)                                                                                                                cAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTcc                                 <  5:66583/149‑1 (MQ=255)                                                                                                                 aGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTccc                                <  8:132528/149‑1 (MQ=255)                                                                                                                 aGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTccc                                <  2:269587/149‑1 (MQ=255)                                                                                                                   gCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAg                                                      >  8:169811/1‑125 (MQ=255)                                                                                                                   gCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAg                                                      <  7:169811/125‑1 (MQ=255)                                                                                                                   gCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTcc                                 >  2:207375/1‑146 (MQ=255)                                                                                                                          gggATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTgcg                                     >  1:304640/1‑135 (MQ=255)                                                                                                                          gggATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTgcg                                     <  2:304640/135‑1 (MQ=255)                                                                                                                           ggATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGtt                                                                                      <  5:314711/85‑1 (MQ=255)                                                                                                                           ggATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGtt                                                                                      >  6:314711/1‑85 (MQ=255)                                                                                                                           ggATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTg                      >  5:12022/1‑149 (MQ=255)                                                                                                                           ggATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTg                      >  3:316530/1‑149 (MQ=255)                                                                                                                            gATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTgc                     <  3:128206/149‑1 (MQ=255)                                                                                                                                        cAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTa         <  7:209738/149‑1 (MQ=255)                                                                                                                                          gggCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTAGGAAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAg        >  1:103886/1‑149 (MQ=255)                                                                                                                                          gggCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGt       >  6:70651/1‑149 (MQ=255)                                                                                                                                           ggCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTg      <  7:59206/149‑1 (MQ=255)                                                                                                                                           gccACGGTCCGAGACCGTGGCTGGTCCCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGCTAGGCGTTGTCCCGCTTGCGCTCCCAGACGCCTTGCGCATAACCGGTAggg      >  6:344696/3‑147 (MQ=255)                                                                                                                                            gCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGa     <  1:207375/149‑1 (MQ=255)                                                                                                                                            ccACGGTCCGAGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGACCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTa         >  1:302282/2‑145 (MQ=255)                                                                                                                                            ccACGGTCCGAGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGACCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTa         <  2:302282/144‑1 (MQ=255)                                                                                                                                               cGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTGGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGACca  <  4:316530/149‑1 (MQ=255)                                                                                                                                                  tctcGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGAc    <  8:101968/144‑1 (MQ=255)                                                                                                                                                  tctcGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGAc    >  7:101968/1‑144 (MQ=255)                                                                                                                                                     |                                                                                                                                                GGCGGTGGCTTCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTA‑TCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGACCA  >  NC_002947/5935409‑5935699

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