Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	871,513	G→A	5.1%	A183T (GCA→ACA)	gsiC →	glutathione ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	871,513	0	G	A	5.1%	254.7 / 3.5	79	A183T (GCA→ACA)	gsiC	glutathione ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (39/36);  new base A (2/2);  total (41/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.68e-02

ATGCCAACGCTGTGGCTGACCATAACCAGTATGGTCTGGGCGGTTATATTTGGTATGGCGGCGGGAATTATCGCCGCCGTCTGGCGTAACCGTTGGCCGGATCGATTGAGTATGACCATTGCGGTGTCGGGGATCTCGTTTCCGGCATTTGCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAAAGGGGTGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGGCTTTTTGCTCGGTGGTTCCATCGTTGTGGAGAAAGTTTTCAACTGGCCGGGACTTGGACGCTTACTCGTTGACTCCGTAGAAATGCGTGATTACCC  >  NC_000913/871252‑871781                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                             aTGCCAACGCTGTGGCTGACCATAACCAGTATGGTCTGGGCGGTTATGTTTGGTATGGCGGCGGGAATTATCGCCGCCGTCTGGCGTAACCGTTGGCCGGATCGATTGAGTATGACCATTGCGGTGTCGGGGATCTCGTTTCCGGCATTTGCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGccgccg                                                                                                                                                                                                                                                                  <  2:8849/274‑1 (MQ=255)       aCGCTGTGGCTGACCATAACCAGTATGGTCTGGGCGGTTATATTTGGTATGGCGGCGGGAATTATCGCCGCCGTCTGGCGTAACCGTTTGCCGGATCGATTGAGTATGACCATTGCGGTGTCGGGGATCTCGTTTCCGGCATTTGCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGTACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCTCTACATTTTACCCTCCCTGACGCCCGGCGCGGCAGTGgcccacgtcgtgg                                                                                                                                                                                                                                                            >  2:207112/1‑264 (MQ=255)                             gTATGGTCTGGGCGGTTATATTTGGTATGGCGGCGGGAATTATCGCCGCCGTCTGGCGTAACCGTTGGCCGGATCGATTGAGTATGACCATTGCGGTGTCGGGGATCTCGTTTCCGGCATTTGCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACCTTTTACCCTCCCTGACGCTCGGCGCGGCAGGGGCCGCCGAGGTGGCGCGCTTTACCCGCCcgtctgtt                                                                                                                                                                                                                                     >  2:57236/1‑271 (MQ=255)                                       ggCGGTTATATTTGGTATGGCGGCGGGAATTATCGCCGCCGTCTGGCGTAACCGTTGGCCGGATCGATTGAGTATGACCATTGCGGTGTCGGGGATCTCGTTTCCGGCATTTGCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGc                                                                                                                                                                                                                                                                             >  2:119054/1‑225 (MQ=255)                                           gTTATATTTGGTATGGCGGCGGGAATTATCGCCGCCGTCTGGCGTAACCGTTGGCCGGATCGATTGAGTATGACCATTGCGGTGTCGGGGATCTCGTTTCCGGCATTTGCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTg                                                                                                                                                                                                                                                                         <  4:119053/225‑1 (MQ=255)                                             tataTTTGGTATGGCGGCGGGAATTATCGCCGCCGTCTGGCGTAACCGTTGGCCGGATCGATTGATTATGACCATTGCGGTGTCGGGGATCTCGTTTCCGGCATTTGCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCg                                                                                                                                                                                                                     <  2:196481/275‑1 (MQ=255)                                                          ggggggggATTTAATGCGGCGGTCTGGCGGAACCGTTGGCCGGATCGATTGAGTATGACCATGGGGGTGTCGGGGATCTCGTTTCCGGCATTTGCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCGGACGCTCGGCGCGGCAGTGGCCGCCGGGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCg                                                                                                                                                                                                        <  4:206645/270‑1 (MQ=255)                                                                                  tGGCGTAACCGTTGGCCGGATCGATTGAGTATGACCATTGCGGTGTCGGGGATCTCGTTTCCGGCATTTGCTCTGGGGAGGCTTTTAATTCAGGTATTCTCCGTGGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAAAGGGGTGAg                                                                                                                                                                                <  2:237362/275‑1 (MQ=255)                                                                                             ttGGCCGGATCGATTGAGTATGACCATTGCGGTGTCGGGGATCTCGTTTCCGGCATTTGCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATgg                                                                                                                                                                                                                                                            <  2:399916/188‑1 (MQ=255)                                                                                                         aTTGAGTATGACCATTGCGGTGTCGGGGATCTCGTTTCCGGCATTTGCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAAAGGGGTGAGCGAAACCTGGGTTGTCCTCAAac                                                                                                                                                         <  2:184957/275‑1 (MQ=255)                                                                                                                       ttGCGGTGTCGGGGATCTCGTTTCCGGCATTTGCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTg                                                                                                                                                                                                                                                                         >  6:468986/1‑149 (MQ=255)                                                                                                                       ttGCGGTGTCGGGGATCTCGTTTCCGGCATTTGCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGAt                                                                                                                                                                                                                                                              >  1:131307/1‑160 (MQ=255)                                                                                                                       ttGCGGTGTCGGGGATCTCGTTTCCGGCATTTGCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCTCGTCGTTTGTCGATGTTTTAAGCGAAGATTCTATGCGTACAGCGAGGGCGAAAGGGGTGGGCGAAACCTGGGTTTTCCTCAAACACGGCCTAccgaac                                                                                                                                           >  2:135677/1‑270 (MQ=255)                                                                                                                           ggTGTCGGGGATCTCGTTTCCGGCATTTGCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGgcg                                                                                                                                                                                                                                                          <  3:131307/160‑1 (MQ=255)                                                                                                                              gtCGGGGATCTCGTTTCCGGCATTTGCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGCAGATTATATTCTTACCGCGGGGGCGAAAGGGGTGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCgataa                                                                                                                                    >  2:309101/1‑273 (MQ=255)                                                                                                                                  gggATCTCGTTTCCGGCATTTGCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGAt                                                                                                                                                                                                                                                              >  5:163402/1‑149 (MQ=255)                                                                                                                                       ctcGTTTCCGGCATTTGCTCTGGGGATGCTTTTAATTCCGGTATTCTCCGGTGAACTGGGCTGGCTGCCTACCGTAGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTGGGCCCGGCAGTGACCGCAGTGTTGGCGCGCTTTACGCGCGGTACGTTTGGCGATGTTTTAAGCTACGATTCTATAGCTACCGGGCGGGAGAAACGGGTGGTGGAAAGCATGGTTGTCCACAAACCAGGGCACCGGAACGCGTTGAAACCGGAAg                                                                                                                           >  4:162206/1‑275 (MQ=255)                                                                                                                                          gTTTCCGGCATTTGCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCt                                                                                                                                                                                                                                                      >  7:417095/1‑149 (MQ=255)                                                                                                                                          gTTTCCGGCATTTGCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCt                                                                                                                                                                                                                                                      >  5:438729/1‑149 (MQ=255)                                                                                                                                           tttCCGGCATTTGCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCtt                                                                                                                                                                                                                                                     >  7:220388/1‑149 (MQ=255)                                                                                                                                            ttCCGGCATTTGCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCttt                                                                                                                                                                                                                                                    <  7:83490/149‑1 (MQ=255)                                                                                                                                             tCCGGCATTTGCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAAAGGGGTGCGCGAAACCTGGTTTGTCCTCAAACACGGGCTACGTAACGCGATGGTCCCGGTAATGACCa                                                                                                                     >  2:133921/1‑275 (MQ=255)                                                                                                                                              ccGGCATTTGCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAAGGGGGTGAGCGACACCTGGGTTGTCCTCAAACACGGGCTACGTACCGCGATGATCCCGGTAGTGAccct                                                                                                                    >  1:60649/1‑273 (MQ=255)                                                                                                                                                 gCATTTGCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCcg                                                                                                                                                                                                                                               <  8:220388/149‑1 (MQ=255)                                                                                                                                                 gCATTTGCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGGTGTTTTAAGCGAAGATTATATGCATACCGCGAGGGCAAAAGGGGTGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCatga                                                                                                                  <  2:484569/274‑1 (MQ=255)                                                                                                                                                    tttGCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCcgcgc                                                                                                                                                                                                                                            >  7:195414/1‑149 (MQ=255)                                                                                                                                                       gCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCttt                                                                                                                                                                                                                                                    <  8:148184/138‑1 (MQ=255)                                                                                                                                                       gCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCttt                                                                                                                                                                                                                                                    >  7:148184/1‑138 (MQ=255)                                                                                                                                                                    ttttAATTCAGTGATTCTCCGTTGAACTGGGCTAGCGGCCTACCGTGGGAACAGACAGTTGGCAGCACTACCTTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAAAGGGGTGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGGCttt                                                                                              <  2:383126/275‑1 (MQ=255)                                                                                                                                                                      ttAATTCAGGTATTCTCCGTTGAGCTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAAAGGGGTGAGCGAAACCTggg                                                                                                                                                                     <  1:62640/202‑1 (MQ=255)                                                                                                                                                                      ttAATGCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAAAGGGGTGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCgat                                                                                                                                      <  2:488179/233‑1 (MQ=255)                                                                                                                                                                        aaTTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTaa                                                                                                                                                                                                                        <  6:390824/149‑1 (MQ=255)                                                                                                                                                                         aTTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAg                                                                                                                                                                                                                       >  7:244727/1‑149 (MQ=255)                                                                                                                                                                          ttCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGc                                                                                                                                                                                                                      >  8:262871/1‑149 (MQ=255)                                                                                                                                                                          ttCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGc                                                                                                                                                                                                                      >  7:372601/1‑149 (MQ=255)                                                                                                                                                                          ttCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGc                                                                                                                                                                                                                      >  6:164975/1‑149 (MQ=255)                                                                                                                                                                             aGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGaa                                                                                                                                                                                                                   >  6:284701/1‑149 (MQ=255)                                                                                                                                                                                   tctcCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATtat                                                                                                                                                                                                             <  6:438729/149‑1 (MQ=255)                                                                                                                                                                                    ttcCGTTGAACTGGGCTGGCTGCCTCCCGTGGGAGCAGACAGTTGGCAGCACTACATTATACCCTCCCTGACGCTCGGCGCGGCAGTGGCGGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAAAGGGGTGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGGCTTTTTGCTCGGTGGTTCCa                                                                              <  2:165588/274‑1 (MQ=255)                                                                                                                                                                                    ctcCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATtata                                                                                                                                                                                                            >  7:237676/1‑149 (MQ=255)                                                                                                                                                                                    ctcCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATtata                                                                                                                                                                                                            >  5:133894/1‑149 (MQ=255)                                                                                                                                                                                             aCTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGaa                                                                                                                                                                                                                   <  5:400198/133‑1 (MQ=255)                                                                                                                                                                                             aCTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGaa                                                                                                                                                                                                                   >  6:400198/1‑133 (MQ=255)                                                                                                                                                                                                  gctggctgCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAgg                                                                                                                                                                                              <  6:163402/149‑1 (MQ=255)                                                                                                                                                                                                    tggctgCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGc                                                                                                                                                                                            <  8:372601/149‑1 (MQ=255)                                                                                                                                                                                                      gctgCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAAAGGGGTGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGGCTTTTTGCTCGGTGGTTCAATCGTTGTGGAGCAAGttt                                                            >  1:136855/1‑275 (MQ=255)                                                                                                                                                                                                           cTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAAAGGGGTGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGGCTTTTTGCTCGGTGGTTCCATCGTTGTGGAGAAAGTTTTCaa                                                        <  2:449545/274‑1 (MQ=255)                                                                                                                                                                                                             aCCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAAAGGGGt                                                                                                                                                                                   >  8:343489/1‑149 (MQ=255)                                                                                                                                                                                                             aCCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAAAGGGGt                                                                                                                                                                                   <  5:164975/149‑1 (MQ=255)                                                                                                                                                                                                              ccGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGcgtcgt                                                                                                                                                                                                                                       >  6:164263/1‑96 (MQ=255)                                                                                                                                                                                                              ccGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGcgtcgt                                                                                                                                                                                                                                       <  5:164263/96‑1 (MQ=255)                                                                                                                                                                                                                 tGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAAAGGGGTGAGc                                                                                                                                                                               >  6:393201/1‑149 (MQ=255)                                                                                                                                                                                                                  ggggacAGACTGGTGGCAGACCTACAGTTTCTCCTCCCTGACGCGCGGGGCGGCAGTGGCCGCCGTGATGGCGAGCTTTGCCCGCGCGGCGTTTGTCGATGTTGTAAGCGAAGATTATATGCGTACAGCGAGGGCGAAAGGGGTGAGCGAAACCTGTGTTGTCCTCAAACACGGGCGACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGGCTTTTTGCTCGGTGGTTCCATCGTTGTGGAGAAAGTTTTCAACTGGCCgg                                                <  4:181580/270‑1 (MQ=255)                                                                                                                                                                                                                      gCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGaaa                                                                                                                                                                                        >  8:394344/1‑135 (MQ=255)                                                                                                                                                                                                                      gCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGaaa                                                                                                                                                                                        <  7:394344/135‑1 (MQ=255)                                                                                                                                                                                                                        aGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAAAGGGGTGAGCGAAACCt                                                                                                                                                                        >  6:361536/1‑149 (MQ=255)                                                                                                                                                                                                                        aGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAAAGGGGTGAGCGAAACCt                                                                                                                                                                        <  5:468986/149‑1 (MQ=255)                                                                                                                                                                                                                         gACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAAAGGGGTGAGCGAAACCTg                                                                                                                                                                       <  5:284701/149‑1 (MQ=255)                                                                                                                                                                                                                          aCAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAAAGGGGTGAGCGAAACCTgg                                                                                                                                                                      <  7:343489/149‑1 (MQ=255)                                                                                                                                                                                                                             gTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATTGCGCGCTTTACCCGCGCGTCGTTTGTGGATGCTTTAAGCGAAGATTATATGCCTACCGCGACGGCGAAAGGGGTGAGCGAAACCTGGGTTGTCCTCAAAGCCGGGCTACGTAACGCCAGGATCCCGGTAGGTACCATGATGGGCCATTCGGTTGGATTTGTGCTCTAGGGGTCGATCGTTGTTGGGCAACCTTTTAACTTGCTGGGCCTTgggctc                                     >  4:126938/1‑271 (MQ=255)                                                                                                                                                                                                                             gTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAAAGGGGTGAGCGAAACCTGGGtt                                                                                                                                                                   <  8:237676/149‑1 (MQ=255)                                                                                                                                                                                                                             gTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAAAGGGGTGAGCGAAACCTGGGtt                                                                                                                                                                   <  6:133894/149‑1 (MQ=255)                                                                                                                                                                                                                               tGGCAGCACTACATTTTACCCTCCCTGACTCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAAAGGGGTGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGGCTTTTTGCTCGGTGGTTCCATCGTTGTGGAGAAAGTTTTCAACTGGCCGGGACTTGGACGCt                                    <  2:128453/274‑1 (MQ=255)                                                                                                                                                                                                                                       cTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAAAGGGGTGAGCGAAACCTGGGTTGTCCTCAAac                                                                                                                                                         >  6:301921/1‑149 (MQ=255)                                                                                                                                                                                                                                       cTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAAAGGGGTGAGCGAAACCTGGGTTGTCCTCAAac                                                                                                                                                         >  7:110892/1‑149 (MQ=255)                                                                                                                                                                                                                                        tACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAAAGGGGTGAGCGAAACCTGGGTTGTCCTCAAaca                                                                                                                                                        <  8:195414/149‑1 (MQ=255)                                                                                                                                                                                                                                            ttttACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAAAggg                                                                                                                                                                                     <  7:455539/116‑1 (MQ=255)                                                                                                                                                                                                                                            ttttACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAAAggg                                                                                                                                                                                     >  8:455539/1‑116 (MQ=255)                                                                                                                                                                                                                                             tttACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAAAGGGGTGAGCGAAACCTGGGTTGTCCTCAAACACGGGc                                                                                                                                                   <  8:417095/149‑1 (MQ=255)                                                                                                                                                                                                                                               tACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAAAGGGGTGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGGCTTTTTGCTCGGTGGTTcaa                                                                              >  1:123091/1‑214 (MQ=255)                                                                                                                                                                                                                                                aCCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAAAGGGGTGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTAc                                                                                                                                                >  7:241524/1‑149 (MQ=255)                                                                                                                                                                                                                                                    tccctGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAAAGGGGTGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTaa                                                                                                                                            >  8:185384/1‑149 (MQ=255)                                                                                                                                                                                                                                                     ccctGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAAAGGGGTGAGCGCAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTGACAGTTTGGCTTTTTTACCGGTGGTTCCATTGTTTGGCAAAAATTTTTCCACGTGCCGGTGCTTTGAGGCTTACTCGTTGAATcgtgacaaagt             >  4:365855/1‑265 (MQ=255)                                                                                                                                                                                                                                                           cGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAAAGGGGTGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCgatg                                                                                                                                     <  7:185384/149‑1 (MQ=255)                                                                                                                                                                                                                                                              tCGGCGCGACAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAAAGGGGTGAGCGAAACCTgg                                                                                                                                                                      >  7:48603/1‑113 (MQ=255)                                                                                                                                                                                                                                                              tCGGCGCGACAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAAAGGGGTGAGCGAAACCTgg                                                                                                                                                                      <  8:179647/113‑1 (MQ=255)                                                                                                                                                                                                                                                              tCGGCGCGACAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAAAGGGGTGAGCGAAACCTgg                                                                                                                                                                      >  7:179647/1‑113 (MQ=255)                                                                                                                                                                                                                                                              tCGGCGCGACAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAAAGGGGTGAGCGAAACCTgg                                                                                                                                                                      <  8:48603/113‑1 (MQ=255)                                                                                                                                                                                                                                                                ggCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGATTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAAAGGTGTGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGGCTTTTTGCTCGGTGGTTCCATCGTTGTGGAGAAAGTTTTCAACTGGCCGGGACTTGGACGCTTACTCGTTGACTCCGTAGAAATGCGTGATTAccc  >  2:378900/1‑275 (MQ=255)                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                             ATGCCAACGCTGTGGCTGACCATAACCAGTATGGTCTGGGCGGTTATATTTGGTATGGCGGCGGGAATTATCGCCGCCGTCTGGCGTAACCGTTGGCCGGATCGATTGAGTATGACCATTGCGGTGTCGGGGATCTCGTTTCCGGCATTTGCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTTTGTCGATGTTTTAAGCGAAGATTATATGCGTACCGCGAGGGCGAAAGGGGTGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGGCTTTTTGCTCGGTGGTTCCATCGTTGTGGAGAAAGTTTTCAACTGGCCGGGACTTGGACGCTTACTCGTTGACTCCGTAGAAATGCGTGATTACCC  >  NC_000913/871252‑871781

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