Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,552,306	A→G	5.1%	E224G (GAA→GGA)	uxuA →	mannonate hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,552,306	0	A	G	5.1%	245.0 / 2.9	78	E224G (GAA→GGA)	uxuA	mannonate hydrolase
Reads supporting (aligned to +/- strand):  ref base A (36/37);  new base G (2/2);  total (38/40)
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.33e-01

GTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGTTTCCACCATTGAAGATATGCAGTGGATGGTTGATACCGTAAACAGCATGGCAAACGGTTTTACCATGTGCACCGGT  >  NC_000913/4552159‑4552454                                                                                                                                                    |                                                                                                                                                     gTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTgaa                                                                                                                                                     <  3:109208/149‑1 (MQ=255)    aCATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTgaagaa                                                                                                                                                  <  7:385060/149‑1 (MQ=255)    aCATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTgaagaa                                                                                                                                                  <  1:392845/149‑1 (MQ=255)    aCATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTgaagaa                                                                                                                                                  <  2:308630/149‑1 (MQ=255)     catcatTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATAATTCCAGTTGCTGaagaag                                                                                                                                                 >  2:294705/1‑149 (MQ=255)      atcatTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGt                                                                                                                                                >  3:22593/1‑149 (MQ=255)      atcatTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGt                                                                                                                                                >  5:205075/1‑149 (MQ=255)       tcatTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTc                                                                                                                                               >  6:50255/1‑149 (MQ=255)        catTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCg                                                                                                                                              >  7:105225/1‑149 (MQ=255)        catTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCg                                                                                                                                              >  5:296651/1‑149 (MQ=255)            gCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGgcg                                                                                                                                           <  7:371782/148‑1 (MQ=255)             cTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGgcgt                                                                                                                                          <  3:27252/148‑1 (MQ=255)             cTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGgcgt                                                                                                                                          <  5:50255/148‑1 (MQ=255)                 tCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTa                                                                                                                                     <  1:116681/149‑1 (MQ=255)                   ttCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATg                                                                                                                                   >  3:348784/1‑149 (MQ=255)                     ccGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATgg                                                                                                                                  <  2:28105/148‑1 (MQ=255)                          cgcgGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTTCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTc                                                                                                                            <  8:194430/149‑1 (MQ=255)                            cgGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCAc                                                                                                                          <  8:119347/149‑1 (MQ=255)                               aagaagGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCg                                                                                                                       >  6:329386/1‑149 (MQ=255)                                      cTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATc                                                                                                                <  5:82124/149‑1 (MQ=255)                                      cTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATc                                                                                                                <  4:348784/149‑1 (MQ=255)                                           cccTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATccgccg                                                                                                           >  1:215771/1‑149 (MQ=255)                                            ccTCGACCAGTTCCTTAAACAAATGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCCCCAGGACGATCcgccgc                                                                                                          >  6:92001/1‑149 (MQ=255)                                                gACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATTTCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGACAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCgcgc                                                                                                        <  8:185085/147‑1 (MQ=255)                                                        ccGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTc                                                                                               <  4:225232/148‑1 (MQ=255)                                                                  cTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCTGACGATCCGCCGCGCCCGATCCTCGGCCTGCcgcg                                                                                    <  2:215771/149‑1 (MQ=255)                                                                   tGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTACAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCAGCGCCCGATCCTCGGCCTGCcgcg                                                                                    >  3:359251/1‑148 (MQ=255)                                                                     gAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCAt                                                                                 >  6:160387/1‑149 (MQ=255)                                                                      aGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCAtt                                                                                <  7:71906/149‑1 (MQ=255)                                                                        cTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGt                                                                              >  7:211571/1‑149 (MQ=255)                                                                         tGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGtt                                                                             >  3:309077/1‑149 (MQ=255)                                                                          gTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGttt                                                                            >  7:133171/1‑149 (MQ=255)                                                                          gTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGttt                                                                            >  2:364246/1‑149 (MQ=255)                                                                           tACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGTTTc                                                                           >  4:294585/1‑149 (MQ=255)                                                                           tACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGTTTc                                                                           <  5:219648/149‑1 (MQ=255)                                                                           tACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCAGTTCACCCGGACGATCCGCAGCGCCCGATCCTCGGCCTGCCGCGCATAGTTTc                                                                           >  3:78423/1‑149 (MQ=255)                                                                              aaaGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGTCCGATCCTCGGCCTGCCGCGCATTGTTTccac                                                                        >  1:181123/1‑149 (MQ=255)                                                                              aaaGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGTTTccac                                                                        <  6:97836/149‑1 (MQ=255)                                                                                  atatCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGTTTCCACCAtt                                                                    >  4:367878/1‑149 (MQ=255)                                                                                   tatCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGACTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGTTTCCACCAtt                                                                    >  1:278655/1‑148 (MQ=255)                                                                                     tCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTAGGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGTTTCCACCATTGaa                                                                 <  1:294705/149‑1 (MQ=255)                                                                                       gACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCATTTGCTGCAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGTTTCCACCATTGAAg                                                                <  5:92001/148‑1 (MQ=255)                                                                                       gACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGATGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGTTTCCGCCATTg                                                                   >  5:143034/1‑145 (MQ=255)                                                                                        aCAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGGAGAAGTCGgcgtgcgt                                                                                                                                      <  8:173749/77‑1 (MQ=255)                                                                                        aCAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGGAGAAGTCGgcgtgcgt                                                                                                                                      <  8:179102/77‑1 (MQ=255)                                                                                        aCAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGGAGAAGTCGgcgtgcgt                                                                                                                                      >  7:173749/1‑77 (MQ=255)                                                                                        aCAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGGAGAAGTCGgcgtgcgt                                                                                                                                      >  7:179102/1‑77 (MQ=255)                                                                                        aCAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTc                                                                                               <  3:314384/116‑1 (MQ=255)                                                                                        aCAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTc                                                                                               >  4:314384/1‑116 (MQ=255)                                                                                        aCAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGTTTCCACCATTGAAGa                                                               >  1:308142/1‑148 (MQ=255)                                                                                         cAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGTTTCCACCATTGAAGata                                                             <  1:17952/149‑1 (MQ=255)                                                                                          aaaGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGTTTCCACCATTGAAGatat                                                            <  5:160387/149‑1 (MQ=255)                                                                                            aGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGTTTCCACCATTGAAGATATGc                                                          <  8:107349/149‑1 (MQ=255)                                                                                                aaGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGTTTCCACCATTGAAGATATGCAGTg                                                      <  4:22593/149‑1 (MQ=255)                                                                                                 aGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGgcgt                                                                                                                                          <  6:294990/64‑1 (MQ=255)                                                                                                 aGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGgcgt                                                                                                                                          >  5:294990/1‑64 (MQ=255)                                                                                                  gCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGTTTCCACCATTGAAGATATGCAGTGGa                                                    <  1:364246/149‑1 (MQ=255)                                                                                                   cTGCGCGAAAACTTTGCCGTCTTCCTGAAACCGAATAATCCAGTTGCTGAAGAACTCGGCGCGCGTATGGCTGGCCACCCGGACGCTCCGCCGCGCCCGATCCTCGGCCGGCCGCGCCTTGTTTCCCCCATTGAAGATATGCAGTGGac                                                   <  6:142837/149‑2 (MQ=255)                                                                                                     gcgcGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGTTTCCACCATTGAAGATATGCAGTGGATgg                                                 >  8:9918/1‑149 (MQ=255)                                                                                                     gcgcGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGTTTCCACCATTGAAGATATGCAGTGGATgg                                                 >  1:385375/1‑149 (MQ=255)                                                                                                      cgcgAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGTCCGATCCTCGGCCTGCCGCGCATTGTTTCCACCATTGAAGATATGCAGTGGATGGt                                                <  2:181123/149‑1 (MQ=255)                                                                                                      cgcgAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGTTTCCACCATTGAAGATATGCAGTGGATGGt                                                <  3:294585/149‑1 (MQ=255)                                                                                                       gcgAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGTTTCCACCATTGAAGATATGCAGTGGATGGtt                                               >  4:145957/1‑149 (MQ=255)                                                                                                       gcgAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGTTTCCACCATTGAAGATATGCAGTGGATGGtt                                               <  7:164207/149‑1 (MQ=255)                                                                                                       gcgAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGGGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGTTTCCACCATTGAAGATATGCAGTGGATGGtt                                               <  6:205075/149‑1 (MQ=255)                                                                                                         gAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTTAAGAAGTCGGCGTGCGTATGGCAGTTCACCCGGACCATCCGCCGCGCCCGATCCTCGGCCTCCCGCGCATTGTTTCCACCATTCAAGATATGCAGTGGATGGATGa                                             >  6:373149/1‑149 (MQ=255)                                                                                                         gAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGGCGCGCATTGTTTCCACCATTGAAGATATGCAGTGGATGGTTGa                                             <  2:225708/149‑1 (MQ=255)                                                                                                                 tGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGTTTCCACCATTGAAGATATGCAGTGGATGGTTGATACCGt                                       <  5:329386/147‑1 (MQ=255)                                                                                                                    cGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACACGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGTTTCCACCATTGAAGATATGCAGAGGATGGTTGATACCGTaaa                                    >  8:294542/1‑147 (MQ=255)                                                                                                                         tCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGTTTCCACCATTGAAGATATGCAGTGGATGGTTGATACCGTAAACAGCATg                             >  4:352570/1‑149 (MQ=255)                                                                                                                                aGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGTTTCCACCATTGaa                                                                 >  8:59648/1‑106 (MQ=255)                                                                                                                                aGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGTTTCCACCATTGaa                                                                 <  7:59648/106‑1 (MQ=255)                                                                                                                                   gATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGTTTCCACCATTGAAGATATGCAGTGGATGGTTGATACCGTAAACAGCATGGCAAACGGtt                   >  6:314609/1‑149 (MQ=255)                                                                                                                                     ttattCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGTTTCCACCATTGAAGATATGCAGTGGATGGTTGATACCGTAAACAGCATGGCAAACGGttt                  >  3:268134/1‑148 (MQ=255)                                                                                                                                       attCCAGTTGCTGAAGAAGTCGGCGTGCGTATGACTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGTTTCCACCATTGAAGATATGCAGTGGATGGTTGATACCGTAAACAGCATGGCAAACGGTTTTa                <  2:278655/148‑1 (MQ=255)                                                                                                                                          cctgTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCAGCCGGACGATCCGCCGCGCCCGATCCTCTGCCTGCCGCGCATTGTTTCCACCTTTGAAGATATGCAGTGGATGGTTGATACCGTAAACAGCATTGCAAACGGTTTTACCa             <  4:78423/145‑1 (MQ=255)                                                                                                                                          ccAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGTTTCCACCATTGAAGATATGCAGTGGATGGTTGATACCGTAAACAGCATGGCaa                         >  8:110016/1‑136 (MQ=255)                                                                                                                                          ccAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGTTTCCACCATTGAAGATATGCAGTGGATGGTTGATACCGTAAACAGCATGGCaa                         <  7:110016/136‑1 (MQ=255)                                                                                                                                          ccAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGTTTCCACCATTGAAGATATGCAGTGGATGGTTGATACCGTAAACAGCATGGCAAACGGTTTTACCa             <  3:145957/148‑1 (MQ=255)                                                                                                                                           cAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGTTTCCACCATTGAAGATATGCAGTGGATGGTTGATTCCGTAAACAGCATGGCAAACGGTTTTACCAtg           >  7:140811/1‑149 (MQ=255)                                                                                                                                           cAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGTTTCCACCATTGAAGATATGCAGTGGATGGTTGATACCGTAAACAGCATGGCAAATGGTTTTACCAt            <  1:5149/148‑1 (MQ=255)                                                                                                                                           cAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGTTTCCACCATTGAAGATATGCAGTGGATGGTTGATACCGTAAACAGCATGGCAAACGGTTTTACCAtg           >  8:237715/1‑149 (MQ=255)                                                                                                                                              ttGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGTTTCCACCATTGAAGATATGCAGTGGATGGTTGATACCGTAAACAGCATGGCAAACGGTTTTACCATGTGc        <  7:9918/149‑1 (MQ=255)                                                                                                                                                    aagaagTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGTTTCCACCATTGAAGATATGCAGTGGATGGTTGATACCGTAAACAGCATGGCAAACGGTTTTACCATGTGCACCgg   <  8:105225/148‑1 (MQ=255)                                                                                                                                                    aagaagTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCAGAGCATTGTTTCGACCATTGGAGATATGCAGTGGATGGATGATACCGTAAACAGCATGGCAAACGGTTTTACCATGTGCACCGGt  >  5:368528/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTTCACCCGGACGATCCGCCGCGCCCGATCCTCGGCCTGCCGCGCATTGTTTCCACCATTGAAGATATGCAGTGGATGGTTGATACCGTAAACAGCATGGCAAACGGTTTTACCATGTGCACCGGT  >  NC_000913/4552159‑4552454

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