Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,259,693	A→T	8.1%	intergenic (‑397/+26)	psuK ← / ← fruA	pseudouridine kinase/fused fructose‑specific PTS enzymes: IIBcomponent/IIC components

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,259,693	0	A	T	8.1%	204.0 / 8.1	62	intergenic (‑397/+26)	psuK/fruA	pseudouridine kinase/fused fructose‑specific PTS enzymes: IIBcomponent/IIC components
Reads supporting (aligned to +/- strand):  ref base A (31/26);  new base T (3/2);  total (34/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 = 7.10e-01

GAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGCGTAATAGCGCCAG  >  NC_000913/2259547‑2259840                                                                                                                                                   |                                                                                                                                                     gAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAAtt                                                                                                                                                    >  1:224020/1‑149 (MQ=255)      taatCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGc                                                                                                                                              >  4:381388/1‑149 (MQ=255)         tCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCt                                                                                                                                           <  1:206488/149‑1 (MQ=255)           tGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGt                                                                                                                                         >  1:86167/1‑149 (MQ=255)                  gAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAacacac                                                                                                                                  >  4:79611/1‑149 (MQ=255)                         tAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTttat                                                                                                                           <  2:446906/149‑1 (MQ=255)                             cTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACg                                                                                                                       >  3:64158/1‑149 (MQ=255)                                        atgaACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGGCAGAAATTTC‑CTGTCCTGTAac                                                                                                                                      >  2:15079/1‑123 (MQ=38)                                        atgaACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGGCA‑GAAATTTCCTGTCCTGTAac                                                                                                                                      <  1:15079/123‑1 (MQ=38)                                           aaCGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGTCAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACt                                                                                                         <  4:197969/149‑1 (MQ=255)                                           aaCGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACt                                                                                                         >  1:428226/1‑149 (MQ=255)                                           aaCGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACt                                                                                                         <  1:494811/149‑1 (MQ=255)                                            aCGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTg                                                                                                        <  2:329117/149‑1 (MQ=255)                                             cGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGc                                                                                                       <  2:431787/149‑1 (MQ=255)                                              gTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCg                                                                                                      >  1:315405/1‑149 (MQ=255)                                               tttGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGt                                                                                                     <  2:24854/149‑1 (MQ=255)                                                   aaGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCAc                                                                                                 <  4:47915/149‑1 (MQ=255)                                                     gTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTttatt                                                                                                                          >  3:204524/1‑122 (MQ=255)                                                     gTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTttatt                                                                                                                          <  4:204524/122‑1 (MQ=255)                                                     ctttCTACTGAAAGATAAACCCATTTAAATTCATAAAATTACAATCACAAACCGGTCATTCTCTTTGCGCATCTTTAATATTACATTACCGGAAATTT‑CTGCCCTGTCACACATATTTTATTACGCTCCTTCCGCTACTGTGTCCACtt                                                                                               >  4:513571/2‑149 (MQ=255)                                                             tGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGt                                                                                       >  4:109353/1‑149 (MQ=255)                                                             tGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGt                                                                                       >  4:55926/1‑149 (MQ=255)                                                             tGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGt                                                                                       <  4:73035/149‑1 (MQ=255)                                                              gAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGtt                                                                                      <  3:505925/149‑1 (MQ=255)                                                                 aGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCa                                                                                   <  4:338001/149‑1 (MQ=255)                                                                      aaaCCATTTAAATTCATAAAATTACATTTACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAg                                                                              >  4:203656/1‑149 (MQ=255)                                                                       aaCCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAgg                                                                             >  3:260109/1‑149 (MQ=255)                                                                        aCCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACATGACAGGAAATTT‑CTTCCCTGTAACACACATTTTATTATGCTGCATTCGCTACTGCGTCCCCTTCCGGACGTTTCAGGAAGGc                                                                            >  4:45231/1‑149 (MQ=255)                                                                        aCCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGc                                                                            <  1:546082/149‑1 (MQ=255)                                                                           aTTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTcc                                                                                                   <  4:248711/123‑1 (MQ=255)                                                                           aTTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTcc                                                                                                   >  3:248711/1‑123 (MQ=255)                                                                           aTTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATa                                                                         >  2:430444/1‑149 (MQ=255)                                                                            tttAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAg                                                                        <  4:6466/149‑1 (MQ=255)                                                                             ttAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAgg                                                                       <  3:357603/149‑1 (MQ=255)                                                                                 aTTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCaa                                                                   >  3:148388/1‑149 (MQ=255)                                                                                  ttCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCaaa                                                                  <  1:219262/149‑1 (MQ=255)                                                                                    cATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAAcc                                                                >  3:62119/1‑149 (MQ=255)                                                                                     aTAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAAccc                                                               >  2:438811/1‑149 (MQ=255)                                                                                          aTTACATTCACAAACCGGTCTTTCTCTTTGCGCTCATCAAACGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGccac                                                          >  3:370236/1‑149 (MQ=255)                                                                                               aTTCACAACCCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCg                                                     <  2:280244/149‑1 (MQ=255)                                                                                               aTTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCg                                                     >  1:12859/1‑149 (MQ=255)                                                                                                   acaAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAg                                                                              >  1:571690/1‑120 (MQ=255)                                                                                                   acaAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAg                                                                              <  2:571690/120‑1 (MQ=255)                                                                                                    caAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCg                                                >  1:81967/1‑149 (MQ=255)                                                                                                      aaCCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGc                                              >  2:556891/1‑149 (MQ=255)                                                                                                           gTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCaataa                                         >  1:239597/1‑149 (MQ=255)                                                                                                             cATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAAtt                                       <  3:63396/149‑1 (MQ=255)                                                                                                              aTTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTg                                      >  3:186960/1‑149 (MQ=255)                                                                                                                    tttGCGCTCATCAAATGTTACAGGGCAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCa                                <  2:153387/149‑1 (MQ=255)                                                                                                                    tttGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCa                                >  1:575239/1‑149 (MQ=255)                                                                                                                    tttGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCa                                <  2:275898/149‑1 (MQ=255)                                                                                                                       gcgcTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGt                             <  4:348523/149‑1 (MQ=255)                                                                                                                         gcTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTaa                           >  2:278380/1‑149 (MQ=255)                                                                                                                            catcaAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAAccc                        <  1:4321/149‑1 (MQ=255)                                                                                                                                   tGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCg                 <  2:224020/149‑1 (MQ=255)                                                                                                                                     ttACAGGGCAGAAATTTC‑CTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGc               <  1:360185/149‑1 (MQ=255)                                                                                                                                     ttACAGGGCAGAAATTTC‑CTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGc               >  4:504131/1‑149 (MQ=255)                                                                                                                                     ttACAGGGCAGAAATTTC‑CTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGc               >  2:360185/1‑149 (MQ=255)                                                                                                                                     ttACAGGGCAGAAATTTC‑CTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGc               <  3:504131/149‑1 (MQ=255)                                                                                                                                          ggacaggaAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCGTAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGCGTAAt          >  4:320742/1‑149 (MQ=255)                                                                                                                                               ggaAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGCGTAATAgcgc     >  1:359456/1‑149 (MQ=255)                                                                                                                                                gaAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACTCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGCGTAATAGCGcc    >  2:99683/1‑149 (MQ=255)                                                                                                                                                  aaTTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGCGTAATAGCGCCAg  <  1:556891/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     GAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGCGTAATAGCGCCAG  >  NC_000913/2259547‑2259840

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