Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	914,569	T→G	K97Q (AAA→CAA)	ybjE ←	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	914,569	0	T	G	98.3%	216.3 / ‑2.8	60	K97Q (AAA→CAA)	ybjE	putative transporter
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (31/28);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

GGAAAAGTAGCAAAATTAACGTGTATTCACTGGCTTCGGTCGCGTGTTGTAAGAAAGCCAGTCCACTTAGACCAATGGCAAAACCAATCACTACTACGCCGCACAGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCAACAGGTTACTGGCGAGGTTATCGAGAAACGCCAG  >  NC_000913/914422‑914716                                                                                                                                                     |                                                                                                                                                    ggAAAAGTAGCAAAATTAACGTGTATTCACTGGCTTCGGTCGCGTGTTGTAAGAAAGCCAGTCCACTTAGACCAATGGCAAAACCAATCACTACTACGCCGCACAGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTc                                                                                                                                                  <  1:794355/151‑1 (MQ=255)     aaGTAGCAAAATTAACGTGTATTCACTGGCTTCGGTCGCGTGTTGTAAGAAAGCCAGTCCACTTAGACCAATGGCAAAACCAATCACTACTACGCCGCACAGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGc                                                                                                                                              >  1:400316/1‑151 (MQ=255)      aGTAGCAAAATTAACGTGTATTCACTGGCTTCGGTCGCGTGTTGTAAGAAAGCCAGTCCACTTAGACCAATGGCAAAACCAATCACTACTACGCCGCACAGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGc                                                                                                                                              >  1:478393/1‑150 (MQ=255)      aGTAGCAAAATTAACGTGTATTCACTGGCTTCGGTCGCGTGTTGTAAGAAAGCCAGTCCACTTAGACCAATGGCAAAACCAATCACTACTACGCCGCACAGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCt                                                                                                                                             <  1:581151/151‑1 (MQ=255)                 tAACGTGTATTCACTGGCTTCGGTCGCGTGTTGTAAGAAAGCCAGTCCACTTAGACCAATGGCAAAACCAATCACTACTACGCCGCACAGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGtt                                                                                                                                   <  2:746679/150‑1 (MQ=255)                 tAACGTGTATTCACTGGCTTCGGTCGCGTGTTGTAAGAAAGCCAGTCCACTTAGACCAATGGCAAAACCAATCACTACTACGCCGCACAGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGtt                                                                                                                                   <  1:803065/150‑1 (MQ=255)                    cGTGTATTCACTGGCTTCGGTCGCGTGTTGTAAGAAAGCCAGTCCACTTAGCCCAATGGCAAAACCAATCACTACTACGCCGCACAGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTgcg                                                                                                                                <  2:610518/150‑1 (MQ=255)                    cGTGTATTCACTGGCTTCGGTCGCGTGTTGTAAGAAAGCCAGTCCACTTAGACCAATGGCAAAACCAATCACTACTACGCCGCACAGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTgcgc                                                                                                                               <  1:637991/151‑1 (MQ=255)                       gtATTCACTGGCTTCGGTCGCGTGTTGTAAGAAAGCCAGTCCACTTAGACCAATGGCAAAACCAATCACTACTACGCCGCACAGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCAc                                                                                                                            >  1:5214/1‑151 (MQ=255)                                    tCGGTCGCGTGTTGTAAGAAAGCCAGTCCACTTAGACCAATGGCAAAACCAATCACTACTACGCCGCACAGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGc                                                                                                                >  2:595166/1‑150 (MQ=255)                                      ggTCGCGTGTTGTAAGAAAGCCAGTCCACTTAGACCAATGGCAAAACCAATCACTACTACGCCGCACAGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTcc                                                                                                             >  2:148117/1‑151 (MQ=255)                                      ggTCGCGTGTTGTAAGAAAGCCAGTCCACTTAGACCAATGGCAAAACCAATCACTACTACGCCGCACAGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTcc                                                                                                             >  2:566364/1‑151 (MQ=255)                                         cgcgTGTTGTAAGAAAGCCAGTCCACTTAGACCAATGGCAAAACCAATCACTACTACGCCGCACAGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGc                                                                                                          <  2:667779/151‑1 (MQ=255)                                         cgcgTGTTGTAAGAAAGCCAGTCCACTTAGACCAATGGCAAAACCAATCACTACTACGCCGCACAGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGc                                                                                                          <  1:506201/151‑1 (MQ=255)                                          gcgTGTTGTAAGAAAGCCAGTCCACTTAGACCAATGGCAAAACCAATCACTACTACGCCGCACAGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGcc                                                                                                         <  2:898978/151‑1 (MQ=255)                                            gtgtTGTAAGAAAGCCAGTCCACTTAGACCAATGGCAAAACCAATCACTACTACGCCGCACAGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCcac                                                                                                       <  1:478156/151‑1 (MQ=255)                                                 gtAAGAAAGCTAGTCCACTTAGACCAATGGAAAAACCAATCAATACTACGCCGTACAGTTTTAGCGACTCCA‑GCGCCATCCCAATACGCGACGGGAGTTGATCTTGCTGATGGTGGCTGCGCCAGGGCAGGCCTCGCTCCAGCCACAtaa                                                                                                   >  1:731094/1‑148 (MQ=255)                                                 gtAAGAAAGCCAGTCCACTTAGACCAATGGCAAAACCAATCACTACTACGCCGCACAGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCa                                                                                                   <  2:342605/150‑1 (MQ=255)                                                    aGAAAGCCAGTCCACTTAGACCAATGGCAAAACCAATCACTACTACGCCGCACAGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGggc                                                                                               <  1:531616/151‑1 (MQ=255)                                                      aaaGCCAGTCCACTTAGACCAATGGCAAAACCAATCACTACTACGCCGCACAGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGggcg                                                                                              <  2:583148/150‑1 (MQ=255)                                                       aaGCCAGTCCACTTAGACCAATGGCAAAACCAATCACTACTACGCCGCACAGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGggcggc                                                                                            <  2:375149/151‑1 (MQ=255)                                                               ccACTTAGACCAATGGCAAAACCAATCACTACTACGCCGCACAGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTac                                                                                    <  2:134471/151‑1 (MQ=255)                                                                cACTTAGACCAATGGCAAAACCAATCACTACTACGCCGCACAGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTaca                                                                                   <  1:687725/151‑1 (MQ=255)                                                                 aCTTAGACCAATGGCAAAACCAATCACTACTACGCCGCACAGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACCTCAGGGCGGCAATATTaca                                                                                   >  1:561867/1‑150 (MQ=255)                                                                       aCCAATGGCAAAACCAATCACTACTACGCCGCACAGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTaa                                                                             >  1:633535/1‑150 (MQ=255)                                                                        ccAATGGCAAAACCAATCACTACTACGCCGCACAGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTaaa                                                                            >  1:280385/1‑150 (MQ=255)                                                                         cAATGGCAAAACCAATCACTACTACGCCGCACAGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAAt                                                                          >  2:671442/1‑151 (MQ=255)                                                                         cAATGGCAAAACCAATCACTACTACGCCGCACAGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAAt                                                                          <  1:111045/151‑1 (MQ=255)                                                                          aaTGGCAAAACCAATCACTACTACGCCGCACAGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATa                                                                         <  2:385521/151‑1 (MQ=255)                                                                                    ccAATCACTACTACGCCGCACAGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATAc                                                               <  2:73377/151‑1 (MQ=255)                                                                                          actactACGCCGCACAGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACgg                                                         >  2:169026/1‑151 (MQ=255)                                                                                           ctactaCGCCGCACAGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACgg                                                         <  1:663806/150‑1 (MQ=255)                                                                                             actacGCCGCACAGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGc                                                        >  2:704678/1‑149 (MQ=255)                                                                                               tacGCCGCACAGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGTGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGa                                                     >  2:62059/1‑150 (MQ=255)                                                                                                  gccgcACAGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGaataa                                                 <  2:789938/151‑1 (MQ=255)                                                                                                   ccgcACAGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGaataat                                                >  1:209412/1‑151 (MQ=255)                                                                                                   ccgcACAGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGaataat                                                >  2:824707/1‑151 (MQ=255)                                                                                                    cgcACAGGTTTTGCCGATCCAGGCGCCATC‑CAATACGGGACGGGGGGTGTTCTTGCTGAGGGGGGTTTCCTCCCGGGAGGCCTTGGTTCAGCCACATCAGGGGGGCCATATTACCCAGGAAAATAACGGTAATACTGACGGCAGaataa                                                 <  1:655456/149‑1 (MQ=255)                                                                                                     gcACAGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGc                                              >  2:535397/1‑151 (MQ=255)                                                                                                       acaGTTCTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCa                                             <  1:131202/150‑1 (MQ=255)                                                                                                         aGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGa                                           >  2:607787/1‑150 (MQ=255)                                                                                                          gTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAAt                                         >  2:19281/1‑151 (MQ=255)                                                                                                            tttAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATc                                        >  2:307479/1‑150 (MQ=255)                                                                                                                gCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCa                                    >  1:448444/1‑150 (MQ=255)                                                                                                                   aCTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCAAc                                  >  1:754723/1‑149 (MQ=255)                                                                                                                   aCTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCAACAg                                <  2:96341/151‑1 (MQ=255)                                                                                                                    cTCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCAACAgg                               >  1:858611/1‑151 (MQ=255)                                                                                                                     tCCA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCAACAGGt                              >  2:402196/1‑151 (MQ=255)                                                                                                                      ccA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCAACAgg                               <  1:351295/149‑1 (MQ=255)                                                                                                                      ccA‑GCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCAACAGGtt                             >  1:685214/1‑151 (MQ=255)                                                                                                                           cgcCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCAACAGGTTACTg                         >  1:662973/1‑151 (MQ=255)                                                                                                                           cgcCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCAACAGGTTACTg                         >  2:864871/1‑151 (MQ=255)                                                                                                                             ccATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCGCGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCAACAGGTTACTGGc                       <  2:68372/151‑1 (MQ=255)                                                                                                                             ccATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCAACAGGTTACTGGc                       >  1:237619/1‑151 (MQ=255)                                                                                                                                     aTACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCAACAGGTTACTGGCGAGGTTAt               <  1:107181/151‑1 (MQ=255)                                                                                                                                      tACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCAACAGGTTACTGGCGAGGTTATc              >  1:59352/1‑151 (MQ=255)                                                                                                                                        cgcgACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCAACAGGTTACTGGCGAGGTTATCga            <  1:483693/151‑1 (MQ=255)                                                                                                                                           gACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCAACAGGTTACTGGCGAGGTTATCgaga          >  1:890284/1‑150 (MQ=255)                                                                                                                                            aCGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCAACAGGTTACTGGCGAGGTTATCGAGaa         >  1:284186/1‑150 (MQ=255)                                                                                                                                             cGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCAACAGGTTACTGGCGAGGTTATCGAGaaa        <  1:624503/150‑1 (MQ=255)                                                                                                                                                  gttgttCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCAACAGGTTACTGGCGAGGTTATCGAGAAACGCCAg  <  2:482208/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GGAAAAGTAGCAAAATTAACGTGTATTCACTGGCTTCGGTCGCGTGTTGTAAGAAAGCCAGTCCACTTAGACCAATGGCAAAACCAATCACTACTACGCCGCACAGTTTTAGCGACTCCA‑GCGCCATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCAACAGGTTACTGGCGAGGTTATCGAGAAACGCCAG  >  NC_000913/914422‑914716

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