Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,632,711	G→A	R454C (CGT→TGT)	guaB ←	inosine 5'‑monophosphate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,632,711	0	G	A	95.6%	125.2 / ‑2.2	46	R454C (CGT→TGT)	guaB	inosine 5'‑monophosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (24/19);  total (26/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.13e-01

TGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCG  >  NC_000913/2632474‑2632949                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                               tGAGTGTACTGAGAACCGAAGTCCAGAATGAGGCTGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTtcgtcgcg                                                                                                                                                                                                                                    >  1:276271/1‑248 (MQ=255) tGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAAGGTTTTCCGTCATTGACGCGTATTCCGAGGCAAGTGAAACAGATAATATAAAGCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTCCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTtcgtcg                                                                                                                                                                                                                                      <  2:276271/248‑1 (MQ=255)                        cAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCCAGTGCAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTAGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATCCGTACACACTCCGCTTTAGTACACATTTCGTCGATAGTACCACTACCGGTCAGACCCAt                                                                                                                                                                                                            >  2:209743/1‑251 (MQ=255)                          gAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATCCGTACAAACTCCGCTTTAGTCCACAGTTCGTCGCAAGTACCACAGCCGGTCCGACCCATa                                                                                                                                                                                                           >  2:339991/1‑250 (MQ=255)                                       atgcTTATGCATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAAATCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGAGCAGgcc                                                                                                                                                                                             >  1:471472/1‑251 (MQ=255)                                              tGAATGTTTTCCGTCATTGACGCTTATTCCGCGGCAAGCGAAACAGGTAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCACAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCATGCCGCCCATc                                                                                                                                                                                      >  1:145240/1‑251 (MQ=255)                                                aaTGGTTTCCGTCATTGACGCTTATTCCTGGGCAAGTGAAACAGATAATATAAAACGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGTCCAGACGGTAGTTCGGGGACTCTTTAGGAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTCCCACAGCCGGTCAGACCCATACAGGAGTGCAGGCCGCCCATctg                                                                                                                                                                                    <  2:491369/251‑1 (MQ=255)                                                           ccATTGACGCTGAGTCCGAGGCAAGTTAAACAGATAATATAAATCCCCCGCCAGGAAGTCGGGCGAAGAGAGTCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCGTGAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCATGCCGCCCATCTGCTGGTGAATGa                                                                                                                                                                         <  2:144368/250‑1 (MQ=255)                                                             aTTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACACACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGCCCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGAtc                                                                                                                                                                       >  1:330912/1‑251 (MQ=255)                                                             aTTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCAGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGAtc                                                                                                                                                                       >  1:141190/1‑251 (MQ=255)                                                                   gCTTATTCCGAGGCAAGAGAAACAGATAATATATATCGGCCGACATGAAGTCGGGCGAGGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGCACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTc                                                                                                                                                                 >  2:58935/1‑251 (MQ=255)                                                                         tCCGAGGCAAGGGAAACAGATAATATAATACGCCCGACATGAAGTCGGGCGAAGAGAAGCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGGCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACGCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCAGCTGCTGGTGAATGATCTCTTTCAGGCg                                                                                                                                                            <  2:296762/250‑1 (MQ=255)                                                                              ggCAAGTGAAAGAGATAATATAAATCGCCCGACATGAAGACGGGCGAAGCGAACCAGGAGCACAGACGGAACTTCTGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCCACGCTGGTACGTAACAACTCCGCTTTAGGACCCAGTTCGTCGATCGTACCACAGCTGGTCAGACCCCTACAGGAGCGCAGGACGCCCCTATGCCGGTTAATGATCTCTTTCacgccactttt                                                                                                                                                      >  2:12035/1‑241 (MQ=255)                                                                                    ggAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGc                                                                                                                                                <  1:169403/250‑1 (MQ=255)                                                                                     gAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCt                                                                                                                                               >  1:8932/1‑251 (MQ=255)                                                                                        aCAGATAATATAAATCGCCCGAAATGAAGTCGGGCGAAAAGAATCAGGAGCCCAGCCGGTATTTCGGTGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACACCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAAGGATCTCTTTCAGGCGACCTTTATAGGCTAcg                                                                                                                                            <  2:252393/251‑1 (MQ=255)                                                                                         cAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTAcgc                                                                                                                                           >  1:358288/1‑251 (MQ=255)                                                                                               atatAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCtt                                                                                                                                     <  1:318206/251‑1 (MQ=255)                                                                                                 ataAATCGGCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTACCACAGCCGTGCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTc                                                                                                                                    <  1:58935/250‑1 (MQ=255)                                                                                                    gctCGCCCGACATGGAGTCTGGAGAAGAGAATCAGGAGCCCAGAAGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGCACGTGGCTTTCCTGAATGACCGCACCGCTGATACGTCCAAACTCCGCTTTAGTACACAGTTCGTCGCGAGTACCACAGCCGGGCAGACCCATAAAGGAGCGCGGGCCGCCCCGCTGCTGGTGAACGATCTCTTTCACACGTCATTTATCGGCAACGCGACCTTCggtc                                                                                                                                >  2:193526/3‑248 (MQ=255)                                                                                                        gCCCGACATGAAGTCGGGCGAAGATAATCAGGAGCCCAGACGGTATTTCGGGGACTCTTGAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTGCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATCCCtt                                                                                                                            <  2:271593/251‑1 (MQ=255)                                                                                                                       gggCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCttt                                                                                                                                                      <  2:393846/210‑1 (MQ=255)                                                                                                                       gggCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCttt                                                                                                                                                      >  1:393846/1‑210 (MQ=255)                                                                                                                             agagAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGCGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGcgg                                                                                                       >  2:277004/1‑249 (MQ=255)                                                                                                                                gaATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGt                                                                                                    >  1:188484/1‑251 (MQ=255)                                                                                                                                    cAGGAGCCCAGACGGTAGTTCGGGGACTCGTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGTCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATc                                                                                                <  1:461184/251‑1 (MQ=255)                                                                                                                                           ccAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGAc                                                                                                                                                                                                                >  1:84790/1‑132 (MQ=255)                                                                                                                                           ccAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGAc                                                                                                                                                                                                                <  2:84790/132‑1 (MQ=255)                                                                                                                                                gtgtcgtccGGGGTCTCTTTAGTAATGGTAACGGCGTTAACGTGGCTTTCCTGAATGCCCGTACCGCTGATACGTACAAATTCCGCTTTAGTACACAGTTCGTCGATAGTACTACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATa                                                                                    <  2:466218/243‑1 (MQ=255)                                                                                                                                                 ggTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATa                                                                                    >  1:175083/1‑250 (MQ=255)                                                                                                                                                          gggACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAAGGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCagag                                                                          <  1:97366/251‑1 (MQ=255)                                                                                                                                                                    aGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTgg                                                                >  1:478652/1‑251 (MQ=255)                                                                                                                                                                          ggTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGt                                                                                                                                                                                                                                 <  2:252757/84‑1 (MQ=255)                                                                                                                                                                          ggTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGt                                                                                                                                                                                                                                 >  1:252757/1‑84 (MQ=255)                                                                                                                                                                                        cGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAAccc                                            >  2:102710/1‑251 (MQ=255)                                                                                                                                                                                           ggCTTTCCTGAATGCCCGCCCCGCTGATACGTACAACCTCCGCTTTAGTACACAGTTCGGCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATa                                         <  2:8932/251‑1 (MQ=255)                                                                                                                                                                                                     aaTGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAg                                                                                                                                                                                                >  1:429171/1‑90 (MQ=255)                                                                                                                                                                                                     aaTGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAg                                                                                                                                                                                                <  2:429171/90‑1 (MQ=255)                                                                                                                                                                                                         cccGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGCAAGAttt                           >  1:172686/1‑251 (MQ=255)                                                                                                                                                                                                              aCCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTa                        >  1:448730/1‑249 (MQ=255)                                                                                                                                                                                                                       aCGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCTCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCt             <  1:293487/251‑1 (MQ=255)                                                                                                                                                                                                                        cGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTc                                                                                                            >  2:101867/1‑155 (MQ=255)                                                                                                                                                                                                                        cGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTGCCGGCACCAGTTTGTc                                                                                                            <  1:101867/155‑1 (MQ=255)                                                                                                                                                                                                                             aaaCTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGACCATCGCGCCCAAGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTgga          >  2:359986/1‑247 (MQ=255)                                                                                                                                                                                                                              aaCTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATCCCTTCCGGCACCAGTTTGTCGGCAGAGTTCTCGCTCTGGAAATACCGGTCAGAGGAACCTTTGGACATCGCGCCCCGGGAACCCATACCACGGTAAGATTTTTAAGAACGGCCCTGGttacg      >  2:234394/1‑247 (MQ=255)                                                                                                                                                                                                                                  ccGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCg  <  1:393191/251‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                               TGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCG  >  NC_000913/2632474‑2632949

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