Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,788,688	C→T	D196D (GAC→GAT)	yibQ →	divergent polysaccharide deacetylase domain‑containing protein YibQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,788,688	0	C	T	97.7%	139.6 / ‑6.1	43	D196D (GAC→GAT)	yibQ	divergent polysaccharide deacetylase domain‑containing protein YibQ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (26/16);  minor base A (0/1);  total (26/17)
Fisher's exact test for biased strand distribution p-value = 3.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CCCTATGCCGTGGGGATCAACAACCACATGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGCACTGGCGTGAAGGTGATCAAACGGAAGGTATTCCTCGACGATTCGCAAAATGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGACCTGGCGCGTCGCAACGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTAT  >  NC_000913/3788443‑3788931                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     cccTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGgc                                                                                                                                                                                                                                                 >  2:216387/1‑251 (MQ=255)                    acaacCACAAGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGCTCTGGCGCGGCGCAACGGTTCAAcc                                                                                                                                                                                                                              >  2:132737/1‑251 (MQ=255)                              gggtagcaAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAAccattgccattg                                                                                                                                                                                                                    <  2:106921/247‑1 (MQ=255)                                  agcaAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCa                                                                                                                                                                                                                <  2:45050/251‑1 (MQ=255)                                     aaGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATcc                                                                                                                                                                                                             >  2:277768/1‑251 (MQ=255)                                                 aaCCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCTGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAAc                                                                                                                                                                                                 <  1:250386/251‑1 (MQ=255)                                                      gTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTTGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTAc                                                                                                                                                                                            <  1:366287/251‑1 (MQ=255)                                                        ttGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTAcg                                                                                                                                                                                           <  2:55652/250‑1 (MQ=255)                                                            tATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGGGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCgtgc                                                                                                                                                                                      >  2:129777/1‑250 (MQ=255)                                                               gCAGAAAGTGAGTCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCATAACCATCGGTAATACCCAGGCGAGGCGCGCCTCGCAAGGTACTGGTGTGAAGGTGATCAAACGGAAGGTGTTCCGCGACGATTCGCAGAACGAAGCGGACAACCGTGTGCAATTTAATCGCGCAATTGATCTGGCTCGACGCAACGGTTCAACCATTGCTATTGGGCATCTTCACCCTTCAACGGTACGCgtgcttc                                                                                                                                                                                   <  2:290352/251‑5 (MQ=255)                                                                   aaaGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCCATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCgtgctacaac                                                                                                                                                                                >  1:411961/1‑243 (MQ=255)                                                                             aGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTtat                                                                                                                                                                     >  1:151309/1‑251 (MQ=255)                                                                                   tGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCCACGGTACGCGTGCTCCAACAGATGGTTTATAACCTg                                                                                                                                                               >  1:309026/1‑251 (MQ=255)                                                                                         gTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCCTTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTgccgcc                                                                                                                                                          >  2:319765/1‑250 (MQ=255)                                                                                                   ttACTTCCTCGACAGCGTCACCATCGGAAATACCCAGGCGATGCGCGCCGCGCACGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCACCGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGc                                                                                                                                               >  2:51842/1‑251 (MQ=255)                                                                                                   ttACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTa                                                                                                                                                  <  1:82789/248‑1 (MQ=255)                                                                                                      cTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGCTGCGCGCCGCGCAAGGTACTGGCGTGACGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGAGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTgg                                                                                                                                            >  2:116614/1‑251 (MQ=255)                                                                                                            cGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGaa                                                                                                                                        >  2:57433/1‑249 (MQ=255)                                                                                                             gACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCACGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGCCATTACGCTGGTGaacgc                                                                                                                                     >  1:225001/1‑248 (MQ=255)                                                                                                              aCAGCGTAACCATCGGTAATACCCAAGCGCTGCGCGCCGCGCACGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGCCGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGGACGCGTGCTACAACAGATGGTTTCTAACCTGCCGCCAGAAATTACGCTGGTGAAAGcc                                                                                                                                    >  2:307353/1‑251 (MQ=255)                                                                                                                         aTCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCt                                                                                                                         >  1:148962/1‑251 (MQ=255)                                                                                                                            ggTAATACCCAGGCGATGAGCGCCGCGCAAGGCACTGGCGTGACGGTGATCAAACGGAAGGTGTTCCTCGACGATTTGCCGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCACCCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTCCAACAGATGGTTTATAACCTGCCGCCCGACATTCCGCTGGTGAAAGCCAGCAGCTTGCtga                                                                                                                       >  2:275752/1‑250 (MQ=255)                                                                                                                              tAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCtgaatg                                                                                                                    >  1:334818/1‑251 (MQ=255)                                                                                                                                                      gcAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTacac                                                                                           >  1:230097/1‑251 (MQ=255)                                                                                                                                                                ggCGTGAAGGTGAACAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCACCCTCCCCCTTCAACGGTACGCGTGCTACAACAGAGGGTTTATAACCGGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCATGCTGAAGGAACCGCAAGGTT‑ATACTTCTACACCGCCGaaaaa                                                                                 <  2:285170/251‑1 (MQ=255)                                                                                                                                                                 gCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTACACCGCCGAAAAAc                                                                                <  1:216387/251‑1 (MQ=255)                                                                                                                                                                    ggAAGGGGAACAAACGGAAGGTGTTCCCCGGCGAATCGCAGAACCAAAGGGACCTCCGTGTGGAAATTAAACGCGCAAATGAACTGGCGGGGCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTACACCGCCGAAAAACGCt                                                                             <  1:260241/250‑1 (MQ=255)                                                                                                                                                                            aTCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTg                                                                      >  1:65446/1‑250 (MQ=255)                                                                                                                                                                                  cGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCAcc                                                               <  1:57433/251‑1 (MQ=255)                                                                                                                                                                                   ggAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCCGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACcg                                                              >  1:163254/1‑251 (MQ=255)                                                                                                                                                                                             ccTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACACCAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTCCACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAAtccgt                                                    >  1:240472/1‑251 (MQ=255)                                                                                                                                                                                                 gacgaTTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCG‑GGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAAtccgttccg                                                >  1:12818/1‑251 (MQ=255)                                                                                                                                                                                                        cgcagaacgAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGt                                          <  1:198229/242‑1 (MQ=255)                                                                                                                                                                                                          cagaacgAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCtg                                                                                                                        <  2:166462/165‑1 (MQ=255)                                                                                                                                                                                                          cagaacgAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCtg                                                                                                                        >  1:166462/7‑171 (MQ=255)                                                                                                                                                                                                             aacgAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCt                                    <  1:51842/248‑1 (MQ=255)                                                                                                                                                                                                                   gCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCAGTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTGCCGTGGCGTGAAGCTGTGCaa                              <  2:230097/251‑1 (MQ=255)                                                                                                                                                                                                                         aTCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAg                                                                                                                                                                              >  1:104043/1‑102 (MQ=255)                                                                                                                                                                                                                         aTCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCga                         >  1:30684/1‑250 (MQ=255)                                                                                                                                                                                                                         aTCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTCCAACGAGTGTCATTGGGCATCCTCACCCTTCAACGGTACCGGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCTCA‑GGTTGATACTGCGACACCGCCGAATAACGCGGTGCCTTATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCgaa                        <  1:56351/251‑1 (MQ=255)                                                                                                                                                                                                                                tgCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCCGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATCCTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAGGAAACCg                 >  1:430196/1‑251 (MQ=255)                                                                                                                                                                                                                                        aaTCGCGCAATTGATCTGGCGCGACGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAACCCGATAGAAc          <  2:242862/250‑1 (MQ=255)                                                                                                                                                                                                                                               cAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTAt  >  1:157827/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     CCCTATGCCGTGGGGATCAACAACCACATGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGCACTGGCGTGAAGGTGATCAAACGGAAGGTATTCCTCGACGATTCGCAAAATGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGACCTGGCGCGTCGCAACGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTAT  >  NC_000913/3788443‑3788931

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