Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,348,789	C→A	5.5%	intergenic (+21/+33)	yfaE → / ← inaA	ferredoxin involved with ribonucleotide reductase diferric‑tyrosyl radical (Y*) cofactor maintenance/acid‑inducible Kdo/WaaP family putative kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,348,789	0	C	A	5.5%	250.1 / 2.6	72	intergenic (+21/+33)	yfaE/inaA	ferredoxin involved with ribonucleotide reductase diferric‑tyrosyl radical (Y*) cofactor maintenance/acid‑inducible Kdo/WaaP family putative kinase
Reads supporting (aligned to +/- strand):  ref base C (32/36);  new base A (2/2);  total (34/38)
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.60e-01

GCACACGCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATCGCGACGC  >  NC_000913/2348650‑2348930                                                                                                                                            |                                                                                                                                               gCACACGCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTa                                                                                                                                        <  3:339324/148‑1 (MQ=255)      cGCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTc                                                                                                                                  <  7:81333/149‑1 (MQ=255)       gCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCt                                                                                                                                 <  5:257959/149‑1 (MQ=255)         cTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTa                                                                                                                               >  4:198107/1‑149 (MQ=255)            gTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATc                                                                                                                             <  3:152229/148‑1 (MQ=255)             ttGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTg                                                                                                                           >  5:204766/1‑149 (MQ=255)               gCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGc                                                                                                                         >  2:140402/1‑149 (MQ=255)                cAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGc‑c                                                                                                                       <  6:129008/149‑1 (MQ=255)                cAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGc                                                                                                                         <  3:303761/148‑1 (MQ=255)                cAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGc                                                                                                                         <  8:113152/148‑1 (MQ=255)                cAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATAAGACGCATAAGCGTTGTATCTGGca                                                                                                                        <  6:279444/149‑2 (MQ=39)                 aGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGc‑c                                                                                                                       >  7:177992/1‑148 (MQ=255)                   gTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTa                                                                                                                     >  7:153868/1‑148 (MQ=255)                    tCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTAc                                                                                                                    >  1:303669/1‑148 (MQ=255)                      aaGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAAc                                                                                                                 >  4:258303/1‑149 (MQ=255)                             cTGGATTGCCGAACCGTTAGCCTTTATTAAGCCGGGGGAAATTTTGCCCAGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACa                                                                                                          >  2:315340/1‑149 (MQ=255)                              tGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACAt                                                                                                         >  8:97375/1‑149 (MQ=255)                                gATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACAtcg                                                                                                       >  4:350701/1‑147 (MQ=255)                                gATTGCCGAACCGTTAGCCTTTATTCAGCCGCGGTAAATTTTGCCCTGTTGTGGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACAt                                                                                                         <  3:191521/147‑1 (MQ=255)                                     ccGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCgtag                                                                                                  <  8:153868/149‑1 (MQ=255)                                      cGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCgtagt                                                                                                 <  5:219329/149‑1 (MQ=255)                                        aaCCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCGTAGTaa                                                                                               >  5:238281/1‑149 (MQ=255)                                          ccGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCACTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTTGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCGTAGTAATa                                                                                             <  4:202190/149‑1 (MQ=255)                                           cGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTCCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCGTAGTAATAc                                                                                            <  5:337568/149‑1 (MQ=255)                                           cGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGGGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCGTAGTAATAc                                                                                            <  7:336726/149‑1 (MQ=255)                                             ttAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGc                                                                                                                                           <  2:8633/101‑1 (MQ=255)                                             ttAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGc                                                                                                                                           >  1:8633/1‑101 (MQ=255)                                             ttAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCGTAGTAATACGc                                                                                          <  8:177992/149‑1 (MQ=255)                                                  cTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATgg                                                                                                       >  3:278150/1‑131 (MQ=255)                                                  cTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATgg                                                                                                       <  4:278150/131‑1 (MQ=255)                                                  cTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCGTAGTAATACGCTTTCa                                                                                     >  8:17725/1‑149 (MQ=255)                                                   tttATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGAtt                                                                                                             <  4:221025/124‑1 (MQ=255)                                                   tttATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGAtt                                                                                                             >  3:221025/1‑124 (MQ=255)                                                                      ttttGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTg                                                                 <  3:350701/149‑1 (MQ=255)                                                                        ttGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTgg                                                                <  1:140402/148‑1 (MQ=255)                                                                        ttGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTgg                                                                <  8:215823/148‑1 (MQ=255)                                                                            ccTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGc                                                                                                                         >  6:372862/1‑88 (MQ=255)                                                                            ccTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGc                                                                                                                         <  5:372862/88‑1 (MQ=255)                                                                             cTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGGTGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCGAAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTg                                                          >  3:246196/1‑149 (MQ=255)                                                                             cTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCgta                                                                                                   <  5:178583/108‑1 (MQ=255)                                                                             cTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCgta                                                                                                   >  6:178583/1‑108 (MQ=255)                                                                              tgttgtTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTgg                                                                <  7:261431/142‑1 (MQ=255)                                                                              tgttgtTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTgg                                                                >  8:261431/1‑142 (MQ=255)                                                                                   ttgCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCa                                                    <  1:315340/149‑1 (MQ=255)                                                                                     gCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTg                                                          >  6:248652/1‑141 (MQ=255)                                                                                     gCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTg                                                          <  5:248652/141‑1 (MQ=255)                                                                                             aaaGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTc                                          <  3:198107/149‑1 (MQ=255)                                                                                             aaaGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTc                                          <  2:303669/149‑1 (MQ=255)                                                                                              aaGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTc                                          >  3:195907/1‑148 (MQ=255)                                                                                                ggCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAg                                       >  7:132313/1‑149 (MQ=255)                                                                                                    atatTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTgg                                   >  8:232946/1‑149 (MQ=255)                                                                                                       ttGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATAtt                                              >  5:76360/1‑135 (MQ=255)                                                                                                       ttGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATAtt                                              <  6:76360/135‑1 (MQ=255)                                                                                                         gAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCgta                                                                                                   >  5:253255/1‑80 (MQ=255)                                                                                                         gAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCgta                                                                                                   <  6:253255/80‑1 (MQ=255)                                                                                                          aaaTCGAGGTGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCGAAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGaa                              <  4:246196/148‑1 (MQ=255)                                                                                                          aaaTCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAg                             >  7:26950/1‑149 (MQ=255)                                                                                                               gagaTGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCAt                         <  8:26950/148‑1 (MQ=255)                                                                                                               gagaTGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATg                        <  3:258303/149‑1 (MQ=255)                                                                                                                agaTGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGa                       <  8:132313/149‑1 (MQ=255)                                                                                                                agaTGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGa                       >  3:228125/1‑149 (MQ=255)                                                                                                                       aaTTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCg                >  5:289951/1‑149 (MQ=255)                                                                                                                         ttGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTg                                                                 >  2:240747/1‑98 (MQ=255)                                                                                                                         ttGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTg                                                                 <  1:240747/98‑1 (MQ=255)                                                                                                                            gTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGAt                      >  1:381969/1‑138 (MQ=255)                                                                                                                            gTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCACCTGGCGGAAGTCATGAt                      <  2:381969/138‑1 (MQ=255)                                                                                                                            gTTGTAG‑GCCAGATAAGACGCATAAGCGTTGTATCTGGCACTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCa                                                    <  6:100892/108‑1 (MQ=18)                                                                                                                            gTTGTAG‑GCCAGATAAGACGCATAAGCGTTGTATCTGGCACTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCa                                                    >  5:100892/1‑108 (MQ=18)                                                                                                                             ttGTAG‑GCCAGATAAGACGCATAAGCGTTGTATCTGGCACTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTAt          >  4:281283/1‑149 (MQ=37)                                                                                                                                  gTGCCAGATACAGCGGGTATGCGTCTTAGCTGGG‑CTACAACGATTACATGGCGTAGTAATACGCTTTCGCCTGCTGCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATCGCGa     >  8:317997/1‑149 (MQ=255)                                                                                                                                  gTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATCGCGa     <  7:232946/149‑1 (MQ=255)                                                                                                                                      cAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATCGCGACGc  >  8:275128/1‑148 (MQ=255)                                                                                                                                            |                                                                                                                                               GCACACGCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGC‑CTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATCGCGACGC  >  NC_000913/2348650‑2348930

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