Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	260,466	C→T	5.5%	H27Y (CAT→TAT)	proB →	gamma‑glutamate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	260,466	0	C	T	5.5%	236.9 / 2.9	72	H27Y (CAT→TAT)	proB	gamma‑glutamate kinase
Reads supporting (aligned to +/- strand):  ref base C (36/32);  new base T (2/2);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

CAACCGACGACAGTCCTGCTAAAACGTTCGTTTGATATCATTTTTCCTAAAATTGAATGGCAGAGAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCT  >  NC_000913/260319‑260611                                                                                                                                                    |                                                                                                                                                   cAACCGACGACAGTCCTGCTAAAACGTTCGTTTGATATCATTTTTCCTAAAATTGAATGGCAGAGAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCa                                                                                                                                                   >  1:53627/1‑149 (MQ=255)  aaCCGACGACAGTCCTGCTAAAACGTTCGTTTGATATCATTTTTCCTAAAATTGAATGGCAGAGAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCat                                                                                                                                                  >  3:168726/1‑149 (MQ=255)  aaCCGACGACAGTCCTGCTAAAACGTTCGTTTGATATCATTTTTCCTAAAATTGAATGGCAGAGAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCat                                                                                                                                                  >  4:207858/1‑149 (MQ=255)  aaCCGACGACAGTCCTGCTAAAACGTTCGTTTGATATCATTTTTCCTAAAATTGAATGGCAGAGAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCat                                                                                                                                                  >  1:352644/1‑149 (MQ=255)           cAGTCCTGCTAAAACGTTCGTTTGATATCATTTTTCCTAAAATTGAATGGCAGAGAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGaa                                                                                                                                         <  1:235414/149‑1 (MQ=255)            aGTCCTGCTAAAACGTTCGTTTGATATCATTTTTCCTAAAATTGAATGGCAGAGAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAAc                                                                                                                                        >  4:267957/1‑149 (MQ=255)              tCCTGCTAAAACGTTCGTTTGATATCATTTTTCCTAAAATTGAATGGCAGAGAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGTTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACtt                                                                                                                                      <  2:69355/149‑1 (MQ=255)              tCCTGCTAAAACGTTCGTTTGATATCATTTTTCCTAAAATTGAATGGCAGAGAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACtt                                                                                                                                      >  1:351710/1‑149 (MQ=255)               ccTGCTAAAACGTTCGTTTGATATCATTTTTCCTAAAATTGAATGGCAGAGAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGccc                                                                                                                                                    >  2:304410/1‑134 (MQ=255)               ccTGCTAAAACGTTCGTTTGATATCATTTTTCCTAAAATTGAATGGCAGAGAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGccc                                                                                                                                                    <  1:304410/134‑1 (MQ=255)               ccTGCTAAAACGTTCGTTTGATATCATTTTTCCTAAAATTGAATGGCAGAGAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTg                                                                                                                                     <  2:439163/149‑1 (MQ=255)                 tGCTAAAACGTTCGTTTGATATCATTTTTCCTAAAATTGAATGGCAGAGAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGtt                                                                                                                                   >  4:464971/1‑149 (MQ=255)                  gCTAAAACGTTCGTTTGATATCATTTTTCCTAAAATTGAATGGCAGAGAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTc                                                                                                                                  <  3:462381/149‑1 (MQ=255)                    tAAAACGTTCGTTTGATATCATTTTTCCTAAAATTGAATGGCAGAGAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGc                                                                                                                                >  1:417857/1‑149 (MQ=255)                    tAAAACGTTCGTTTGATATCATTTTTCCTAAAATTGAATGACAGAGAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGc                                                                                                                                >  2:34985/1‑149 (MQ=255)                               tttGATATCATTTTTCCTAAAATTGAATGGCAGAGAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCa                                                                                                                     >  1:87653/1‑149 (MQ=255)                               tttGATATCATTTTTCCTAAAATTGAATGGCAGAGAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCa                                                                                                                     >  3:48625/1‑149 (MQ=255)                                  gATATCATTTTTCCTAAAATTGAATGGCAGAGAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGtt                                                                                                                  >  3:306153/1‑149 (MQ=255)                                    tatCATTTTTCCTAAAATTGAATGGCAGAGAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTAc                                                                                                                >  1:74786/1‑149 (MQ=255)                                     atCATTTTTCCTAAAATTGAATGGCAGAGAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACa                                                                                                               >  3:297946/1‑149 (MQ=255)                                      tCATTTTTCCTAAAATTGAATGGCAGAGAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACAt                                                                                                              <  4:32070/149‑1 (MQ=255)                                      tCATTTTTCCTAAAATTGAATGGCAGAGAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACAt                                                                                                              >  4:453041/1‑149 (MQ=255)                                         tttttCCTAAAATTGAATGGCAGAGAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATgcc                                                                                                           >  4:219236/1‑149 (MQ=255)                                                tAAAATTGAATGGCAGAGAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGccc                                                                                                                                                    <  2:108439/101‑1 (MQ=255)                                                tAAAATTGAATGGCAGAGAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGccc                                                                                                                                                    >  1:108439/1‑101 (MQ=255)                                                tAAAATTGAATGGCAGAGAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATTGTTGAACGTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGc                                                                                                    >  2:41910/1‑149 (MQ=255)                                                tAAAATTGAATGGCAGAGAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATgcccggcggc                                                                                                    >  1:104138/1‑142 (MQ=255)                                                  aaaTTGAATGGCAGAGAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCAt                                                                                                  <  2:53627/149‑1 (MQ=255)                                                   aaTTGAATGGCAGAGAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATc                                                                                                 >  3:34036/1‑149 (MQ=255)                                                       gAATGGCAGAGAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGAt                                                                                             <  3:13664/149‑1 (MQ=255)                                                        aaTGGCAGAGAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGAtt                                                                                            <  2:357309/149‑1 (MQ=255)                                                           ggCAGAGAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGtt                                                                                         <  1:83849/149‑1 (MQ=255)                                                           ggCAGAGAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGtt                                                                                         <  3:333173/149‑1 (MQ=255)                                                           ggCAGAGAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGtt                                                                                         >  4:377474/1‑149 (MQ=255)                                                           ggCAGAGAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGtt                                                                                         <  4:530122/149‑1 (MQ=255)                                                              agcgaATCATGAGTGACAGCCAGACGCTGGTGCTAAATCTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGCGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTAtt                                                                                      <  1:256112/146‑1 (MQ=255)                                                              agagAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTAtt                                                                                      <  2:79455/149‑1 (MQ=255)                                                              agagAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTAtt                                                                                      >  2:423759/1‑149 (MQ=255)                                                              agagAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTAtt                                                                                      >  3:230853/1‑149 (MQ=255)                                                               gagaATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCTATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATtg                                                                                     >  2:394291/1‑149 (MQ=255)                                                               gagaATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCTATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATtg                                                                                     >  2:158298/1‑149 (MQ=255)                                                               gagaATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGCTTGTTATtg                                                                                     >  1:191515/1‑149 (MQ=255)                                                                agaATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATtgt                                                                                    <  4:48625/149‑1 (MQ=255)                                                                agaATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATtgt                                                                                    <  4:74011/149‑1 (MQ=255)                                                                agaATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATtgt                                                                                    <  3:453041/149‑1 (MQ=255)                                                                agaATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATtgt                                                                                    <  1:105274/149‑1 (MQ=255)                                                                          gTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGaa                                                                                                                                         >  1:369795/1‑86 (MQ=255)                                                                          gTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGaa                                                                                                                                         <  2:369795/86‑1 (MQ=255)                                                                          gTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGgc                                                                          <  2:74786/149‑1 (MQ=255)                                                                          gTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGgc                                                                          >  2:218879/1‑149 (MQ=255)                                                                                   aGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATcgcc                                                                >  3:301822/1‑149 (MQ=255)                                                                                          ggtggtAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCCGGCGGCGATcgcgcccgac                                                          <  2:353705/149‑8 (MQ=255)                                                                                              gtAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGa                                                     <  2:417857/149‑1 (MQ=255)                                                                                              gtAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGa                                                     >  3:570097/1‑149 (MQ=255)                                                                                                 aaaCTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATc                                                                                                 >  2:39383/1‑103 (MQ=255)                                                                                                 aaaCTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATc                                                                                                 <  1:39383/103‑1 (MQ=255)                                                                                                  aaCTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCAc                                                 >  3:39261/1‑149 (MQ=255)                                                                                                     tCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTg                                              >  4:189128/1‑149 (MQ=255)                                                                                                     tCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTg                                              >  2:199496/1‑149 (MQ=255)                                                                                                      cGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTgg                                             <  2:91389/149‑1 (MQ=255)                                                                                                      cGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTgg                                             <  4:39261/149‑1 (MQ=255)                                                                                                      cGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTgg                                             <  4:227475/149‑1 (MQ=255)                                                                                                      cGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTgg                                             <  2:351710/149‑1 (MQ=255)                                                                                                      cGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGATCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATAGCTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTgg                                             <  4:359695/149‑1 (MQ=255)                                                                                                       ggCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATtgtg                                                                                   >  2:459797/1‑111 (MQ=255)                                                                                                       ggCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATtgtg                                                                                   <  1:459797/111‑1 (MQ=255)                                                                                                        gCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGt                                           >  1:345692/1‑149 (MQ=255)                                                                                                              gtgtGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCg                                     >  1:551666/1‑149 (MQ=255)                                                                                                                   cTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACt                                <  3:141920/149‑1 (MQ=255)                                                                                                                       cAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCa                            <  4:168726/149‑1 (MQ=255)                                                                                                                        aGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAg                           >  3:370996/1‑149 (MQ=255)                                                                                                                        aGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAg                           >  2:397158/1‑149 (MQ=255)                                                                                                                               tCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCAt                    <  1:199496/149‑1 (MQ=255)                                                                                                                                     cGTCTGAACCGTGCCTATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTc              <  1:394291/149‑1 (MQ=255)                                                                                                                                     cGTCTGAACCGTGCCTATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACATCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTc              <  1:158298/149‑1 (MQ=255)                                                                                                                                                 gCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAActgct  <  1:397158/149‑1 (MQ=255)                                                                                                                                                   ccATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATcgccgc                                                              >  4:317054/1‑87 (MQ=255)                                                                                                                                                   ccATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATcgccgc                                                              <  3:317054/87‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   CAACCGACGACAGTCCTGCTAAAACGTTCGTTTGATATCATTTTTCCTAAAATTGAATGGCAGAGAATCATGAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCT  >  NC_000913/260319‑260611

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