Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,163,979	G→A	5.2%	L240L (CTG→CTA)	thiK →	thiamine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,163,979	0	G	A	5.2%	244.6 / 2.4	77	L240L (CTG→CTA)	thiK	thiamine kinase
Reads supporting (aligned to +/- strand):  ref base G (32/41);  new base A (2/2);  total (34/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01

ATGGTGATATCGCGCTGGAACTGGCGGCGGTGTGGGTGGAAAATACTGAACAGCACCGGCAATTGGTCAATGACTATGCCACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTGGCGGCAGCTATTAATAAAACAATAAGGAGAGCAGTGTGGGTCCAGTAATGTTGGATGTCGAAGG  >  NC_000913/1163834‑1164123                                                                                                                                                  |                                                                                                                                                 aTGGTGATATCGCGCTGGAACTGGCGGCGGTGTTGGTGGAAAATACTGAACAGCACCGGCAATTGGTCAATGACTATGCCACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATg                                                                                                                                               <  2:257268/149‑1 (MQ=255)  tGGTGATATCGCGCTGGAACTGGCGGCGGTGTGGGTGGAAAATACTGAACAGCACCGGCAATTGGGCAAGGACTATGCCACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGc                                                                                                                                              >  1:352680/1‑149 (MQ=255)            gcgcTGGAACTGGCGGCGGTGTGGGTGGAAAATACTGAACAGCACCGGCAATTGGTCAATGACTATGCCACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAgg                                                                                                                                    <  2:426057/149‑1 (MQ=255)            gcgcTGGAACTGGCGGCGGTGTGGGTGGAAAATACTGAACAGCACCGGCAATTGGTCAATGACTATGCCACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAg                                                                                                                                     <  5:22612/148‑1 (MQ=255)             cgcTGGAACTGGCGGCTGTGTGGGTGGAAAATACTGAACAGCACCGGCTATTGGTCAATGACTATGCCACTCTCTCTAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAggg                                                                                                                                   <  3:252608/149‑1 (MQ=255)              ggtGGAACGGGCGGCGGGGTGGGGGGAAAATACTGAACAGCCCCGGCAATTGGTCAATGACTATGCCACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGGTCAAAGCAGggt                                                                                                                                  <  6:348904/147‑1 (MQ=255)                tGGAACTGGCGGCGGTGTGGGTGGAAAATACTGAACAGCACCGGCAATTGGTCAATGACTATGCCACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGggtgg                                                                                                                                >  7:158587/1‑149 (MQ=255)                     cTGGCGGCGGTGTGGGTGGAAAATACTGAACAGCACCGGCAATTGGTCAATGACTATGCCACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGa                                                                                                                           <  8:312140/149‑1 (MQ=255)                     cTGGCGGCGGTGTGGGTGGAAAATACTGAACAGCACCGGCAATTGGTCAATGACTATGCCACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGa                                                                                                                           <  5:293210/149‑1 (MQ=255)                        gcggcgGTGTGGGTGGAAAATACTGAACAGCACCGGCAATTGGTCAATGACTATGCCACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTa                                                                                                                        <  8:158587/149‑1 (MQ=255)                            cggTGTGGGTGGAAAATACTGAACAGCACCGGCAATTGGTCAATGACTATGCCACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGc                                                                                                                    >  1:246681/1‑149 (MQ=255)                               ggtgggtggAAAATACTGAACAGCACCGGCAATTGGTCACTGACTATGCCACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTgg                                                                                                                 <  6:384135/148‑1 (MQ=255)                                     tggAAAATACTGAACAGCACCGGCAATTGGTCAATGACTATGCCACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACa                                                                                                            >  8:210184/1‑148 (MQ=255)                                       ggaaaTACTGAACAGCACCGGCAATTGGTCAATGACTATGCCACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAAc                                                                                                         >  1:427871/3‑149 (MQ=255)                                           aTACTGAACAGCACCGGCAATTGGTCAATGACTATGCCACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGc                                                                                                     <  3:278506/149‑1 (MQ=255)                                              cttaaCAGCACCGGCAATTGGTCAATGACTATGCCACTCGCTCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCTGGCGATTCTTTCCCTGGCTGCTGTTGCTCAAAGCATGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGAt                                                                                                  <  6:109302/146‑1 (MQ=255)                                               tGAACAGCACCGGCAATTGGTCAATGACTATGCCACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGa                                                                                                   <  2:8402/147‑1 (MQ=255)                                               tGAACAGCACCGGCAATTGGTCAATGACTATGCCACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGAt                                                                                                  >  8:262436/1‑148 (MQ=255)                                               tGAACAGCACCGGCAATTGGTCAATGACTATGCCACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGAt                                                                                                  >  5:420059/1‑148 (MQ=255)                                                 aaCAGCACCGGCAATTGGTCAATGACTATGCCACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATc                                                                                                 <  7:262436/147‑1 (MQ=255)                                                         cGGCAATTGGTCAATGACTATGCCACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGa                                                                                                                                                 >  6:329113/1‑91 (MQ=255)                                                         cGGCAATTGGTCAATGACTATGCCACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGa                                                                                                                                                 <  5:329113/91‑1 (MQ=255)                                                          ggCAATTGGTCAATGACTATGCCACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCTCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGa                                                                                                              <  2:390485/125‑1 (MQ=255)                                                          ggCAATTGGTCAATGACTATGCCACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCTCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGa                                                                                                              >  1:390485/1‑125 (MQ=255)                                                          ggCAATTGGTCAATGACTATGCCACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAAt                                                                                           <  2:62281/144‑1 (MQ=255)                                                          ggCAATTGGTCAATGACTATGCCACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAAt                                                                                           >  1:62281/1‑144 (MQ=255)                                                           gCAATTGGTCAATGACTATGCCACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCa                                                                                     <  1:243124/149‑1 (MQ=255)                                                              aTTGGTCAATGACTATGCCACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTAATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGc                                                                                  >  6:448798/1‑149 (MQ=255)                                                              aTTGGTCAATGACTATGCCACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTAATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGc                                                                                  >  4:418485/1‑149 (MQ=255)                                                                     aaTGACTATGCCACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGa                                                                           >  1:296695/1‑149 (MQ=255)                                                                        gACTATGCCACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGa                                                                        >  1:331808/1‑149 (MQ=255)                                                                         aCTATGCCACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGAGCAACAATTTATCAGGCTGGCCGATGa                                                                        <  8:287393/148‑1 (MQ=255)                                                                           tATGCCACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATcaac                                                                                              >  7:229386/1‑124 (MQ=255)                                                                           tATGCCACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATcaac                                                                                              <  8:229386/124‑1 (MQ=255)                                                                             tGCCACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGAGGACACCt                                                                   >  2:116802/1‑149 (MQ=255)                                                                               ccACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTgg                                                                 >  7:344469/1‑149 (MQ=255)                                                                                cACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTggc                                                                <  2:246681/149‑1 (MQ=255)                                                                                 aCTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGacac                                                                     >  3:422265/1‑143 (MQ=255)                                                                                 aCTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGacac                                                                     <  4:422265/143‑1 (MQ=255)                                                                                         gAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTGGCGGCAGCTAt                                                       <  2:427871/149‑1 (MQ=255)                                                                                         gAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTGGCGGCAGCTAt                                                       <  7:209051/149‑1 (MQ=255)                                                                                         gAAGATTAATCCGGCGCAATTAAGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAAAATATCAGCCTGGCCGATGACACCTGGCGGCAGCTAt                                                       >  8:165508/1‑149 (MQ=255)                                                                                          aaGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTGGCGGCAGCTAtt                                                      <  8:344469/149‑1 (MQ=255)                                                                                           aGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTGTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTTGCCGATGACACCTGGCGGCAGCTATta                                                     <  6:424570/149‑1 (MQ=255)                                                                                             aTTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTGGCGGCAGCTATtaa                                                    >  3:370671/1‑148 (MQ=255)                                                                                                tATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTGGCGGCAGCTATtaataa                                                 <  1:116802/148‑1 (MQ=255)                                                                                                  tCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTGGCGGCAGCTATTAATaaaa                                               >  7:74399/1‑148 (MQ=255)                                                                                                  tCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTAATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTGGCGGCAGCTATTAATaaaa                                               <  3:418485/148‑1 (MQ=255)                                                                                                  tCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTAATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTGGCGGCAGCTATTAATaaaa                                               <  5:448798/148‑1 (MQ=255)                                                                                                     ggCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTGGCGGCAGCTATTAATAAAACAAt                                           <  6:420059/149‑1 (MQ=255)                                                                                                        gcAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGTTGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTGGCGGCAGCTATTAATAAAACAAt                                           >  3:346403/1‑146 (MQ=255)                                                                                                        gcAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTGGCGGCAGCTATTAATAAAACAATAAg                                        <  4:247280/149‑1 (MQ=255)                                                                                                          aaTTCTGGCGTCCGGTCAGGCGTTGGTTTCCCTGGCTGCTGAGGCTCAAAGCATGGTGGTATGAGTACCGCTGGCGACAAACCGGCGATCAACTATTTATCTGGCTGGCCGATGACACCTGGCGGCAGCTATTAATAAAACAATAAGga                                      <  2:352680/149‑1 (MQ=255)                                                                                                          aaTTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGTTGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTGGCGGCAGCTATTAATAAAACAATaa                                         <  4:346403/146‑1 (MQ=255)                                                                                                           aTTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTGGCGGCAGCTATTAATAAAACAATAAGgag                                     <  4:398532/149‑1 (MQ=255)                                                                                                            ttATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTGGCGGCAGCTATTAATAAAACAATAAGga                                      >  7:65905/1‑147 (MQ=255)                                                                                                               tGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTGGCGGCAGCTATTAATAAAACAATAAGGagag                                   >  5:416337/1‑147 (MQ=255)                                                                                                                   gtcaggtcagGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTGGCGGCAGCTATTAATAAAACAATAAGGAGAGCAgtgt                             <  8:65905/149‑1 (MQ=255)                                                                                                                      aggtcaggCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTGGCGGCAGCTATTAATAAAACAATAAGGAGAGCAGTGTgg                           <  8:74399/148‑1 (MQ=255)                                                                                                                        gtcaggCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTGGCGGCAGCTATtaata                                                  <  2:59995/123‑1 (MQ=255)                                                                                                                        gtcaggCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTGGCGGCAGCTATtaata                                                  >  1:59995/1‑123 (MQ=255)                                                                                                                           aggCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAAGTTATCAGGCTGGCCGATGACACCTGGCGGCAGCTAGTAATAAAACAATAAGGAGAGCAGTGTGGGTCCAg                     >  4:333728/1‑149 (MQ=255)                                                                                                                              cGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGTTCAACAATTTATCAGGCTGGCCGATGACACCTGGCGGCAGCTATTAATAAAACAATAAGGAGAGCAGTGTGGGTCCAGTaa                  <  7:210184/149‑1 (MQ=255)                                                                                                                              cGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTGGCGGCAGCTATta                                                     >  7:335097/1‑114 (MQ=255)                                                                                                                              cGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTGGCGGCAGCTATta                                                     <  8:335097/114‑1 (MQ=255)                                                                                                                                aTGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTGGCGGCAGCTATTAATAAAACAATAAGGAGAGCAGTGTGGGTCCAGTAATg                >  3:86224/1‑149 (MQ=255)                                                                                                                                   gTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTGGCGGCAGCTATTAATAAAACAATAAGGAGAGCAGTGTGGGTCCAGTAATGtt              <  4:370671/148‑1 (MQ=255)                                                                                                                                    tttCCCTGGCTGCTGATGCTCTAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAAGCGTCGATCAAGAATTTATCAGGCTGGCCGATGACACGTGGCGGCAGCTATTAATAAAACAATAAGTAGAGCAGTGTGGTTCCAGTAATGTTg             <  7:74094/148‑1 (MQ=255)                                                                                                                                    tttCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTGGCGGCAGCTATTAATAAAACAATAAGGAGAGCAGTCTGGGTCCAGTAATGTTg             >  8:188044/1‑148 (MQ=255)                                                                                                                                    tttCCCTGGCTGCTGATGCTCAAAGCAGCGTCGTTTGAGAACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTGGCGGCAGCTATTAATAAAACAATAAGGAGAGCAGTGTGGGTCAAGTAATGTTg             >  8:74094/1‑148 (MQ=255)                                                                                                                                        ccTCGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTGGCGGCAGCTATTAATAAAACAATAAGGAGAGCAGTGTTGGTCCAGTAATGTTGGATg         <  8:92478/148‑1 (MQ=255)                                                                                                                                         cTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTGGCGGCAGCTATTAATAAAACAATAAGGAGAGCAGTGTGGGTCCAGTAATGTTgg            <  3:318881/144‑1 (MQ=255)                                                                                                                                         cTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTGGCGGCAGCTATTAATAAAACAATAAGGAGAGCAGTGTGGGTCCAGTAATGTTgg            >  4:318881/1‑144 (MQ=255)                                                                                                                                         cTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTGGCGGCAGCTATTAATAAAACAATAAGGAGAGCAGTGTGGGTCCAGTAATGTTGGATGTc       >  3:208769/1‑149 (MQ=255)                                                                                                                                         cTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTGGCGGCAGCTATTAATAAAACAATAAGGAGAGCAGTGTGGGTCCAGTAATGTTGGATGTc       >  3:202160/1‑149 (MQ=255)                                                                                                                                           ggCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTGGCGGCAGCTATTAATAAAACAATAAGGAGAGCAGTGTGGGTCCAGTAATGTTGGATGTCg      >  6:310780/1‑148 (MQ=255)                                                                                                                                            gctgctGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTGGCGGCAGCTATTAATAAAACAATAAGGAGAGCAGTGTGGGTCCAGTAATGTTGGATGTCGaa    <  2:296695/149‑1 (MQ=255)                                                                                                                                              tgctgATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTGGCGGCAGCTATTAATAAAACAATAAGGAGAGCAGTGTGGGTCCAGTAATGTTGGATGTCGAAgg  >  4:203490/1‑149 (MQ=255)                                                                                                                                                 tgATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATc                                                                                                 <  6:89611/51‑1 (MQ=255)                                                                                                                                                 tgATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATc                                                                                                 >  5:89611/1‑51 (MQ=255)                                                                                                                                                 tgATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTGGCGGCAGCTATTAATAAAACAATAAGGAGAGCAgtgt                             <  7:331647/119‑1 (MQ=255)                                                                                                                                                 tgATGCTCAAAGCAGGGTGGTTTGAGAACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTGGCGGCAGCTATTAATAAAACAATAAGGAGAGCAgtgt                             >  8:331647/1‑119 (MQ=255)                                                                                                                                                  |                                                                                                                                                 ATGGTGATATCGCGCTGGAACTGGCGGCGGTGTGGGTGGAAAATACTGAACAGCACCGGCAATTGGTCAATGACTATGCCACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTGGCGGCAGCTATTAATAAAACAATAAGGAGAGCAGTGTGGGTCCAGTAATGTTGGATGTCGAAGG  >  NC_000913/1163834‑1164123

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