Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,937,048:1	+C	8.8%	coding (1465/1776 nt)	fucI →	L‑fucose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,937,048	1	.	C	8.8%	346.5 / 30.8	91	coding (1465/1776 nt)	fucI	L‑fucose isomerase
Reads supporting (aligned to +/- strand):  ref base . (38/45);  new base C (3/5);  total (41/50)
Fisher's exact test for biased strand distribution p-value = 7.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGC  >  NC_000913/2936902‑2937196                                                                                                                                                    |                                                                                                                                                     aCTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACa                                                                                                                                                   >  2:358847/1‑150 (MQ=255)    tCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATcc                                                                                                                                                >  2:438663/1‑150 (MQ=255)        tGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTcaa                                                                                                                                            >  1:725347/1‑150 (MQ=255)        tGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGT‑CATGCACATCCttcg                                                                                                                                            <  2:374707/150‑4 (MQ=255)         ggcggcGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTcaac                                                                                                                                           <  2:1243949/150‑1 (MQ=255)         ggcggcGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGT‑CATGCACATCCttcgg                                                                                                                                           <  2:716098/150‑5 (MQ=255)          gcggcgGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTTTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGT‑CATGCACATCCttcggc                                                                                                                                          <  1:480399/150‑6 (MQ=255)          gcggcgGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGT‑CATGCACATCCttcggc                                                                                                                                          <  1:475062/150‑6 (MQ=255)           cggcggTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTcaacaa                                                                                                                                         <  2:394435/150‑1 (MQ=255)                gTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCa                                                                                                                                    <  2:880645/150‑1 (MQ=255)                   aCTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCa                                                                                                                                 >  1:1194396/1‑150 (MQ=255)                    ctctTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCaa                                                                                                                                >  1:701552/1‑150 (MQ=255)                         cccGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCaa                                                                                                                           >  1:1296236/1‑150 (MQ=255)                           cGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAAcc                                                                                                                         <  1:358847/150‑1 (MQ=255)                               tCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGc                                                                                                                     <  2:1250807/150‑1 (MQ=255)                                ccTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGcc                                                                                                                    <  2:338412/150‑1 (MQ=255)                                ccTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGcc                                                                                                                    >  1:308569/1‑150 (MQ=255)                                ccTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGcc                                                                                                                    >  1:265233/1‑150 (MQ=255)                                    aCCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAacc                                                                                                                <  1:1470292/150‑1 (MQ=255)                                     ccGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTTCTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAacca                                                                                                               >  2:444371/1‑150 (MQ=255)                                     ccGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAacca                                                                                                               <  1:573762/150‑1 (MQ=255)                                       gAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAaccacc                                                                                                             >  2:138495/1‑150 (MQ=255)                                       gAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAaccacc                                                                                                             <  2:945064/150‑1 (MQ=255)                                        aaGGCGGCGTCCCCTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGAGGAATTGAAGAAGGATGTGCATG‑ACATCATCAACAAACGCACCAACTCAACCTGGCCAAACAACt                                                                                                            >  2:287441/1‑150 (MQ=255)                                          ggcggcGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTgg                                                                                                          <  1:1012269/150‑1 (MQ=255)                                            cggcgTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTcaac                                                                                                                                           <  1:482011/115‑1 (MQ=255)                                            cggcgTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGtt                                                                                                        <  1:636524/150‑1 (MQ=255)                                            cggcgTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATAGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTcaac                                                                                                                                           >  2:482011/1‑115 (MQ=255)                                             ggcggcCCGTTCACCATGACTCGTGGCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGttt                                                                                                       <  2:1351247/150‑1 (MQ=255)                                                  cccGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAAc                                                                                                                                       <  1:163284/113‑1 (MQ=255)                                                  cccGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAAc                                                                                                                                       >  2:163284/1‑113 (MQ=255)                                                    cGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACcgc                                                                                                >  2:695633/1‑150 (MQ=255)                                                      ttCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCgcgc                                                                                              <  2:943522/150‑1 (MQ=255)                                                       tCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCCCCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGcc                                                                                             <  1:312556/150‑1 (MQ=255)                                                        cACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCt                                                                                            <  1:1002339/150‑1 (MQ=255)                                                         aCCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTc                                                                                           >  2:1019806/1‑150 (MQ=255)                                                         aCCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTc                                                                                           <  2:98205/150‑1 (MQ=255)                                                          ccATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCa                                                                                          >  2:919381/1‑150 (MQ=255)                                                              gACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCgg                                                                                      <  2:701552/150‑1 (MQ=255)                                                              gACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCgg                                                                                      <  2:866315/150‑1 (MQ=255)                                                                cTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTa                                                                                    <  1:787823/150‑1 (MQ=255)                                                                cTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTa                                                                                    <  2:283124/150‑1 (MQ=255)                                                                  cGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTaaa                                                                                  >  1:1353072/1‑150 (MQ=255)                                                                  cGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTaaa                                                                                  >  2:1116330/1‑150 (MQ=255)                                                                   gtgtCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAg                                                                                 <  1:695633/150‑1 (MQ=255)                                                                    tgtCAACATCATCAAAGGCCTGGGACCGGTACTGGAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGTATGTGCATG‑AGATGTTCAATAAACGCATCATCTAAAACTGTCCAAGCACCTGTTTTTCACTGATCCTCTCCTGTAATgg                                                                                <  2:572491/150‑1 (MQ=255)                                                                     gtCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAggg                                                                               <  2:534851/150‑1 (MQ=255)                                                                     gtCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAggg                                                                               <  2:344138/150‑1 (MQ=255)                                                                     gtCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAggg                                                                               <  1:763829/150‑1 (MQ=255)                                                                      tCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGc                                                                              <  2:1020661/150‑1 (MQ=255)                                                                       cAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGcc                                                                             <  1:723697/150‑1 (MQ=255)                                                                        aaCATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGc                                                                              <  1:663610/148‑1 (MQ=255)                                                                        aaCATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGc                                                                              >  2:663610/1‑148 (MQ=255)                                                                         aCATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGt                                                                           <  1:1301555/150‑1 (MQ=255)                                                                          catcatCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACGAACTCAACCTGGCCATCCACCTGGTTTGGACCGCGCCTCACCGGTAAAGGGCCGtt                                                                          <  2:847293/150‑1 (MQ=255)                                                                          catcatCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGtt                                                                          >  1:1054326/1‑150 (MQ=255)                                                                             catcaAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTAc                                                                       >  2:168584/1‑150 (MQ=255)                                                                               tcaAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACgg                                                                     >  1:346132/1‑150 (MQ=255)                                                                                 aaaGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGAt                                                                   >  1:624839/1‑150 (MQ=255)                                                                                    ggCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGAtgtg                                                                >  2:78531/1‑150 (MQ=255)                                                                                    ggCCTGGGACCGGTACTGAAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTCGCACCGCGCCTCACCGGTAAATGGCCGTTTACGGAtgtg                                                                <  2:588590/150‑1 (MQ=255)                                                                                           gACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCgg                                                         <  2:1194396/150‑1 (MQ=255)                                                                                            aCCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGt                                                        <  1:168584/150‑1 (MQ=255)                                                                                               ggTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAAt                                                     >  2:472687/1‑150 (MQ=255)                                                                                                 tACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATgg                                                   <  1:919381/150‑1 (MQ=255)                                                                                                 tACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATgg                                                   <  2:448929/150‑1 (MQ=255)                                                                                                  aCTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGc                                                  <  1:28599/150‑1 (MQ=255)                                                                                                     gCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAACCGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGaa                                               >  1:1452325/1‑150 (MQ=255)                                                                                                          tCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTggg                                          >  2:882051/1‑150 (MQ=255)                                                                                                              ggAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGgcgc                                      <  1:43530/150‑1 (MQ=255)                                                                                                               gAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCt                                     <  2:624839/150‑1 (MQ=255)                                                                                                                 aGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTaa                                   >  1:1001153/1‑150 (MQ=255)                                                                                                                  ggCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCAGCGCGCCTGACCGGTAAAGGGCCGTTTTCGGATGTGTACTCGGTAATGGCGAACTTGGGCGCTAAc                                  <  2:725347/150‑1 (MQ=255)                                                                                                                   gCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAAcc                                 >  1:209468/1‑150 (MQ=255)                                                                                                                       gAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATgg                             <  2:308569/150‑1 (MQ=255)                                                                                                                         gCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACg                                                                      <  1:666362/107‑1 (MQ=255)                                                                                                                         gCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACg                                                                      >  2:666362/1‑107 (MQ=255)                                                                                                                         gCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATgggg                           <  1:381090/150‑1 (MQ=255)                                                                                                                            tGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTc                        >  1:1269574/1‑150 (MQ=255)                                                                                                                                aTTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGAc                    >  2:1340978/1‑150 (MQ=255)                                                                                                                                    gttgaGGATGT‑CATGCACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATc                >  1:850185/5‑150 (MQ=255)                                                                                                                                    ccGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATc                >  1:1372340/1‑150 (MQ=255)                                                                                                                                     ttgaGGATGT‑CATGCACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCg               >  1:449661/4‑150 (MQ=255)                                                                                                                                      tgaGGATGT‑CATGCACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCgg              >  1:851546/3‑150 (MQ=255)                                                                                                                                      gAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCgg              <  2:1225145/150‑1 (MQ=255)                                                                                                                                       gaGGATGT‑CATGCACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGc             <  2:449661/149‑1 (MQ=255)                                                                                                                                          gATGTGCATG‑ACATCCTCAACAAACGTACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCAc          >  1:432207/1‑150 (MQ=255)                                                                                                                                               gCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTgg     >  1:1215848/1‑150 (MQ=255)                                                                                                                                                cATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGgc    >  1:193552/1‑150 (MQ=255)                                                                                                                                                 aTG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCCCCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCCCGTTGgcg   >  2:312668/1‑150 (MQ=255)                                                                                                                                                  tG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGgcgc  >  2:54923/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                     ACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATG‑ACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGC  >  NC_000913/2936902‑2937196

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