Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,807,043	T→A	65.8%	L225Q (CTG→CAG)	pfkB →	6‑phosphofructokinase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,807,043	0	T	A	65.8%	92.6 / 74.2	84	L225Q (CTG→CAG)	pfkB	6‑phosphofructokinase II
Reads supporting (aligned to +/- strand):  ref base T (15/13);  new base A (28/27);  total (44/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 = 1.00e+00

AATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAA  >  NC_000913/1806900‑1807188                                                                                                                                                |                                                                                                                                                   aaTTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGT‑c                                                                                                                                              >  7:228528/1‑149 (MQ=255) aaTTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑c                                                                                                                                              >  7:53600/1‑149 (MQ=255)      gTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCATTGCGCTTGTGAATCGCGAACTCACCCAGCCGCACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCTTTTCCCAGGGT‑CCACaa                                                                                                                                         <  4:13221/149‑1 (MQ=255)      gTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACaa                                                                                                                                         <  8:345947/149‑1 (MQ=255)      gTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATAGCGAACTCACCCCGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCGTGGT‑ccca                                                                                                                                           >  7:273902/1‑145 (MQ=255)         aCATCGAGTTGGTTAAGCCTAACCAATAAGATCTCAGTGGGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGTAAATCTTTAATAGCGGCTATGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACATgg                                                                                                                                       <  4:175475/148‑1 (MQ=255)         aCATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAg                                                                                                                                        >  3:309895/1‑147 (MQ=255)         aCATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAg                                                                                                                                        >  1:90075/1‑147 (MQ=255)           aTCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAGGAgc                                                                                                                                    >  2:269499/1‑149 (MQ=255)              gAGTTGGATAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGCCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTGTCCCAGGGT‑CCACAAGGAGCGCt                                                                                                                                 <  8:309972/149‑1 (MQ=255)               aGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCCCCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGT‑CCACAAGGAGCGCTg                                                                                                                                >  2:295540/1‑149 (MQ=255)                gTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAGGAGCGCTgg                                                                                                                               <  5:31872/149‑1 (MQ=255)                   ggTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGCTGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAACCGGGTTGTCGTTTCCCTGGGT‑CCACAAGGAGCGCTggggg                                                                                                                            >  2:46010/1‑147 (MQ=255)                   ggTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGT‑CCACAAGGAGCGCTGGgtg                                                                                                                            >  4:250882/1‑149 (MQ=255)                    gTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAGGAGCGCTGGgtgt                                                                                                                           >  8:61607/1‑149 (MQ=255)                     ttAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCCGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGT‑CCACAAGGAGCGCTGGGTGtt                                                                                                                          >  6:257683/1‑149 (MQ=255)                     ttAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGTCCTCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAGGAGCGCTGGGTGtt                                                                                                                          >  5:97836/1‑149 (MQ=255)                     ttAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAGGAGCGCTGGGTGtt                                                                                                                          >  7:138699/1‑149 (MQ=255)                     ttAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAGGAGCGCTGGGTGtt                                                                                                                          >  5:78190/1‑149 (MQ=255)                       aaGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAGGAGCGCTGGGTGTTg                                                                                                                         <  8:53600/148‑1 (MQ=255)                        aGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAGGAGCGCTGGGTGTTGa                                                                                                                        <  7:61607/148‑1 (MQ=255)                         gCCTAACCAAAAAGAACTCAGGGCGCTGGGTAATCGCGAACTCACCCCGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTAGTTTCCCTGGGT‑ACACAAGGAGCGCTGGGTGTTGATa                                                                                                                      <  3:12729/149‑1 (MQ=255)                            tAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAGGAGCGCTCGGTGTTGATAGt                                                                                                                    >  1:176157/1‑148 (MQ=255)                             aaCCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTg                                                                                                                   >  7:327517/1‑148 (MQ=255)                                cataaaGAACTCAGTGAGCTGGTGAATCGCGAACTCACCCCTCCGGACGATGTCCGCAAATCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAATCGGGTTGTCGTTTCCCTGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGaaaa                                                                                                               <  3:230627/146‑1 (MQ=255)                                cAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGaaaa                                                                                                               >  7:262756/1‑149 (MQ=255)                                cAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGAAGGAAAACGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGaaaa                                                                                                               >  2:23003/1‑149 (MQ=255)                                 aaaaaGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCCGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAAc                                                                                                              >  7:57000/1‑149 (MQ=255)                                 aaaaaGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAAc                                                                                                              >  1:266271/1‑149 (MQ=255)                                  acaaGACCTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCTACCGCGTGGGCGTGTCCCAGGGTACCAC‑AGCAGCGCGGGGGGTTGATAGTGatcact                                                                                                             >  8:239513/3‑144 (MQ=255)                                  aaaaGAACTCAGTGCGCTGGTGAATCGCGAACTCCCCCAGCCGGACGATGTCCGCAAAGCCGCGAAGGAAAACGTTAATAGCGTCAAGGCCAAACGGGTTGTCGTTGCCCAGGGT‑CCACACGGAGCGCTGGGTGTTGAAATTGaaaaat                                                                                                             >  8:239514/1‑147 (MQ=255)                                    aaGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGt                                                                                                           <  4:309895/149‑1 (MQ=255)                                          tCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGT‑CCACAAGGAGCGCTGGGTGTTGa                                                                                                                        >  6:163216/1‑130 (MQ=255)                                          tCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGT‑CCACAAGGAGCGCTGGGTGTTGa                                                                                                                        <  5:163216/130‑1 (MQ=255)                                                   tGGTGAATCGCGAACTCACCCCGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGT‑CCCCAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGcca                                                                                            >  7:166956/1‑149 (MQ=255)                                                        aaTCGCGAACTCACCCCGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGGCGTTTCCCAGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCcaccac                                                                                        <  3:142418/148‑1 (MQ=255)                                                        aaTCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCaccacc                                                                                       >  5:237023/1‑149 (MQ=255)                                                        aaTCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAACGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAGGAGCGCTGGGTGGTGATAGTGAAAACTGTATTCAGGTGGTGCCaccacc                                                                                       <  1:323891/149‑1 (MQ=255)                                                          tCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCgg                                                                                     <  6:97836/149‑1 (MQ=255)                                                             cgAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGTCCAAACCGGTTGTCGTTTCCCAGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGCTGCCACCACCGGTg                                                                                   >  4:296480/1‑148 (MQ=255)                                                                aCTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGaaaa                                                                               <  6:128120/149‑1 (MQ=255)                                                                  tCACCCCGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGc                                                                             >  4:196370/1‑149 (MQ=255)                                                                  tCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGc                                                                             >  7:175413/1‑149 (MQ=255)                                                                      ccAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCaga                                                                         <  2:266271/149‑1 (MQ=255)                                                                      ccAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCaga                                                                         <  6:236041/149‑1 (MQ=255)                                                                       cAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGaaaa                                                                                                               <  6:250392/110‑1 (MQ=255)                                                                       cAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGaaaa                                                                                                               >  5:250392/1‑110 (MQ=255)                                                                             gACGATGTCCGCAAAGCCGCGCAGGAAAGCGTTACTAGCGGCAAGGCCAAACGGGTTCGCGTTTCACAGGGT‑CCACAAGGAGCGCTGGGGGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGt                                                                  <  7:239514/149‑1 (MQ=255)                                                                             gACGATGACAGCAACGCCGCGCAGGCAATCGTTAATAGCTGCAAGGCCAATCGGGTTGTCGTTTCCCAGTGT‑CCACAACGAGCGCTGGGTGTTGCTAGGGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGGACCgg                                                                  <  7:239513/149‑2 (MQ=255)                                                                              aCGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGtt                                                                 <  5:257683/149‑1 (MQ=255)                                                                                      cGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTgg                                                         >  8:29815/1‑149 (MQ=255)                                                                                      cGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTg                                                          <  8:57000/148‑1 (MQ=255)                                                                                        cAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAg                                                                                                                                        >  1:102681/1‑68 (MQ=255)                                                                                        cAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAg                                                                                                                                        <  2:102681/68‑1 (MQ=255)                                                                                         aaaGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGa                                                      >  4:267681/1‑149 (MQ=255)                                                                                            gCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAg                                                   <  8:138699/149‑1 (MQ=255)                                                                                               gcgcAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCAt                                                <  6:237023/149‑1 (MQ=255)                                                                                                cgcAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATg                                               <  8:327517/149‑1 (MQ=255)                                                                                                   aGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTc                                            >  4:72607/1‑149 (MQ=255)                                                                                                     gAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCgg                                          >  5:253188/1‑149 (MQ=255)                                                                                                     gAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGGGAAAACTGTATTCAGGTGTTGCCACCACCGGTGAAAAGCAAGAGTACCGTTGGCGCTGGTGACCGCTTGGTCgg                                          >  7:255782/1‑149 (MQ=255)                                                                                                         tCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAACCTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGGCGgcgc                                       >  4:5712/1‑148 (MQ=255)                                                                                                         tCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCaga                                                                         >  3:181048/1‑114 (MQ=255)                                                                                                         tCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCaga                                                                         <  4:181048/114‑1 (MQ=255)                                                                                                           gTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGa                                     <  3:72607/148‑1 (MQ=255)                                                                                                           gTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGAGCTGGTGACAGCATGGTCGGCGCGAt                                    >  2:272948/1‑149 (MQ=255)                                                                                                              aaTAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGa                                  <  3:61305/148‑1 (MQ=255)                                                                                                              aaTAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGGGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGa                                  <  1:269499/148‑1 (MQ=255)                                                                                                               aTAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACTACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGaca                                <  3:250882/149‑1 (MQ=255)                                                                                                               aTAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGT‑CCACAAGGAGCGCTGGGGGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACTACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGaca                                <  1:46010/149‑1 (MQ=255)                                                                                                               aTAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGaca                                >  4:163263/1‑149 (MQ=255)                                                                                                               aTAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGaca                                <  8:175413/149‑1 (MQ=255)                                                                                                                  gCGGCAATCCCTAACGCGTTGTCGTTTCCCAGGGT‑CCACAAGGAGCGCTGGGTGTTGATATTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACAGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTg                             <  6:7167/149‑1 (MQ=255)                                                                                                                     gCAAGGCCAAACGGGTTGTCGTTTCCCTGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGaaa                          <  1:295540/149‑1 (MQ=255)                                                                                                                     gCAAGGCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCCCCGGTGATAAGCCAGAGTACCGTTGGCGCTGGTGACGGCATGGTCGGCGCGCTGACACTGaaa                          >  3:257820/1‑149 (MQ=255)                                                                                                                        aGGCCAAACGGGTTGTCGTTTCCCTGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTg                                                          >  8:223533/1‑114 (MQ=255)                                                                                                                        aGGCCAAACGGGTTGTCGTTTCCCTGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTg                                                          <  7:223533/114‑1 (MQ=255)                                                                                                                         ggCCAAACGGGTTGTCGTTTCCCAGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGGGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTgg                      <  1:23003/149‑1 (MQ=255)                                                                                                                                cGGGTTGTCGTTTCCCTGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGGGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGaaaa               <  8:273902/149‑1 (MQ=255)                                                                                                                                cGGGTTGTCGTTTCCCAGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGaaaa               >  1:178272/1‑149 (MQ=255)                                                                                                                                    ttGTCGTTTCCCAGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGgc                                         <  8:46956/119‑1 (MQ=255)                                                                                                                                       tCGTTTCCCAGGGT‑CCACAAGGAGAGCTTGGAGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCctc         >  8:77079/1‑148 (MQ=255)                                                                                                                                              ccTGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTg    >  8:5355/1‑146 (MQ=255)                                                                                                                                              ccTGGGT‑CCACAAGGAGCGCGGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGaa  <  3:196370/148‑1 (MQ=255)                                                                                                                                                tGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGaa  >  7:186345/1‑146 (MQ=255)                                                                                                                                                |                                                                                                                                                   AATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGT‑CCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAA  >  NC_000913/1806900‑1807188

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