Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,809,376	C→T	N279N (AAC→AAT)	waaA →	KDO transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,809,376	0	C	T	93.1%	184.4 / NA	59	N279N (AAC→AAT)	waaA	KDO transferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (32/22);  minor base G (0/4);  total (32/27)
Fisher's exact test for biased strand distribution p-value = 3.52e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTACTGTATGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGG  >  NC_000913/3809139‑3809622                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                       gttaaCCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGGCATCCAGAACGCGTCCCGGATGCGATTAAGCTTGTCCGCCAgg                                                                                                                                                                                                                                           <  2:70360/247‑1 (MQ=255) ggtaaCCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAgg                                                                                                                                                                                                                                           <  1:122616/247‑1 (MQ=255)   taaCCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCt                                                                                                                                                                                                                                         >  2:331572/3‑251 (MQ=255)      ccGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGACGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTgg                                                                                                                                                                                                                                       >  1:256107/1‑250 (MQ=255)           aCCGGTAGCCTGAAATTCGATATTTCTGGAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGTTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTaat                                                                                                                                                                                                                                 <  2:243210/251‑2 (MQ=255)           aCCGGTAGCCTGAAATTCGATATTGCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAAGCTTGTCCGCCAGGCTGGACTaa                                                                                                                                                                                                                                  <  1:161967/250‑1 (MQ=255)             cGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCt                                                                                                                                                                                                                               >  1:150687/1‑251 (MQ=255)                     tGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATcac                                                                                                                                                                                                                        >  1:268791/1‑250 (MQ=255)                      gAAATTCGATATTTCTGTAACGCCGCCGTTGGCTGCTACAGCCGTGACGCTGCGAAGCCCGTGGGCACCACACCGCCCGGTATGGATTGCCGCCAGCACTCACGAAGGCGCAGAGAGTGTGGTGATCGCCGCACATCCGGCATTGTTACAGCACTTCCCGAGTTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGCTGCGATAAGTCTTGCCCGCCAGGCTGGACTAAGCTATATAacac                                                                                                                                                                                                                      >  2:71778/1‑251 (MQ=255)                             gATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCAGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGACGGCGAAGAGCGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGCACGCTTCCCGGATGCGCTTACTCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTc                                                                                                                                                                                                               >  2:383994/1‑251 (MQ=255)                                       tAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTc                                                                                                                                                                                                     >  1:144006/1‑251 (MQ=255)                                                aGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAg                                                                                                                                                                                                             <  2:302735/234‑1 (MQ=255)                                                aGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAg                                                                                                                                                                                                             >  1:302735/1‑234 (MQ=255)                                                        gctAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCa                                                                                                                                                                                    <  1:416733/251‑1 (MQ=255)                                                           aaaGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAgg                                                                                                                                                                                  <  2:416088/250‑1 (MQ=255)                                                                    aCGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATGGTTACAGCAATTCGCTAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAAGCTTGTCCGCCATGCTGGACTAAGCTATATCACACGCTCTTCAGGTGAAGTGCCCTCCACCAGCACGCAGGTTGTGGtggg                                                                                                                                                                        <  2:429876/251‑4 (MQ=255)                                                                           gccgccACTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGAGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGAGGGCGATACg                                                                                                                                                                 >  2:172815/1‑251 (MQ=255)                                                                              gccAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGGGATACGATg                                                                                                                                                              >  2:181337/1‑251 (MQ=255)                                                                                      gCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGtt                                                                                                                                                      <  2:268791/251‑1 (MQ=255)                                                                                                    ggTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTacg                                                                                                                                        >  1:179341/1‑249 (MQ=255)                                                                                                     gTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTacgg                                                                                                                                       <  1:305821/251‑4 (MQ=255)                                                                                                      tATGGATTGCCACCAGCACTCACGAAGGCGCAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCGCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGCTACGATGGGCGGGTTGATGTTACTGTACGGc                                                                                                                                      >  2:80493/1‑251 (MQ=255)                                                                                                       aTGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCg                                                                                                                                                          <  1:323928/230‑1 (MQ=255)                                                                                                       aTGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCg                                                                                                                                                          >  2:323928/1‑230 (MQ=255)                                                                                                         ggATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCAtt                                                                                                                                   <  1:424813/251‑1 (MQ=255)                                                                                                             tGCCACCAGCAATCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGGGGGCGATACGATGGGCGTGTTGATGTTACTGTACGGCCTTGCCg                                                                                                                               >  2:393255/1‑251 (MQ=255)                                                                                                               ccaccaGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGAt                                                                                                                             >  1:159359/1‑251 (MQ=255)                                                                                                                     gCACTCACGAAGGCGAAGATAGTGTGGTGATCGCCGCACATCAGGCATTGCTACAGCAACTCCCGAATTTATTGCTCATCCTGGTACCCCGGCAACCAGAACGCTGCCTGGATGCGATTAAACGTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCTCCAGCACGCAGGTGGTGGTGGGCGATACGATGGGCGAGTTGATATTACTGTACGTCATTGCCGATCTCGcc                                                                                                                       <  2:152495/251‑1 (MQ=255)                                                                                                                      cACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCt                                                                                                                      >  1:314245/1‑251 (MQ=255)                                                                                                                           aCGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTTGTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGCACGGCATTGCCGATCTCGCCTttgtt                                                                                                                 <  1:242195/251‑1 (MQ=255)                                                                                                                                       agagTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTAGTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTGGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCGTTGTTGGCGGTTCACTg                                                                                                     <  2:179341/251‑1 (MQ=255)                                                                                                                                           tgtgGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTg                                                                                                 >  2:125539/1‑251 (MQ=255)                                                                                                                                               gtgATCGCCGCACATCAGGCATTGTGACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCTCCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACg                                                                                             <  2:256107/251‑1 (MQ=255)                                                                                                                                                 gATCGCCGCACATCAGGCATTGTTCCAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACgtg                                                                                           >  2:296621/1‑251 (MQ=255)                                                                                                                                                      ccgcACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTggg                                                                                      >  1:258740/1‑251 (MQ=255)                                                                                                                                                       cgcACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGc                                                                                     >  2:221749/1‑251 (MQ=255)                                                                                                                                                       cgcACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTCGTAACCCGTCATCCAGAACGCTTCCCGGCTGCGATTCATCTTGTCAGCCAGGCTGGACTAAGCTATATCACCAGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGGTGAACGTGGTGGGc                                                                                     >  1:415969/1‑251 (MQ=255)                                                                                                                                                                  cATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGc                                                                                     >  2:131547/1‑240 (MQ=255)                                                                                                                                                                  cATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGc                                                                                     <  1:131547/240‑1 (MQ=255)                                                                                                                                                                       ttACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGc                                                                     >  1:267692/1‑251 (MQ=255)                                                                                                                                                                         aCAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTg                                                                   >  1:259727/1‑251 (MQ=255)                                                                                                                                                                             cAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGc                                                               >  2:252935/1‑251 (MQ=255)                                                                                                                                                                                 tCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGcaca                                                            <  2:159359/250‑1 (MQ=255)                                                                                                                                                                                  cccGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCATCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCTCTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCacac                                                           <  2:144006/250‑1 (MQ=255)                                                                                                                                                                                         ttATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTcc                                                   >  1:68380/1‑251 (MQ=255)                                                                                                                                                                                          tATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCg                                                  >  1:219828/1‑251 (MQ=255)                                                                                                                                                                                          tATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAAGCTTGTCCGCCAGGCTGGACTAAGCTATATCTCACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCTCACGCTATTcc                                                   <  2:190323/250‑1 (MQ=255)                                                                                                                                                                                           aTTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCgg                                                 >  1:375838/1‑251 (MQ=255)                                                                                                                                                                                                tCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCt                                                                    >  2:373584/1‑227 (MQ=255)                                                                                                                                                                                                tCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCt                                                                    <  1:373584/227‑1 (MQ=255)                                                                                                                                                                                                      gggTACCCCGTCATCCAGAAAGCTTCCCGGATTCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCCCTCCAGGGGAAGGCCCCTCGACCAGCACGCAGGTTGTGGTGGGCGATTCGATGGGCGAGTTGATGTTACTGGTCGGCATGTCCGATGCCGTCTTTGTTGGCGATTCTCTGGTTGAGCGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACATGCTATTCCGGGATTGATgggg                                      <  2:415969/250‑1 (MQ=255)                                                                                                                                                                                                          aCCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGc                                  >  1:277970/1‑251 (MQ=255)                                                                                                                                                                                                            cccGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCAt                                <  2:348174/251‑1 (MQ=255)                                                                                                                                                                                                             ccGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATa                               <  2:430504/251‑1 (MQ=255)                                                                                                                                                                                                                     cagAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAAc                       <  2:258740/249‑1 (MQ=255)                                                                                                                                                                                                                       gAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACtt                     >  1:131849/1‑251 (MQ=255)                                                                                                                                                                                                                              tCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGCCTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGGGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGATGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATCCTTTTAACTTTAAAGa               <  2:40555/250‑1 (MQ=255)                                                                                                                                                                                                                                         aTTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGcgcgcg   <  2:122361/251‑1 (MQ=255)                                                                                                                                                                                                                                          ttAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCgg  >  2:8376/1‑251 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                       GGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTACTGTATGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGG  >  NC_000913/3809139‑3809622

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