Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,689,728	A→T	7.0%	F178L (TTT→TTA) ‡	bcsZ ←	endo‑1,4‑D‑glucanase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,689,728	0	A	T	7.0%	237.5 / 5.2	92	F178L (TTT→TTA) ‡	bcsZ	endo‑1,4‑D‑glucanase
Reads supporting (aligned to +/- strand):  ref base A (35/50);  new base T (5/2);  total (40/52)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.14e-01

TCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAAC  >  NC_000913/3689582‑3689873                                                                                                                                                   |                                                                                                                                                  gcTTGCTCATAGCGCACCCCGGCTGGCGAAAAACCGTTCGGCGCGCGTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGGCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTaaa                                                                                                                                                 <  1:103926/148‑1 (MQ=255)   tttCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAAc                                                                                                                                               <  1:537961/149‑1 (MQ=255)   tttCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAAc                                                                                                                                               >  2:364923/1‑149 (MQ=255)     tctcATAGCGCACCCAGTCTGGCGAAAAACATTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGCTCCACGGCGCGCCAAAGCGGTTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGc                                                                                                                                             >  1:41583/1‑149 (MQ=255)      ctcATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGGTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGcc                                                                                                                                            <  1:502590/149‑1 (MQ=255)       tcATAGCGCACCCAGTCTGGCGAAAACCCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCa                                                                                                                                           <  1:723221/149‑1 (MQ=255)        cATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGTGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGCGTTTTAACCCCAg                                                                                                                                          <  2:959616/149‑1 (MQ=255)         aTAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGc                                                                                                                                         >  1:707140/1‑149 (MQ=255)            gcgcACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTTGTTTCGAGCATCGTGGTCCACGGAGCGCCAAAGCGGGTGAAATACTGCCCCATCGTCGCAGGCAGGTAGCTGGGATTAAAACGcaacctaa                                                                                                                                      >  1:1028155/1‑142 (MQ=255)            gcgcACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGCGTTTTAACCCCAGCtac                                                                                                                                      <  1:873264/149‑3 (MQ=255)               cACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGGGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTAt                                                                                                                                   >  1:209598/1‑149 (MQ=255)               cACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTAt                                                                                                                                   >  1:658360/1‑149 (MQ=255)                  ccAGTCTGGCGAAAAACCTTTCGGGCCGGTTTCCAGCAATAATCGTTGATGGGTTGCGCGCAGCGTGGTTCTCGGCGCGCCAAAGCGGGTGAAATACTTCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCt                                                                                                                                <  2:947621/149‑1 (MQ=255)                   cAGTCTGTCGAACAACCTTTCGGCGAGGTTTCCAGAAATAAACGTTGATTGGGTTCGCGCAGCGTGGTCAACGGCGCGAAAAAGCGGTTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTc                                                                                                                               <  1:947446/149‑1 (MQ=255)                    aggcTGGCGAAAAACCTTTCGGGGCGGATGCCAGCAATAACCGTGGATTGGTCTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCa                                                                                                                              <  1:774229/146‑1 (MQ=255)                    aGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCa                                                                                                                              <  1:781468/149‑1 (MQ=255)                      tCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGc                                                                                                                            <  1:142933/149‑1 (MQ=255)                      gcTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGc                                                                                                                            <  1:1001175/148‑1 (MQ=255)                          gCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCaaac                                                                                                                        >  1:516255/1‑148 (MQ=255)                          gCGAAAAACCTTTCGGGCCGGTTTCCAGCAATAAACGTTGATTGGTATCGCGCAGCGTGCTCCACGGCGCGCCAAAGCGGGAGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCACCTGTTATACTCAGCaaaa                                                                                                                        >  2:321583/1‑149 (MQ=255)                           cGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAAc                                                                                                                       <  1:1019035/149‑1 (MQ=255)                            gAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAAcc                                                                                                                      >  2:1010667/1‑149 (MQ=255)                                  ccTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACttt                                                                                                                <  2:812664/149‑1 (MQ=255)                                     ttCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAACCGCCAGCTGTTATCCTCCGCAACACCCCCTTTGcc                                                                                                             >  2:945088/1‑149 (MQ=255)                                       cGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCg                                                                                                           <  1:403401/149‑1 (MQ=255)                                         gtgCGGCTTCCACCAATAAACGTTGATTGGTTGCGCGCAGCGTGGGCCACGGCTCTCCAAAGCGGGTGAAATACTGCGCCAGCGTCTGCGGCAGGTAGCTGGGGTTATAACGCCAGCTGTTATCCTCAGAAAAACCCACTTTGCCCGGt                                                                                                         <  2:903946/147‑1 (MQ=255)                                         gggCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGt                                                                                                         <  1:748247/149‑1 (MQ=255)                                          ggCGGTTTCCAGCAATAAACGTGGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTTGGGTTAACACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTa                                                                                                        >  1:736749/1‑149 (MQ=255)                                           gCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTaa                                                                                                       >  1:404258/1‑149 (MQ=255)                                            cGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGTAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTaac                                                                                                      <  2:645323/149‑1 (MQ=255)                                            cGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTaac                                                                                                      <  1:474711/149‑1 (MQ=255)                                            cGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGACAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTaac                                                                                                      <  1:474698/149‑1 (MQ=255)                                               tttCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTaacaac                                                                                                   <  2:441178/149‑1 (MQ=255)                                               tttCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGACGGCGGCAGGTAGCTGGGGTTAATCCGCCAGCTGTTATCCTCAGCCAACCCCACTTGGCCCGGTaacacc                                                                                                   >  1:740934/1‑147 (MQ=255)                                                 tCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCTCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAACCGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTTACAACAt                                                                                                 >  2:511152/1‑149 (MQ=255)                                                  ccAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCAGTTATCCTCAGCAAAACCCACTTTGCCCGGTAacaccaag                                                                                                >  1:420365/1‑147 (MQ=255)                                                       aaTAAACGTTGATTGGTTTCGCGCAGCGTTGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAAcc                                                                                           <  2:404258/149‑1 (MQ=255)                                                       aaTAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGGCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAAcc                                                                                           <  1:220054/149‑1 (MQ=255)                                                         tAAACGTGGATTGGTGTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGGGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAAccca                                                                                         <  2:149703/149‑1 (MQ=255)                                                           accGTTGAATGGATTCGCGCAGCGTGGTCCCCGGCGCGACAACGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGGTAAAAAGCCAGCTGTGAGCCTCCGCAACCCCCACTTTGGCCGGTCACTCCAAGGCAccccgc                                                                                       >  2:174783/3‑149 (MQ=255)                                                            aCGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAAccccgcc                                                                                      >  1:778395/1‑149 (MQ=255)                                                               ttGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAAccccgcccag                                                                                   >  1:543554/1‑149 (MQ=255)                                                               ttGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAAcccagcccag                                                                                   <  2:225998/149‑1 (MQ=255)                                                                tGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAATCCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAgg                                                                                  <  1:169999/149‑1 (MQ=255)                                                                  aTTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCTGCTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTATAAACGACAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGAAACAACATGGAACCCAGCCGAGGCa                                                                                >  2:570318/1‑149 (MQ=255)                                                                       tttCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAACCCCACTTTGCCCGGTAACCACATGGAACCCAGCCCCGACCCCGTc                                                                           >  2:857728/1‑149 (MQ=255)                                                                       tttCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTc                                                                           <  1:43978/149‑1 (MQ=255)                                                                       tttCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTc                                                                           <  1:846052/149‑1 (MQ=255)                                                                         tCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCATAACCCACTTTTCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTcac                                                                         <  2:272331/149‑1 (MQ=255)                                                                         tCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTcac                                                                         <  1:986970/149‑1 (MQ=255)                                                                          agcgcAGCGTGGTCCACGGCGCGCCAAAGCTGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGGTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTcacc                                                                        <  2:207206/148‑1 (MQ=255)                                                                           gcgcAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAACCGCCAGCTGTTATCCTCAGCAAACCCCACTTTGCCCGGTAACAACATGGACCCCAGCCCAGGCACCGTcaccc                                                                       >  2:33999/1‑148 (MQ=255)                                                                                gCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTtcc                                                                  <  2:707140/149‑1 (MQ=255)                                                                                  gTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTtcctc                                                                <  1:460507/149‑1 (MQ=255)                                                                                     gTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATACTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCcg                                                             >  1:556101/1‑149 (MQ=255)                                                                                               cgcCAAATAGTGTGAAATTCTGCGCCAGAGTCTGCGGCAGGTATCTGGTGTTAAAACGCCATCTTTTATCCTCAGCAAAACCCACTTTGCCCGGTAAAAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCAGCGCTATACGt                                                   <  1:23523/149‑1 (MQ=255)                                                                                                 ccAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCCGCAAACCCCCCTGTGCCCGGTAACAACAGGGACCCCACCCCAGGCACAGTCACCACTTCCTCCCGCGCGATActgtt                                                 >  2:160510/1‑145 (MQ=255)                                                                                                 ccAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGttt                                                 >  2:212577/1‑149 (MQ=255)                                                                                                  cAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCCCTTTGCCCGGTAACAACATGGAACCCAGCCCCGGCACCGTCCCCACTTCCTCCCGCGCGATACGtttt                                                >  1:633808/1‑149 (MQ=255)                                                                                                       cGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCaa                                           <  1:229734/149‑1 (MQ=255)                                                                                                        gggTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGATGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAAc                                          <  2:761364/149‑1 (MQ=255)                                                                                                         ggTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAACCCCACTTTGCCCGGTAACAACATGGAACCCAGCCCCGGCACCGTCCCCCCTTCCTCCCGCGCGCTACGTTTTAGCCAcg                                         >  1:439462/1‑149 (MQ=255)                                                                                                         ggTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCTACTTGGCCCGGTAACAACATGGAACCCAGCCCCGGCGCCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAAcg                                         >  2:106833/1‑149 (MQ=255)                                                                                                         ggTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAAcg                                         <  1:877381/149‑1 (MQ=255)                                                                                                            gAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGCACCCAGCCCAGGCACCGTCTCCACTTCCTCCCGCGCGAGACGTTTTAGCAACgcgc                                      >  1:414527/1‑149 (MQ=255)                                                                                                               aTACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGc                                   <  2:664284/149‑1 (MQ=255)                                                                                                                tACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAACTCCCACTTTGCCCGGTACCAACATGGATGCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGGTTTAGCAACGCGCTGcc                                  >  1:52037/1‑149 (MQ=255)                                                                                                                   tGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACGTCCTCCCGCGCGATACGTTTTAGTAACGCGCTGCCGAt                               >  1:984908/1‑149 (MQ=255)                                                                                                                      gcCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACACCGCCCAGGCACCGTCCCCACTGCCTCCCGCGCGGTACCTGTTAGCGACGCGCTGCCGaaggc                            >  2:332073/1‑145 (MQ=255)                                                                                                                       ccATCGTCGTCGGCAGGTAGCTGGGGTAAAAACGCCAGCTGATATCCTCAGCAAAACCCACTTTGCCCTGTAACATCATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCg                           <  1:672900/149‑1 (MQ=255)                                                                                                                       ccAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCg                           <  1:206116/149‑1 (MQ=255)                                                                                                                       ccAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCg                           <  2:658360/149‑1 (MQ=255)                                                                                                                       ccAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTAACCTCAGCAAAACCCACTTTCCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCg                           <  1:590749/149‑1 (MQ=255)                                                                                                                       ccAGAGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCTGCAAAACCCACTTTGCCCGGTTACAACATGGAACCCAGCCGAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACTCGCTGCCGATGTCg                           <  1:321583/149‑1 (MQ=255)                                                                                                                         aGCGTCGGCGGCAGGTAGCGGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGACACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGt                         <  1:857728/149‑1 (MQ=255)                                                                                                                               ggcggcAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCCCCGGCCCCCCTTCCTCCCGCGCGATACGGTGTAGCAACGCGCTGCCAATGTCGGGAGAACg                   >  1:824034/1‑149 (MQ=255)                                                                                                                                 cggcAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACgct                 <  2:367010/149‑1 (MQ=255)                                                                                                                                    ataacaGCTGGCGTTTTAACCCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACgctgct              >  2:622844/6‑149 (MQ=255)                                                                                                                                      ggTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAACCCCAATTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGctct            <  2:452828/149‑1 (MQ=255)                                                                                                                                      aacaGCTGGCGTTTTAACCCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCCCTTCCTCCCGCGCGATACGTTTTAGCCACGCGCTGCCGATGTCGGTATAACGCTGctct            >  2:950072/4‑149 (MQ=255)                                                                                                                                        tAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCttt          <  1:348460/149‑1 (MQ=255)                                                                                                                                        tAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCttt          <  1:815595/149‑1 (MQ=255)                                                                                                                                        caGCTGGCGTTTTAACCCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACAGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGTTGTCGGTATAACGCTGCTCttt          >  2:146476/2‑149 (MQ=255)                                                                                                                                         aGCTTGGGTTAAAACGCCAGCTGTTATCCTCAGCAACAGCCACTTTGCCCGGAAACATCATGGAACCCAGCCCAGGCACCGTCACCACTTCGTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTc         <  2:210982/149‑1 (MQ=255)                                                                                                                                         aGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTc         <  1:831775/149‑1 (MQ=255)                                                                                                                                         aGCTGGCGTTTTAACCCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGACCCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTc         >  1:883011/1‑149 (MQ=255)                                                                                                                                          gCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAACCCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACAGTCACCCCTTCCTCCCGCGCGATACTTTTTAGCAACGCGCTGCCGATGGCGGTATAGCGCTGCTCTTTcc        >  1:129544/1‑149 (MQ=255)                                                                                                                                            gggggCTAAAACGCCAGCTGTTACCCTCATCAAAACCCACTGTGCCCGGTAACAACATGGAACCCAGCCCAGGGACCGTCACCACTTCCTCCCGCGCGCTACGTGTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCcac      <  2:825047/148‑1 (MQ=255)                                                                                                                                             ggggTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCcaca     >  2:297908/1‑149 (MQ=255)                                                                                                                                               ggTTAAAACGCCAGCTGTTATCCTCAGCAAACCCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACaaa   <  2:152090/149‑1 (MQ=255)                                                                                                                                               ggTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCGGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACaaa   <  2:660700/149‑1 (MQ=255)                                                                                                                                                gTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACGGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAAc  >  2:88938/1‑149 (MQ=255)                                                                                                                                                gTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCCCCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAAc  >  1:721858/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                  TCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAAC  >  NC_000913/3689582‑3689873

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