Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,937	A→G	P68P (CCA→CCG)	thrC →	threonine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,937	0	A	G	82.9%	64.2 / 3.8	36	P68P (CCA→CCG)	thrC	threonine synthase
Reads supporting (aligned to +/- strand):  ref base A (3/3);  new base G (15/14);  total (18/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-01

CGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATT‑CAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCA‑CAGGAAATCCTGGAAGAGCGCGTGCGCGC‑GGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGTGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTA  >  NC_000913/3695‑4177                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   cGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATT‑CAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCG‑CAGGAAAt                                                                                                                                                                                                                                           <  1:194377/251‑1 (MQ=255)                   caaccgtgaatGGAAATCAAAAAGAAACTCTAACTTCTGAACGAACACAACGATAAGGTCATCTTTGCGCTAGCCGTAACCCAGTGGTTGGGTAAAACTCAGGGTCTGTTGTTTCCTCACGTTCTGCCGTAATT‑CAGCTTGAATGAAATTGATGATAGGCTGAAGCTGGATTTTGTCACCCGCAGTGCTAAGATCCGCTCGGCGTGTATTGGTTATGAAATCTCT‑CAGTAATTCCTGAAATAGCGCGTgcg                                                                                                                                                                                                                         <  2:253234/241‑1 (MQ=255)                              ggAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATT‑CAGCCTGACCGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCG‑CAGGAAATCCTGGAAGAGCGCGTGCGCGC‑GGCGTTTg                                                                                                                                                                                                             >  1:368723/1‑251 (MQ=255)                                           aaaCTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATT‑CAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCA‑CAGGAAATCCTGGAAGAGCGCGTGCGCGC‑GGCGTTTGCCTTCCCGGCTcc                                                                                                                                                                                                >  1:209003/1‑251 (MQ=255)                                              ctctACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGGATGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGCATT‑CAGCCTGACTGAAATTGATGAGATGCTGCAGCTGGATTTTGTCACCCGCAGCGCGAAGATACTCTCGGCGTTTATTGGTGATGCAATCCCG‑CAGGAAATCATGGCAGAGCGCGTGCGCGC‑GGCATTTGCCTTCACGGCACCGGt                                                                                                                                                                                             >  2:337411/1‑251 (MQ=255)                                                ccaCAATCTGAAAGATCACAATGAGCAGGGCAGCGTTGGGCAAGCCGTAGCCCAGGGGTTGGGCAAAAAGCAGGGGCGGTTTTTTCCGCACGACCTGCCGGAATT‑CAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCG‑CAGGAAATCCTGGAAGAGCGCGTGCGCGC‑GGCGTTTGCATTCCCGGCTCCGGt                                                                                                                                                                                             <  2:144861/247‑1 (MQ=255)                                                             gATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATT‑CAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCG‑CAGGAAATCCTGGAAGAGCGCGTGCGCGC‑GGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGaa                                                                                                                                                                                >  1:149506/1‑249 (MQ=255)                                                              aTCACAATGTGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGTGGATGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGACGGAATTACA‑CCTGACTGAAATTGCAGAAATGCTGAAGCTGGATTTTGTCAACCGCCGTGCGCAGATCCTCTCGGCGTTTATTAGTGATGAAATCCCG‑TAGGAAATCCTGCATGAGCCCCTGAGCGA‑GGCGTTTGCCTTGCGGGCTCCGTTCGCCaaggtggcaggc                                                                                                                                                                             >  2:287415/1‑241 (MQ=255)                                                              aTCACAATGAGAAGGTCCGATTTGCGCAAGCCGTACCACAGGGGTTAGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTACCGGAATT‑CAGACTGAATGAAATGAATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGGGCGAAGATCCTCTCTGCGTTTATTGGTGGTGAAATCACG‑CGGGAAATCCTGGAAAAGCGCGTGC‑CGCTGGCGTTTCCATTCCCGACTCCGGTCGCCAATggtgcaagg                                                                                                                                                                             >  2:287844/1‑243 (MQ=255)                                                                      gAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGAGGTTGGGCAAGAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATT‑CAGCCTGACTGAACTTGCTGAGATGCTGAAGCTGGTTTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGATATGAAATCCCG‑CAGGAAATACTGGAAGAGCGCGTGCGCGC‑GGCGTTTGACTTCACTGCTCCGGTCGCCAATGTTGGAAGCGaggaagg                                                                                                                                                                     >  2:325085/1‑245 (MQ=255)                                                                                                    cAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATT‑CAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCG‑CAGGAAATCCTGGAAGAGCGCGTGCGCGC‑GGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCt                                                                                                                                       >  1:304745/1‑251 (MQ=255)                                                                                                     aGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATT‑CAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCG‑CAGGAAATCCTGGAAGAGCGCGTGCGCGC‑GGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTg                                                                                                                                      >  1:147381/1‑251 (MQ=255)                                                                                                           tGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGCCCTGCCGGAATT‑CAGCCTGTCTGAAATTGAAGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCG‑CAGGAAATCCTGGAAGAGCGCGTGCGCgc‑gg                                                                                                                                                                                                                   <  2:281213/168‑1 (MQ=255)                                                                                                           tGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATT‑CAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCG‑CAGGAAATCCTGGAAGAGCGCGTGCGCgc‑gg                                                                                                                                                                                                                   >  1:281213/1‑168 (MQ=255)                                                                                                            gggCAAAAATCAGGGGCTGGTTTTTCCGCACGACCTGCCGGAATT‑CAGCCTGACTGAAATTGATGAGATGCGGAAGCTGGATTTTGTCCCCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCG‑CAGGAAATCCTGGAAGAGCGCGTGCGCGC‑GGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTggtatgt                                                                                                                                <  2:207984/250‑6 (MQ=255)                                                                                                                     tCAGGAGCTGTTTTTTCCGAAAGCCCTGCCGGAATT‑CAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGCTTTTGTCACCCGCAGTGAGAAGCTCCTCTAGGCGTTTCTTGGAGTTGAAATCCCG‑CAGGAAATCCTGGAATAGCGCGTGCGCAC‑GGCGTTTGCCTTCCAGGATCCGGTCGACAATGTTGAACGCGATGTCGGTTGTCTGGACTTGTTCCACGTGCAAACGCTGGCATTTACAGAttt                                                                                                                        >  2:167430/1‑249 (MQ=255)                                                                                                                         gggCTGTTTTTTCCGCACGACCTGCCGGAATT‑CAGCCTGACTGCAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGGTGAAATCCCA‑CAGGAAATCCTGCACGAGCGCGTGCGCGC‑GGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAACGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAACGATTTcgtctg                                                                                                                  >  1:22408/1‑247 (MQ=255)                                                                                                                              gTTTTTTCCGCACGACCTGCCGGAATT‑CAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCTG‑CAGGAAATCCTGGAAGAGCGCGTGCGCGC‑GGCGTTTGTCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCg                                                                                                               <  1:106323/249‑1 (MQ=255)                                                                                                                                tttttCCGCACGACCTGCCGGAAGT‑CAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCTCCCGCAGGGCGAAGCTCCTCTCGGCGTTTATTGGTGATGAAATCCCG‑CAGGAAATCCTGGAAGAGCGCGTGCGCGC‑GGCGTTTGCCTTCCCGGCGCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCtt                                                                                                            <  2:52712/250‑1 (MQ=255)                                                                                                                                   ttCCGCACGACCTGCCGGAATT‑CAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCG‑CAGGAAATCCTGGAAGAGCGCGTGCGCGC‑GGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTAt                                                                                                         <  1:223812/250‑1 (MQ=255)                                                                                                                                     ccGCACGACCTGCCGGAATT‑CAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCG‑CAGGAAATCCTGGAAGAGCGCGTGCGCGC‑GGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGt                                                                                                                                                                        <  1:382294/185‑1 (MQ=255)                                                                                                                                     ccGCACGACCTGCCGGAATT‑CAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCG‑CAGGAAATCCTGGAAGAGCGCGTGCGCGC‑GGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGt                                                                                                                                                                        >  2:382294/1‑185 (MQ=255)                                                                                                                                             ccTGCCGGAATT‑CAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCG‑CAGGAAATCCTGGAAGAGCGCGTGCGCGC‑GGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCTCTTTCTGGCTCAAATGc                                                                                              >  2:121456/1‑251 (MQ=255)                                                                                                                                                              attactgaAATTGATGGGATGCTGAAGATGGATTTTCTCACCCGCCGTGCGAAGCTCCTCTCGTCGTTTATTGGTG‑CGGAATCCCAGTATGAAATCCTGGAAGAGCGCGTGCGCGG‑GGCGTTGGCTTCCCCGGCTCCGGTCGCCAATGTTGAATGCGATGTCGGTTGTCTGGAATTTTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGTCCCATATTGCggg                                                                              <  2:258904/248‑1 (MQ=255)                                                                                                                                                                                     gAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCG‑CAGGAAATCCTGGAAGAGCGCGTGCGCGC‑GGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGACCAAAAGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTgg                                                       >  1:236159/1‑250 (MQ=255)                                                                                                                                                                                                         gctgTGCGAATATCCTCTCGGCGTGTATTGTTTATGAAACCCCG‑CTGGAAATCCTGGAATAGCGCGTGCGCTC‑GGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTTTTCCACGTGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCATTGCCCATTCTGACCGCGTCCTCCGGTGATAcc                                   <  2:141496/248‑1 (MQ=255)                                                                                                                                                                                                          cagtgtgaagattcgtgcggggTTTATTGGGGATGATACCCCG‑GAGGAATGCCTGGGAGAGCGCGTGCGCGC‑TGCGTTTGCTTGTCCGGCACCGGTGGCCAATGGTGAAAGCGATGTCGGTTGTCTGTAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCGCAAATGCTGACCCATATGGCGTGCGATAATCCATTGACCAGTCGGACCTGGACCTCCTGTTATACCg                                  <  2:245948/230‑1 (MQ=255)                                                                                                                                                                                                              gCGAAGACCCTCTCTGCGTTGATGGGTGATGAAAGCCCG‑CAGGAAATCCTGGAAGAGCTCGGGCGCGC‑GGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTGGTCTGGAATGGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCTACCTCCGGTGATACCGGAg                               <  2:33102/250‑1 (MQ=255)                                                                                                                                                                                                               cGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCG‑CAGGAAATCCTGGAAGAGCGCGTGCGCGC‑GGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCg                             <  1:325085/251‑1 (MQ=255)                                                                                                                                                                                                                    aTCCTCGCGGCGTTTAGTGGTGATGAAATCCCG‑CAGGAAATCCTGGAAGAGCGCGTGCGCGC‑GGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCAGATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCGCCGGTGATACCGGAGCGGCAGt                        <  2:149506/251‑1 (MQ=255)                                                                                                                                                                                                                           cggggTTTTTTGGTGACGAAAGCCCG‑CAGGAACTTCTGGAAGAGCTCGTGCGCGC‑GGCGTTTGCTGTCCCGGCTCCGGGCGCCAATGTTGAAAGCTATGTCGGTTGTCTGGAATTGTTCCTCGGGCCAACGCTGGCATTTAAAGATTTCTGCGGGCTCTTTATGGCACAAATGCTGACCCATGTTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATCCCGGAGCGGCAGTGGCTCAt                 <  2:137976/247‑1 (MQ=255)                                                                                                                                                                                                                           cGGCGTTTATTGGTTATGAAATCCCA‑CAGGAAATCCTGGAAGAGCGCGGGCTCGC‑GGCGTTTGCCTTCCCGGCTTCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGTGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCAt                 <  2:209003/251‑1 (MQ=255)                                                                                                                                                                                                                               gTTTATTGGTGATGAAATCCCG‑CAGGAAATCCTGGAAGAGCGCGTGCGCGC‑GGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAACTGCTGACCCATATTGCGGGCGATAAGCCATTGACCATTCTTAACGCGACCTCCGGTGATACCGGAGCGGCAGTTGCTCATGCtt             >  1:275149/1‑251 (MQ=255)                                                                                                                                                                                                                                 ttATTGGTGATGAAATCCCA‑CAGGAAATCCTGGAAGAGCGCGTGCGCGC‑GGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCCTATTGCGGGTGATCAGCCCGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCCCATGCTTTc           >  1:132830/1‑251 (MQ=255)                                                                                                                                                                                                                                    ttGGTGATGAAATCCCG‑CAGGAAATCCTGGAAGAGCGCGTGCGCGC‑GGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGGGGCTCATGCTTTCTAc        >  1:308764/1‑251 (MQ=255)                                                                                                                                                                                                                                          atgaAATCCCA‑CAGGAAATCCTGGAAGAGCGCGGGCGCGC‑GGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGTGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTa  <  2:132830/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   CGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATT‑CAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCA‑CAGGAAATCCTGGAAGAGCGCGTGCGCGC‑GGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGTGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTA  >  NC_000913/3695‑4177

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