Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	538,646	C→T	6.8%	A338A (GCC→GCT)	ybbW →	putative allantoin transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	538,646	0	C	T	6.8%	184.8 / 2.6	59	A338A (GCC→GCT)	ybbW	putative allantoin transporter
Reads supporting (aligned to +/- strand):  ref base C (26/29);  new base T (2/2);  total (28/31)
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.78e-01

ACTACGGCGCTGATACCTGGAACGTGCTGGATATTGTTCAGCGTTGGGACAGCCTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTGGAAATTAATGGAAAATCAGGACAGCATTTATCTTTTCCTCGATATTATCGGCGGAATGCTTGGTCCGGT  >  NC_000913/538498‑538792                                                                                                                                                     |                                                                                                                                                   aCTACGGCGCTGATACCTGGAACGTGCTGGATATTGTTCAGCGTTGGGACAGCCTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTgcc                                                                                                                                                    <  8:166690/149‑1 (MQ=255) aCTACGGCGCTGATACCTGGAACGTGCTGGATATTGTTCAGCGTTGGGACAGCCTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTgcc                                                                                                                                                    >  4:59743/1‑149 (MQ=255)   tACGGCGCTGATACCTGGAACGTGCTGGATATTGTTCAGCGTTGGGACAGCCTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTgccgc                                                                                                                                                  >  2:27158/1‑149 (MQ=255)       gcgcTGATACCTGGAACGTGCTGGATATTGTTCAGCGTTGGGACAGCCTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCAtt                                                                                                                                              >  6:200082/1‑149 (MQ=255)                 cTGGAACGTGCTGGATATTGTTCAGCGTTGGGACAGCCTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACa                                                                                                                                     >  1:24702/1‑148 (MQ=255)                 cTGGAACGTGCTGGATATTGTTCAGCGTTGGGACAGCCTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACa                                                                                                                                     >  6:182157/1‑148 (MQ=255)                           cTGGATATTGTTCAGCGTTGGGACAGCCTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCta                                                                                                                          <  2:529/149‑1 (MQ=255)                           cTGGATATTGTTCAGCGTTGGGACAGCCTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCt                                                                                                                           <  6:102007/148‑1 (MQ=255)                           cTGGATATTGTTCAGCGTTGGGACAGCCTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCt                                                                                                                           <  2:4935/148‑1 (MQ=255)                           cTGGATATTGTTCAGCGTTGGGACAGCCTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCAAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCta                                                                                                                          >  6:361963/1‑149 (MQ=255)                                tatTGTTCAGCGTTGGGACAGCCTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATaaa                                                                                                                      >  8:89240/1‑148 (MQ=255)                                  ttGTTCAGCGTTGGGACAGCCTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAAc                                                                                                                   <  7:89240/149‑1 (MQ=255)                                     ttCAGCGTTGGGACAGCCTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGc                                                                                                                >  6:198446/1‑149 (MQ=255)                                     ttCAGCGTTGGGACAGCCTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGc                                                                                                                >  6:203356/1‑149 (MQ=255)                                     ttCAGCGTTGGGACAGCCTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGc                                                                                                                <  2:24702/149‑1 (MQ=255)                                     ttCAGCGTTGGGACAGCCTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGc                                                                                                                <  5:200082/149‑1 (MQ=255)                                         gCGTTGGGACAGCCTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCTGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAACCTGACCTATAAAAACGGCGTAc                                                                                                            <  3:24331/149‑1 (MQ=255)                                            ttGGGACAGCCTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGa                                                                                                         <  1:27158/149‑1 (MQ=255)                                            ttGGGACAGCCTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGa                                                                                                         >  2:352204/1‑149 (MQ=255)                                                  cAGCCTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCa                                                                                                   <  1:352204/149‑1 (MQ=255)                                                   aGCCTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCAGTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAg                                                                                                  >  5:20921/1‑149 (MQ=255)                                                     ccTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGta                                                                                                <  3:59743/149‑1 (MQ=255)                                                        gTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGtatta                                                                                             <  5:329861/149‑1 (MQ=255)                                                           cGCCTCGTTCTTAGCGTTACTGTTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCAAAGACAAAACTGACCTATAAAAAAGGCGTACTGATTGCCAGTATTATc                                                                                           <  2:140267/148‑1 (MQ=255)                                                             ccTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGc                                                                                        >  6:132865/1‑149 (MQ=255)                                                                         cGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCTGTAATATTAATCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGc                                                                            <  1:301217/149‑1 (MQ=255)                                                                               tGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAACCGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTg                                                                       >  4:343576/1‑148 (MQ=255)                                                                               tGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTg                                                                       >  3:124176/1‑148 (MQ=255)                                                                                 gTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCTGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTGGa                                                                     >  8:275940/1‑148 (MQ=255)                                                                                 gTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTGGa                                                                     >  2:352396/1‑148 (MQ=255)                                                                                 gTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTGGa                                                                     >  7:268635/1‑148 (MQ=255)                                                                                     ttCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCCGTTATCAGAGTGCCGCCATTGCACCGACAAAACTGACCTATACAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTGGAAAtt                                                                 >  5:354611/1‑148 (MQ=255)                                                                                          atgaCAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCTGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTGGAAATTAAt                                                              <  7:275940/146‑1 (MQ=255)                                                                                          atgaCAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTGGAAATTAAt                                                              <  5:132865/146‑1 (MQ=255)                                                                                          atgaCAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTGGAAATTAAt                                                              <  2:284230/146‑1 (MQ=255)                                                                                          atgaCAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTGGAAATTAATg                                                             <  5:203356/147‑1 (MQ=255)                                                                                          atgaCAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTGGAAATTAATg                                                             <  5:198446/147‑1 (MQ=255)                                                                                               aaCTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGc                                                                                        >  1:139581/1‑115 (MQ=255)                                                                                               aaCTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGc                                                                                        <  2:139581/115‑1 (MQ=255)                                                                                                 cTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTGGAAATTaa                                                               <  7:347104/138‑1 (MQ=255)                                                                                                 cTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTGGAAATTaa                                                               >  8:347104/1‑138 (MQ=255)                                                                                                 cTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTGGAAATTAATGGAAAATCAg                                                    <  8:165044/149‑1 (MQ=255)                                                                                                    tctcCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTGGAAATTAATGGAAAATCAGGa                                                  >  5:340690/1‑148 (MQ=255)                                                                                                        cACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTGGAAATTAATGGAAAATCAGGACAGCa                                             >  5:251139/1‑149 (MQ=255)                                                                                                         aCTCACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTTCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTGGAAATTAATGGAAAATCAGGACAGCAt                                            <  6:86372/149‑1 (MQ=255)                                                                                                         aCTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTGGAAATTAATGGAAAATCAGGACAGCAt                                            <  2:141091/149‑1 (MQ=255)                                                                                                         aCTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTGGAAATTAATGGAAAATCAGGACAGCAt                                            >  6:237400/1‑149 (MQ=255)                                                                                                         aCTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTGGAAATTAATGGAAAATCAGGACAGCAt                                            >  6:6372/1‑149 (MQ=255)                                                                                                                 gACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTGGAAATTAATGGAAAATCAGGACAGCATTTATCttt                                    <  3:343576/149‑1 (MQ=255)                                                                                                                  aCCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTGGAAATTAATGGAAAATCAGGACAGCATTTATCtttt                                   >  8:81833/1‑149 (MQ=255)                                                                                                                      gTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTGGAAATTAATGGAAAATCAGGACAGCATTTATCTTTTCCTc                               >  5:213685/1‑149 (MQ=255)                                                                                                                         atatTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTGGAAATTAATGGAAAATCAGGACAGCATTTATCTTTTCCTCGa                             <  6:251139/148‑1 (MQ=255)                                                                                                                          tattatTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTGGAAATTAATGGAAAATCAGGACAGCATTTATCTTTTCCTCGata                           <  5:237400/149‑1 (MQ=255)                                                                                                                           attattCCAGCCGGTTATCAGATTGCCGCCATTCCACCGACAAAACTCACCTATAAAAACGGCGTACTGAATGCCAGTATTAACAGCTAGCTGATCTGCCCGTGGAAATTAATGGAAAATCAGCACAGCATTTATCTTTTCCACGatat                          >  5:299806/1‑149 (MQ=255)                                                                                                                               ttCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTGGAAATTAATGGAAAATCAGGACAGCATTTATCTTTTCCTCGATATTATc                      >  3:161057/1‑149 (MQ=255)                                                                                                                                  cagcagGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTGGAAATTAATGGAAAATCAGGACAGCATTTATCTTTTCCTCGATATTATcggc                   >  2:105280/1‑149 (MQ=255)                                                                                                                                  cAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTGGAAATTAATGGAAAATCAGGACAGCATTTATCTTTTCCTCGATATTATcggc                   <  7:81833/149‑1 (MQ=255)                                                                                                                                  cAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTGGAAATTAATGGAAAATCAGGACAGCATTTATCTTTTCCTCGATATTATcgg                    <  5:182157/148‑1 (MQ=255)                                                                                                                                    gCCGGTTATCAGATTGCTGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTGGAAATTAATGGAAAATCAGGACAGCATTTATCtttt                                   <  2:62310/131‑1 (MQ=255)                                                                                                                                    gCCGGTTATCAGATTGCTGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTGGAAATTAATGGAAAATCAGGACAGCATTTATCtttt                                   >  1:62310/1‑131 (MQ=255)                                                                                                                                           aTCAGATTGCCGCCATTACACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTGGAAATTAATGGAAAATCAGGACAGCATTTATCTTTTCCTCGATATTATCGGCGGAATGCt           <  8:326424/148‑1 (MQ=255)                                                                                                                                              aGATTGCCGCATTTGTACCGTCAAAAGTGAACTATAATCACGGCGTACTGATTGCCAGTATTATCAGGTTGCTGATCTGCCCGTGGAAATTAATGGAAAATCAGGACAGCATTTATCTTTTCCTCGATATTATCGGCGGAATGCTTgg        <  6:299806/148‑1 (MQ=255)                                                                                                                                                    ccgccATTGCCCCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTGGAAATTAATGGAAAATCAGGACAGCATTTATCTTTTCCTCGATATTATCGGAGGAATGCTTGGTCCGGt  >  1:318603/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                   ACTACGGCGCTGATACCTGGAACGTGCTGGATATTGTTCAGCGTTGGGACAGCCTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTGGAAATTAATGGAAAATCAGGACAGCATTTATCTTTTCCTCGATATTATCGGCGGAATGCTTGGTCCGGT  >  NC_000913/538498‑538792

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