Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	871,607	G→A	S214N (AGC→AAC)	gsiC →	glutathione ABC transporter membrane subunit GsiC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	871,607	0	G	A	97.8%	150.4 / NA	46	S214N (AGC→AAC)	gsiC	glutathione ABC transporter membrane subunit GsiC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (20/25);  minor base T (0/1);  total (20/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GTATGACCATTGCGGTGTCGGGGATCTCGTTTCCGGCATTTGCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAA‑AGGGGTGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGGCTTTTTGCTCGGTGGTTCCATCGTTGTGGAGAAAGTTTTCAACTGGCCGGGACTTGGACGCTTACTCGTTGACTCCGTAGAAATGCGTGATTACCCGGTGATTCAGGCGGAAATTCTGCTTTTCTCGCTGGAATTTATTCTTATCAACTTAGTGGTGGATGTG  >  NC_000913/871361‑871848                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  gtgtgcacATTGGGGGGTCGTGGGTCTCGTTGCCGGCATTTGCTCTGGGGATGATTTGAATTCCGGTATTTTCGGTTGAAAGGGGCGTGATGCCTCCCGTGGGTGCGGACCGGTGGCAGCCCTACATTTCCCCCTCCTCGACGCTCGGCGCGGCAGTGGCCTCCTTGATGGCGCGCTTAACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACTTCGAGCGCGAA‑AGGGGTGAACGaa                                                                                                                                                                                                                                               <  2:372050/244‑1 (MQ=255)              ggTGTCGGGGATCTCGGTTCCGGCATTTGCTATGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACAGCGAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTcc                                                                                                                                                                                                                                  >  2:353598/1‑251 (MQ=255)                    ggggATCTCGTTTCCGGCATTTGCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAac                                                                                                                                                                                                                            <  1:182751/251‑1 (MQ=255)                            cGTTTCCGGCATTTGCTCTGGGTATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCTAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACgg                                                                                                                                                                                                                        <  2:35779/247‑1 (MQ=255)                                                     gCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTg                                                                                                                                                                                            >  1:278862/1‑250 (MQ=255)                                                     gCTTTTAATTCAGGTATTCTCCGTTGAACGGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTCCCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGa                                                                                                                                                                                           <  2:391570/251‑1 (MQ=255)                                                       ttttAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGACCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAA‑AGGGGTTAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGGAACGCGATGATCCCGGTAGTGAcc                                                                                                                                                                                         >  2:226039/1‑251 (MQ=255)                                                               cAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATgg                                                                                                                                                                                  <  1:116277/250‑1 (MQ=255)                                                                aGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGAGGGGc                                                                                                                                                                                >  1:325164/1‑251 (MQ=255)                                                                aGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGAGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGGTTTAAGCGGAGGTTATATGCGTACCGCGAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAGCGAGATGATCCCGGTAGTGACCATGATGGGc                                                                                                                                                                                >  2:316203/1‑251 (MQ=255)                                                                     ttCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTAc                                                                                                                                                                            >  2:87450/1‑250 (MQ=255)                                                                        tctgTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGtt                                                                                                                                                                        >  2:375915/4‑251 (MQ=255)                                                                           gTTGAAATGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg                                                                                                                                                                     <  1:226039/251‑1 (MQ=255)                                                                                aCTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTCCCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGGCtttt                                                                                                                                                                <  1:201369/251‑1 (MQ=255)                                                                                        ggctgCTTACCGGGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTGCCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGGCTTTTTGCTCg                                                                                                                                                          <  2:367060/249‑1 (MQ=255)                                                                                            gCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTGTGCCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGGCTTTTTGCTCggtgg                                                                                                                                                      <  2:83394/249‑1 (MQ=255)                                                                                                       gAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGGCTTTTTGCTCGGTGGTTCCATCGTtgtg                                                                                                                                         >  1:155094/1‑251 (MQ=255)                                                                                                           aGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACACGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGGCTTTTTGCTCGGTGGTTCCATCGTTGTgg                                                                                                                                        >  1:208565/1‑248 (MQ=255)                                                                                                             gcAGTTGGCAGCCATACATTTTACCCTCCCTTACGCGCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTGGTCGAGGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGGCTTTTTGCTCGGTGGTTCCATCGTTGTGGAGaaa                                                                                                                                   <  2:4955/250‑1 (MQ=255)                                                                                                              cAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGCAGCTTATATGCGTACCGCGAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGGCTTTTTGCTCGGTGGTTCCCTCGTTGTGGAGaaa                                                                                                                                   >  1:4955/1‑250 (MQ=255)                                                                                                                       gcaCTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGGCTTTTTGCTCGGTGGTTCCATCGTTGTGGAGAAAGTTTTCAACt                                                                                                                         <  1:7725/251‑1 (MQ=255)                                                                                                                                tttACCCGCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTCATGGCGCCCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGGGACCATGATGGGCTTACAGTTTGGCTTTTTGCTCGGTGGTTCCATCGTTGTGGAGAAAGTTTTCATCTGGCCgggcc                                                                                                                >  2:52488/1‑249 (MQ=255)                                                                                                                                           tGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATggg                                                                                                                                                                                 >  1:71381/1‑175 (MQ=255)                                                                                                                                           tGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATggg                                                                                                                                                                                 <  2:71381/175‑1 (MQ=255)                                                                                                                                                                    cgTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGGCTTTTTGCTCGGTGGTTCCATCGTTGTGGAGAAAGTTTTCAACTGGCCGGGACTTGGACGCTTACTCGTTGACTCCGTAGAAATGCGTg                                                                            <  1:87450/251‑1 (MQ=255)                                                                                                                                                                     gTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGGCTTTTTGCTCGGTGGTTCCATCGTTGTGGAGAAAGTTTTCAACTGGCCGGGACTTGGACGCTTACTCGTTGACTCCGTAGAAATGCGTg                                                                            >  1:260144/1‑250 (MQ=255)                                                                                                                                                                       gATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGGCTTTTTGCTCGGTGGTTCCATCGTTGTGGAGAAAGTTTTCAACTGGCCGGGACTTGGACGCTTACTCGTTGACTCCGTAGAAATGCGTGAt                                                                          <  2:247219/250‑1 (MQ=255)                                                                                                                                                                             ggncTTTACCCGCGCGTCGTTTGGCGCTGTTTTAAGCGAAGATGATATTCGTACCGCGATAGCTAACACGGG‑GATCGAATCCTGGGTTGTCCTCAACCACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACTGTTTGGCTGTTTGCTCGGTGTTTCAATCGTTGTGGAGAAAGTTTCCACCGTGCCGGGTCTTGGCCGCTTACTCGTTGTCTCCGTAGAAATGCGTGATTatcct                                                                    <  2:421925/247‑5 (MQ=255)                                                                                                                                                                             gcgcTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGGCTTTTTGCTCGGTGGTTCCATCGTTGTGGAGAAAGTTTTCAACTGGCCGGGACTTGGACTCTTACTCGTTGACTCCGTAGAAATGCGTGATTACCCgg                                                                   <  1:77267/251‑1 (MQ=255)                                                                                                                                                                                         cgcgTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGt                                                                                                                                                                         >  1:196960/1‑137 (MQ=255)                                                                                                                                                                                         cgcgTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGt                                                                                                                                                                         <  2:196960/137‑1 (MQ=255)                                                                                                                                                                                              cgtTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGGCTTTTTGCTCGGTGGTTCCATCGTTGTGGAGAAAGTTTTCAACTGGCCGGGACTTGGACGCTTACTCGTTGACTCCGTAGAAATGCGTGATTACCCGGTGATTCAGGCGGAAAtt                                                  >  1:36722/1‑251 (MQ=255)                                                                                                                                                                                                tttGTCGTTGTTTTAAGCGACGATTATATGCGTACCGCGAGGGCGTA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTAAAGTTTGGCTTTTTGATCGGTGGTTCCATCGTTGTGGAGAAAGTTTTAAACTGGCCGGGACTTGGACGCTTACTCGTTGACTCCGTAGAAATGCGTGATTACCCGGTGATTCAGGCGGACATTCt                                                >  2:5854/1‑251 (MQ=255)                                                                                                                                                                                                   gTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGGCTTTTTGCTCGGTGGTTCCATCGTTGTGGAGAAAGTTTTCAACTGGCCGGGACTTGGACGCTTACTCGTTGACTCCGTAGAAATGCGTGATTACCCGGTGATTCAGGCGGAAATTCTGCt                                             >  1:365726/1‑251 (MQ=255)                                                                                                                                                                                                      gATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGGCTTTTTGCTCGGTGGTTCCATCGTTGTGGAGAAAGTTTTCAACTGGCCGGGACTTGGACGCTTACTCGTTGACTCCGTAGAAATGCGTGATTACCCGGTGATTCAGGCGGAAATTCTGCtttt                                          <  1:330263/251‑1 (MQ=255)                                                                                                                                                                                                             tAAGCGAAGATTATATGCGTACCGCGAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGGCTTTTTGCTCGGTGGTTCCATCGTTGTGGAGAAAGTTTTCAACTGGCCGGGACTTGGACGCTTACTCGTTGACTCCGTAGAAATGCGTGATTACCCGGTGATTCAGGCGGAAATTCTGCTTTTCTCGCt                                    <  1:64074/250‑1 (MQ=255)                                                                                                                                                                                                              aaGCGAAGATTATATGCGTACCGCGAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGGCTTTTTGCTCGGTGGTTCCATCGTTGTGGAGAAAGTTTTCAACTGGCCGGGACTTGGACGCTTACTCGTTGACTCCGTAGAAATGCGTGATTACCCGGTGATTCAGGCGGAAATTCTGCTTTTCTCGCTg                                   <  1:52488/250‑1 (MQ=255)                                                                                                                                                                                                               aGCGAAGATTATATGCGTACCGCGAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGGCTTTTTGCTCGGTGGTTCCATCGTTGTGGAGAAAGTTTTCAACTGGCCGGGACTTGGACGCTTACTCGTTGACTCCGTAGAAATGCGTGATTACCCGGTGATTCAGGCGGAAATTCTGCTTTTCTCGCTGGa                                 <  1:336536/251‑1 (MQ=255)                                                                                                                                                                                                               aGCGAAGATTATATGCGTACCGCGAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGGCTTTTTGCTCGGTGGTTCCATCGTTGTGGAGAAAGTTTTCAACTGGCCGGGACTTGGACGCTTACTCGTTGACTCCGTAGAAATGCGTGATTACCCGGTGATTCAGGCGGAAATTCTGCTTTTCTCGCTGGa                                 <  2:123968/251‑1 (MQ=255)                                                                                                                                                                                                                gCGAAGATTATATGCGTACCGCGAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGGCTTTTTGCTCGGTGGTTCCATCGTTGTGGAGAAAGTTTTCAACTGGCCGGGACTTGGACGCTTACTCGTTGACTCCGTAGAAATGCGTGATTACCCGGTGATTCAGGCGGAAATTCTGCTTTTCTCGCTGGa                                 >  1:386720/1‑250 (MQ=255)                                                                                                                                                                                                                   aaGATTATATGCGTACCGCGAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGGCTTTTTGCTCGGTGGTTCCATCGTTGTGGAGAAAGTTTTCAACTGGCCGGGACTTGGACGCTTACTCGTTGACTCCGTAGAAATGCGTGATTACCCGGTGATTCAGGCGGAAATTCTGCTTTTCTCGCTTGAAttt                             <  1:99401/251‑1 (MQ=255)                                                                                                                                                                                                                   aaGATTATATGCGTACCGCGAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGGCTTTTTGCTCGGTGGTTCCATCGTTGTGGAGAAAGTTTTCAACTGGCCGGGACTTGGACGCTTACTCGTTGACTCCGTAGAAATGCGTGATTACCCGGTGATTCAGGCGGAAATTCTGCTTTTCTCGCTGGAAttt                             <  2:155094/251‑1 (MQ=255)                                                                                                                                                                                                                           atGCGTACCGCGAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGGCTTTTTGCTCGGTGGTTCCATCGTTGTGGAGAAAGTTTTCAACTGGCCGGGACTTGGACGCTTACTCGTTGACTCCGTAGAAATGCGTGATTACCCGGTGATTCAGGCGGAAATTCTGCTTTTCTCGCTGGAATTTATTCTTAt                     >  1:88391/1‑251 (MQ=255)                                                                                                                                                                                                                             gCGTACCGCGAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGGCTTTTTGCTCGGTGGTTCCATCGTTGTGGAGAAAGTTTTCAACTGGCCGGGACTTGGACGCTTACTCGTTGACTCCGTAGAAATGCGTGATTACCCGGTGATTCAGGCGGAAATTCTGCTTTTCTCGCTGGAATTTATTCTTATc                    <  1:316203/250‑1 (MQ=255)                                                                                                                                                                                                                                      gAGGGCGAA‑AGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGGCTTTTTGCTCGGTGGTTCCATCGTTGTGGAGAAAGTTTTCAACTGGCCGGGACTTGGACGCTTACTCGTTGACTCCGTAGAAATGCGTGATTACCCGGTGATTCAGGCGGAAATTCTGCTTTTCTCGCTGGAATTTATTCTTATCAACTTAgtgg          >  2:83466/1‑251 (MQ=255)                                                                                                                                                                                                                                              a‑aGGGGTGAACGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGGCTTGTTGCTCGGTGGTTCCATCGTTGTGGAGAAAGTTTTCAACTGGCCGGGACTTGGACGCTTACTCGTTGACTCCGTAGAAATGCGTGATTACCCGGTGATTCAGGCGGAAATTCTGCTTTTCTCGCTGGAATTTATTCTTATCAACTTAGTGGTGGAtgtg  <  2:150321/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  GTATGACCATTGCGGTGTCGGGGATCTCGTTTCCGGCATTTGCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAA‑AGGGGTGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGGCTTTTTGCTCGGTGGTTCCATCGTTGTGGAGAAAGTTTTCAACTGGCCGGGACTTGGACGCTTACTCGTTGACTCCGTAGAAATGCGTGATTACCCGGTGATTCAGGCGGAAATTCTGCTTTTCTCGCTGGAATTTATTCTTATCAACTTAGTGGTGGATGTG  >  NC_000913/871361‑871848

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