Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,719,845	G→A	5.5%	R208H (CGC→CAC)	yfiP →	DTW domain protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,719,845	0	G	A	5.5%	237.4 / 3.2	72	R208H (CGC→CAC)	yfiP	DTW domain protein
Reads supporting (aligned to +/- strand):  ref base G (32/36);  new base A (2/2);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

TTTCCGTCGATCTTTCCCGGCTTTCTGCCTATCGCCTGCGTGAAGCCCAGGCTGAAGGCCAATATTGTACTGCCGAGGTAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGAACAGTTAGAAAGCGTTTAAAATCATTCGGTCACTTCTGCGGGAGACCGGTATGAGTCAGCGAGGACTGGAAGCACTACTGCGACCAAAA  >  NC_000913/2719698‑2719991                                                                                                                                                    |                                                                                                                                                   tttCCGTCGATCTTTCCCGGCTTTCTGCCTATCGCCTGCGTGAAGCCCAGGCTGAAGGCCAATATTGTACTGCCGAGGTAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGc                                                                                                                                                   <  3:446872/149‑1 (MQ=255)  ttCCGTCGATCTTTCCCGGCTTTCTGCCTATCGCCTGCGTGAAGCCCAGGCTGAAGGCCAATATTGTACTGCCGAGGTAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCt                                                                                                                                                  >  4:569501/1‑149 (MQ=255)   tCCGTCGATCTTTCCCGGCTTTCTGCCTATCGCCTGCGTGAAGCCCAGGCTGAAGGCCAATATTGTACTGCCGAGGTAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCtt                                                                                                                                                 <  3:240881/149‑1 (MQ=255)     cGTCGATCTTTCCCGGCTTTCTGCCTATCGCCTGCGTGAAGCCCAGGCTGAAGGCCAATATTGTACTGCCGAGGTAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTa                                                                                                                                               >  2:499325/1‑149 (MQ=255)       tCGATCTTTCCCGGCTTTCTGCCTATCGCCTGCGTGAAGCCCAGGCTGAAGGCCAATATTGTACTGCCGAGGTAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCACTTTaaa                                                                                                                                             >  1:589732/1‑149 (MQ=255)        cGATCTTTCCCGGCTTTCTGCCTATCGCCTGCGTGAAGCCCAGGCTGAAGGCCAATATTGTACTGCCGAGGTAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTaaaa                                                                                                                                            <  3:43136/149‑1 (MQ=255)        cGATCTTTCCCGGCTTTCTGCCTATCGCCTGCGTGAAGCCCAGGCTGAAGGCCAATATTGTACTGCCGAGGTAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTaaaa                                                                                                                                            <  1:586638/149‑1 (MQ=255)         gATCTTTCCCGGCTTTCTGCCTATCGCCTGCGTGAAGCCCAGGCTGAAGGCCAATATTGTACTGCCGAGGTAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAac                                                                                                                                           >  2:243443/1‑149 (MQ=255)          aTCTTTCCCGGCTTTCTGCCTATCGCCTGCGTGAAGCCCAGGCTGAAGGCCAATATTGTACTGCCGAGGTAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAaca                                                                                                                                          >  3:416797/1‑149 (MQ=255)          aTCTTTCCCGGCTTTCTGCCTATCGCCTGCGTGAAGCCCAGGCTGAAGGCCAATATTGTACTGCCGAGGTAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAaca                                                                                                                                          >  1:413250/1‑149 (MQ=255)          aTCTTTCCCGGCTTTCTGCCTATCGCCTGCGTGAAGCCCAGGCTGAAGGCCAATATTGTACTGCCGAGGTAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCACTTTAAAaca                                                                                                                                          >  4:200688/1‑149 (MQ=255)           tCTTTCCCGGCTTTCTGCCTATCGCCTGCGTGAAGCCCAGGCTGAAGGCCAATATTGTACTGCCGAGGTAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAacac                                                                                                                                         <  3:42224/149‑1 (MQ=255)             tttCCCGGCTTTCTGCCTATCGCCTGCGTGAAGCCCAGGCTGAAGGCCAATATTGTACTGCCGAGGTAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGc                                                                                                                                       <  4:324862/149‑1 (MQ=255)                 ccGGCTTTCTGCCTATCGCCTGCGTGAAGCCCAGGCTGAAGGCCAATATTGTACTGCCGAGGTAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATc                                                                                                                                   >  3:20957/1‑149 (MQ=255)                  cGGCTTTCTGCCTATCGCCTGCGTGAAGCCCAGGCTGAAGGCCAATATTGTACTGCCGAGGTAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCt                                                                                                                                  >  4:211967/1‑149 (MQ=255)                      tttCTGCCTATCGCCTGCGTGAAGCCCAGGCTGAAGGCCAATATTGTACTGCCGAGGTAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCa                                                                                                                              >  4:441722/1‑149 (MQ=255)                          tGCCTATCGCCTGCGTGAAGCCCAGGCTGAAGGCCAATATTGTACTGCCGAGGTAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGaa                                                                                                                          >  1:4462/1‑149 (MQ=255)                           gCCTATCGCCTGCGTGAAGCCCAGGCTGAAGGCCAATATTGTACTGCCGAGGTAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGaaa                                                                                                                         >  1:433372/1‑149 (MQ=255)                            ccTATCGCCTGCGTGAAGCCCAGGCTGAAGGCCAATATTGTACTGCCGAGGTAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGaaaa                                                                                                                        >  2:580327/1‑149 (MQ=255)                                  gCCTGCGTGAAGCCCAGGCTGAAGGCCAATATTGTACTGCCGAGGTAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCaa                                                                                                                  <  2:385013/149‑1 (MQ=255)                                   ccTGCGTGAAGCCCAGGCTGAAGGCCAATATTGTACTGCCGAGGTAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCt                                                                                                                                                  >  4:503190/1‑116 (MQ=255)                                   ccTGCGTGAAGCCCAGGCTGAAGGCCAATATTGTACTGCCGAGGTAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCt                                                                                                                                                  <  3:503190/116‑1 (MQ=255)                                       cGTGAAGCCCAGGCTGAAGGCCAATATTGTACTGCCGAGGTAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCt                                                                                                             <  1:28710/149‑1 (MQ=255)                                         tGAAGCCCAGGCTGAAGGCCAATATTGTACTGCCGAGGTAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTgg                                                                                                           <  4:264166/149‑1 (MQ=255)                                         tGAAGCCCAGGCTGAAGGCCAATATTGTACTGCCGAGGTAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTgg                                                                                                           <  2:433372/149‑1 (MQ=255)                                         tGAAGCCCAGGCTGAAGGCCAATATTGTACTGCCGAGGTAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTgg                                                                                                           <  4:508335/149‑1 (MQ=255)                                         tGAAGCCCAGGCTGAAGGCCAATATTGTACTGCCGAGGTAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTgg                                                                                                           <  2:62168/149‑1 (MQ=255)                                                cAGGCTGAAGGCCAATATTGTACTGCCGAGGTAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCAt                                                                                                    >  1:512546/1‑149 (MQ=255)                                                    cTGAAGGCCAATATTGTACTGCCGAGGTAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATcaca                                                                                                <  2:391605/149‑1 (MQ=255)                                                    cTGAAGGCCAATATTGTACTGCCGAGGTAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATcaca                                                                                                <  3:301033/149‑1 (MQ=255)                                                    cTGAAGGCCAATATTGTACTGCCGAGGTAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATcaca                                                                                                <  4:22769/149‑1 (MQ=255)                                                      gAAGGACAATATTGTACTGCCGAGGTAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCAcagc                                                                                              <  3:302703/149‑1 (MQ=255)                                                          gCCAATATTGTACTGCCGAGGTAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGaa                                                                                          <  1:437383/149‑1 (MQ=255)                                                                         ccGAGGTAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGAACAGTTAGAAAGCGtt                                                                           >  4:105902/1‑149 (MQ=255)                                                                           gAGGTAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCACTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGAACAGTTAGAAAGCGTTTa                                                                         <  2:589732/149‑1 (MQ=255)                                                                            aGGTAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGAACAGTTAGAAAGCGTTTaa                                                                        >  2:440716/1‑149 (MQ=255)                                                                              gTAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGAACAGTTAGAAAGCGTTTaaaa                                                                      <  3:415967/149‑1 (MQ=255)                                                                               tAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGAACAGTTAGAAAGCGTTTAAAAt                                                                     <  3:569501/149‑1 (MQ=255)                                                                                aGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGAACAGTTAGAAAGCGTTTAAAATc                                                                    <  1:440716/149‑1 (MQ=255)                                                                                 gCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGAACAGTTAGAAAGCGTTTAAAATCa                                                                   >  2:18966/1‑149 (MQ=255)                                                                                   cATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGAACAGTTAGAAAGCGTTTAAAATCAtt                                                                 >  1:500455/1‑149 (MQ=255)                                                                                   cATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGAACAGTTAGAAAGCGTTTAAAATCAtt                                                                 <  4:20957/149‑1 (MQ=255)                                                                                    aTCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAacac                                                                                                                                         >  1:169622/1‑76 (MQ=255)                                                                                    aTCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAacac                                                                                                                                         <  2:169622/76‑1 (MQ=255)                                                                                    aTCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGAACAGTTAGAAAGCGTTTAAAATCATTc                                                                >  3:546545/1‑149 (MQ=255)                                                                                    aTCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGAACAGTTAGAAAGCGTTTAAAATCATTc                                                                >  2:1987/1‑149 (MQ=255)                                                                                     tCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGAACAGTTAGAAAGCGTTTAAAATCATTCg                                                               >  1:71204/1‑149 (MQ=255)                                                                                     tCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCACTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGAACAGTTAGAAAGCGTTTAAAATCATTCg                                                               <  3:200688/149‑1 (MQ=255)                                                                                           tGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGAACAGTTAGAAAGCGTTTAAAATCATTCGGTCACt                                                         <  3:345096/149‑1 (MQ=255)                                                                                           tGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGAACAGTTAGAAAGCGTTTAAAATCATTCGGTCACt                                                         <  3:105902/149‑1 (MQ=255)                                                                                             ttAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGAACAGTTAGAAAGCGTTTAAAATCATTCGGTCACTTc                                                       >  3:400293/1‑149 (MQ=255)                                                                                             ttAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGAACAGTTAGAAAGCGTTTAAAATCATTCGGTCACTTc                                                       >  3:105137/1‑149 (MQ=255)                                                                                                 atatGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGAACAGTTAGAAAGCGTTTAAAATCATTCGGTCACTTCTGCg                                                   >  2:362178/1‑149 (MQ=255)                                                                                                            gATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGAACAGTTAGAAAGCGTTTAAAATCATTCGGTCACTTCTGCGGGAGACCGGTa                                        <  1:189425/149‑1 (MQ=255)                                                                                                            gATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGAACAGTTAGAAAGCGTTTAAAATCATTCGGTCACTTCTGCGGGAGACCGGTa                                        <  3:133645/149‑1 (MQ=255)                                                                                                            gATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGAACAGTTAGAAAGCGTTTAAAATCATTCGGTCACTTCTGCGGGAGACCGGTa                                        <  4:105137/149‑1 (MQ=255)                                                                                                               aCCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGAACAGTTAGAAAGCGTTTAAAATCATTCGGTCACTTCTGCGGGAGACCGGTATGa                                     <  1:580327/149‑1 (MQ=255)                                                                                                                   gggCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGAACAGTTAGAAAGCGTTTAAAATCATTCGGTCACTTCTGCGGGAGACCGGTATGAGTCa                                 <  4:543682/149‑1 (MQ=255)                                                                                                                     gcggcggcAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGAACAGTTAGAAAGCGTTTAAAATCATTCGGTCACTTCTGCGGGAGACCGGTATGAGTCAGc                               <  2:500455/149‑1 (MQ=255)                                                                                                                          ggcAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGAACAGTTAGAAAGCGTTTAAAATCATTCGGTCACTTCTGCGGGAGACCGGTATGAGTCAGCGAGGa                          >  4:10165/1‑149 (MQ=255)                                                                                                                             aGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGAACAGTTAGAAAGCGTTTAAAATCATTCGGTCACTTCTGCGGGAGACCGGTATGa                                     <  2:281374/135‑1 (MQ=255)                                                                                                                             aGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGAACAGTTAGAAAGCGTTTAAAATCATTCGGTCACTTCTGCGGGAGACCGGTATGa                                     >  1:281374/1‑135 (MQ=255)                                                                                                                                    ggCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGAACAGTTAGAAAGCGTTTAAAATCATTCGGTCACTTCTGCGGGAGACCGGTATGAGTCAGCGAGGACTGGAAGCac                <  4:433744/149‑1 (MQ=255)                                                                                                                                     gCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGAACAGTTAGAAAGCGTTTAAAATCATTCGGTCACTTCTGCGGGAGACCGGTATGAGTCAGCGAGGACTGGAAGCact               >  2:440848/1‑149 (MQ=255)                                                                                                                                     gCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGAACAGTTAGAAAGCGTTTAAAATCATTCGGTCACTTCTGCGGGAGACCGGTATGAGTCAGCGAGGACTGGAAGCact               >  3:467171/1‑149 (MQ=255)                                                                                                                                      cGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGAACAGTTAGAAAGCGTTTAAAATCATTCGGTCACTTCTGCGGGAGACCGGTATGAGTCAGCGAGGACTGGAAGCacta              <  3:201936/149‑1 (MQ=255)                                                                                                                                          cATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGAACAGTTAGAAAGCGTTTAAAATCATTCGGTCACTTCTGCGGGAGACCGGTATGAGTCAGCGa                             <  3:23315/130‑1 (MQ=255)                                                                                                                                          cATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGAACAGTTAGAAAGCGTTTAAAATCATTCGGTCACTTCTGCGGGAGACCGGTATGAGTCAGCGa                             >  4:23315/1‑130 (MQ=255)                                                                                                                                            ttttACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGAACAGTTAGAAAGCGTTTAAAATCATTCGGTCACTTCTGCGGGAGACCGGTATGAGTCAGCGAGGACTGGAAGCACTACTGCGa        >  4:470140/1‑149 (MQ=255)                                                                                                                                                 cccGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGAACAGTTAGAAAGCGTTTAAAATCATTCGGTCACTTCTGCGGGAGACCGGTATGAGTCAGCGAGGACTGGAAGCACTACTGCGACCaaa   >  1:339871/1‑149 (MQ=255)                                                                                                                                                 cccGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGAACAGTTAGAAAGCGTTTAAAATCATTCGGTCACTTCTGCGGGAGACCGGTATGAGTCAGCGAGGACTGGAAGCACTACTGCGACCaaa   <  3:340301/149‑1 (MQ=255)                                                                                                                                                  ccGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGAACAGTTAGAAAGCGTTTAAAATCATTCGGTCACTTCTGCGGGAGACCGGTATGAGTCAGCGAGGACTGGAAGCACTACTGCGACCaaaa  >  3:194232/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                   TTTCCGTCGATCTTTCCCGGCTTTCTGCCTATCGCCTGCGTGAAGCCCAGGCTGAAGGCCAATATTGTACTGCCGAGGTAGCCATCGCACTGTTAGATATGGCCGGCGATACCGGGGCGGCGGCAGGGTTAGGCGAGCATTTTACCCGCTTTAAAACACGCTATCTGGCAGGAAAAACGCAACATCTGGGTAGCATCACAGCAGAACAGTTAGAAAGCGTTTAAAATCATTCGGTCACTTCTGCGGGAGACCGGTATGAGTCAGCGAGGACTGGAAGCACTACTGCGACCAAAA  >  NC_000913/2719698‑2719991

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: ‑