Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,938,548	A→T	T361S (ACG→TCG)	fucK →	L‑fuculokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,938,548	0	A	T	98.2%	149.6 / NA	54	T361S (ACG→TCG)	fucK	L‑fuculokinase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (20/33);  total (21/33)
Fisher's exact test for biased strand distribution p-value = 3.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCACGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGAATCAGATTAAAGCCAATATGCTTGATATTCCGGTAAAAGTTCTCGACGACGCCGAAACGACCGTCGCAGGAGCTGCGCTGTTCGGTTGGTATGGCGTAGGGG  >  NC_000913/2938301‑2938791                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    gTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTcgc                                                                                                                                                                                                                                                  >  2:310375/1‑251 (MQ=255)    aaGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACTCTTAATACCTCGCgg                                                                                                                                                                                                                                                <  1:161251/250‑1 (MQ=255)        cAGTACGCCGGTTCCACCTGCGAACTGGCTAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGc                                                                                                                                                                                                                                            >  2:38389/1‑250 (MQ=255)          gTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGGGACGCTTAATACCTCACGGGGGcctt                                                                                                                                                                                                                                         >  2:455323/1‑248 (MQ=255)              gCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCt                                                                                                                                                                                                                                      >  1:637700/1‑250 (MQ=255)                cGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGAAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGAGACGCTTAATACCTCGCGGGGGCATTTCTATc                                                                                                                                                                                                                                   >  1:701220/1‑251 (MQ=255)                 ggtgcctccTGAGAACTGGATAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATcg                                                                                                                                                                                                                                  <  2:222509/244‑1 (MQ=255)                  gTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATcgc                                                                                                                                                                                                                                 <  1:639325/251‑1 (MQ=255)                           gCGAACTGGATAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTg                                                                                                                                                                                                                        <  1:329015/251‑1 (MQ=255)                            cGAACTGGAAAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGAATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGAAGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTAGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCACGCGGGGGCATTTCTATCGCGCGGCGCTgg                                                                                                                                                                                                                       >  2:790191/1‑251 (MQ=255)                             gAACTGGATAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGa                                                                                                                                                                                                                      <  2:137592/251‑1 (MQ=255)                                   gATAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATTCCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGtt                                                                                                                                                                                                                <  2:333476/251‑1 (MQ=255)                                   gATAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTG‑ACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTa                                                                                                                                                                                                               <  1:436157/251‑1 (MQ=255)                                        ccAGGCTGGGTTGTATAAACTAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACCCCCTGGAAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGGGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTg                                                                                                                                                                                                           <  2:252683/251‑1 (MQ=255)                                          aggcaggGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTgcg                                                                                                                                                                                                         >  1:798610/1‑251 (MQ=255)                                                ggTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAAtt                                                                                                                                                                                                    >  1:616610/1‑250 (MQ=255)                                                        aaCCCAGGTATGCAATGGCTTGGAGCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCa                                                                                                                                                                                            <  2:441888/250‑1 (MQ=255)                                                          cccAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATc                                                                                                                                                                                         <  2:805228/251‑1 (MQ=255)                                                            cAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGCTTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAACTGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGAGTAACTGCGCACTTCCAGCGCAATCt                                                                                                                                                                                        >  1:83596/1‑250 (MQ=255)                                                                 aTGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACCCCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGa                                                                                                                                                                                   <  2:147152/250‑1 (MQ=255)                                                                   gCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGc                                                                                                                                                                                <  2:682099/251‑1 (MQ=255)                                                                    cAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTg                                                                                                                                                                                                           <  1:107927/223‑1 (MQ=255)                                                                     aaTGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTTCGCAATTACAGCGCAATCTACAGATGCTg                                                                                                                                                                              >  2:224285/1‑251 (MQ=255)                                                                          ctggcATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGa                                                                                                                                                                            >  2:610113/1‑248 (MQ=255)                                                                            ggcATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGaaaaaa                                                                                                                                                                       <  1:345885/251‑1 (MQ=255)                                                                             gcATTCGGCGTGCTTGAATGGGTTAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTTGCTGGCAAGGAGTGACGCTTAATACCTCGCTGGGTCATTTCTATCGCGCGGCGCTTGAAGAGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAAt                                                                                                                                                                      <  2:92480/251‑1 (MQ=255)                                                                                        gCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACttt                                                                                                                                                            <  2:45277/250‑1 (MQ=255)                                                                                             aaTGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGcc                                                                                                                                                      <  1:290906/251‑1 (MQ=255)                                                                                                          aCTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGt                                                                                                                                         <  1:474312/251‑1 (MQ=255)                                                                                                                 tGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGCACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAg                                                                                                                                                                                                >  2:675840/1‑189 (MQ=255)                                                                                                                 tGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAg                                                                                                                                                                                                <  1:675840/189‑1 (MQ=255)                                                                                                                          gAAACGCCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCGATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGATGaa                                                                                                                          <  1:720016/250‑1 (MQ=255)                                                                                                                            aaCACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAACCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTc                                                                                                                       <  1:494326/251‑1 (MQ=255)                                                                                                                                   tGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGAAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACACTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAacac                                                                                                                >  1:68465/1‑251 (MQ=255)                                                                                                                                      cAAATGTTGATTGAAGAAGCTCGTCTGATCGCTCCTGGCGCGGATGGCGTAAAAATGCATTGTGATTTATGGTCGGTTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTTCAGATGCTGGAAAAAATCGGGCTCTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACAtt                                                                                                             <  2:637700/251‑1 (MQ=255)                                                                                                                                         aTGTTGATTTAAGAAGCTCGTCTGATCGCGCCTGGCGCGGAGGGCGTAAAAATGCAGTGTGATTTATTGTCTTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATtgtg                                                                                                          <  2:746257/251‑1 (MQ=255)                                                                                                                                             tgattgaAGAAGCTCGTCTGATCGAGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGAAt                                                                                                      >  1:353038/1‑251 (MQ=255)                                                                                                                                                                    gcgcCTGGAGCGGATGGCGTAAAAATGCAGTGTGATTTATTTTCGTGGCAGAACGGTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGAATCAGATTAAAGCCAATATGCTTGa                                                                               <  2:701220/251‑1 (MQ=255)                                                                                                                                                                       ccTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGAATCAGATTAAAGCCAATATGCTTGatat                                                                            <  1:551794/251‑1 (MQ=255)                                                                                                                                                                             gcgGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACGGCGCAATTACAGCGCAATCTACAGATGCTGTAAAAAATCGGGCTCTTTAAGTCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGAATCAGATTAAAGCCAATATGCTTGATATTCCGGt                                                                      <  1:743190/251‑1 (MQ=255)                                                                                                                                                                                   ggCGTAAAAATGAAGTGTGATTTATTGTCGTGTCAGAACGCTGTCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACCCATTGTGGAATCAGATTAAAGCCAATATGATTGATATTCAGGTAAAAGt                                                                >  2:806811/1‑251 (MQ=255)                                                                                                                                                                                     cGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGAATCAGATTAAAGCCAATATGCTTGTTATTCCGGTAAAAGttt                                                              >  2:573062/1‑250 (MQ=255)                                                                                                                                                                                     cGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGAATCAGATTAAAGCCAATATGCTTGATATTCCGGTAAAAGTtc                                                              <  1:275443/251‑1 (MQ=255)                                                                                                                                                                                        aaaaaTGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGAATCAGATTAAAGCCAATATGCTTGATATTCCGGTAAAAGTtctc                                                            <  1:224285/250‑1 (MQ=255)                                                                                                                                                                                                                gtgtCAGAACGCTGGCTGGCAAGGAGTGCCGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTCCAGCGCAATCTACCGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGAATCAGATTAAAGCCACTATGCTTGATATTCCGGTGAAAGTTCTCGACGACGCCGAAACGCACGTCGc                                     >  2:617850/1‑249 (MQ=255)                                                                                                                                                                                                                      gAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGAATCAGATTAAAGCCAATATGCTTGATATTCCGGTAAAAGTTCTCGACGACGCCGAAACGACCGTCGCAGGAGCTg                             <  2:45936/251‑1 (MQ=255)                                                                                                                                                                                                                          gctggctgGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGAATCAGATTAAAGCCAATATGCTTGATATTCCGGTAAAAGTTCTCGACGACGCCGAAACGACCGTCGCAGGAGCTGCGCt                         <  2:798610/251‑1 (MQ=255)                                                                                                                                                                                                                              gctggcAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGAATCAGATTAAAGCCAATATGCTTGATATTCCGGTAAAAGTTCTCGACGCAGCCGAAACGACCGTCGCAGGAGCTGCGCTGTTc                     >  1:516891/1‑251 (MQ=255)                                                                                                                                                                                                                              gctggcAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTAAAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGAATCAGATTAAAGCCAATATGCTTGATATTCCGGTAAAAGTTCTCGACGACGCCGAAACGACCGTCGCAGGAGCTGCGCTGTTc                     <  1:455323/251‑1 (MQ=255)                                                                                                                                                                                                                                   cAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCCACACATTGTGGACTCAGATTAAAGCCAATATGCTTGATATTCCGGTAAAAGTTCTCGACGACGCCGAAACGACCGTCGCAGGAGCTGCGCTGTTCGGTTg                >  2:79904/1‑251 (MQ=255)                                                                                                                                                                                                                                      ggAGTGACGCTTAATACCTCTCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCGGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGAATCAGATTAAAGCCAATATGCTTGATATTCCGGTAAAAGTTCTCGACGACGCCGAAACGACCGTCGCAGGAGCTGCGCTGTTCGGTTGGTa             <  2:408858/251‑1 (MQ=255)                                                                                                                                                                                                                                          gccgatctAATACCTCGCGGGGGCATTTCGATCGCGCTGCGCTGGAAGGGTTAACCGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGAATCAGATTAAAGCCAATATGCTTGATATTCCGGTAAAAGTTCTCGGCGACGCCGAAACGACCGTCGCAGGAGCTGCGCTGTTCGGTTGGTATgg          <  2:58988/243‑1 (MQ=255)                                                                                                                                                                                                                                                  aaTACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGACGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGACACACATTGTGGCATCAGTTTAAAGCCAATATGCTTGATATTCCGTTAAAAGTTCTCGACGTCGCCGAACCGACCGTCGCAGGAGCTGCGCTGTTCGGTAGGTCGTGCGTAgggg  >  2:730204/1‑250 (MQ=255)                                                                                                                                                                                                                                                  aaTACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGAATCAGATTAAAGCCAATATGCTTGATATTCCGGTAAAAGTTCTCGACGACGCCGCAACGACCGTCGCAGGAGCTGCGCTGTTCGGTTGGTATgg          >  1:58988/1‑242 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    GTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCACGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGAATCAGATTAAAGCCAATATGCTTGATATTCCGGTAAAAGTTCTCGACGACGCCGAAACGACCGTCGCAGGAGCTGCGCTGTTCGGTTGGTATGGCGTAGGGG  >  NC_000913/2938301‑2938791

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