Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	344,938	G→T	6.0%	G255C (GGT→TGT)	yahL →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	344,938	0	G	T	6.0%	238.4 / 2.8	67	G255C (GGT→TGT)	yahL	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (32/31);  new base T (2/2);  total (34/33)
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.34e-01

TATTACACTTAAAAAATATCTACCCTGTTTTTCCTGAATCACTGGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCATCTGTATATCTGTAT  >  NC_000913/344804‑345076                                                                                                                                       |                                                                                                                                           tattaCACTTAAAAAATATCTACCCTGTTTTTCCTGAATCTCTGGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGAc                                                                                                                                       <  2:771791/140‑1 (MQ=255) tattaCACTTAAAAAATATCTACCCTGTTTTTCCTGAATCACTGGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCTATTAAATGAATGGATATATAATAGAGTGCAGGAATTAGGTGAc                                                                                                                                       <  1:9963/140‑1 (MQ=255) tattaCACTTAAAAAATATCTACCCTGTTTTTCCTGAATCACTGGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGAc                                                                                                                                       <  1:714729/140‑1 (MQ=255) tattaCACTTAAAAAATATCTACCCTGTTTTTCCTGAATCACTGGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGAc                                                                                                                                       <  2:305318/140‑1 (MQ=255) tattaCACTTAAAAAATATCTACCCTGTTTTTCCTGAATCACTGGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGAc                                                                                                                                       <  1:13129/140‑1 (MQ=255) tattaCACTTAAAAAATATCTACCCTGTTTTTCCTGAATCACTGGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGAc                                                                                                                                       <  1:80482/140‑1 (MQ=255)      cacTTAAAAAATATCTACCCTGTTTTTCCTGAATCACTGGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTATTAGAGTGCAGGAATTAGGTGACAGata                                                                                                                                  >  1:1211613/1‑140 (MQ=255)      cacTTAAAAAATATCTACCCTGTTTTTCCTGAATCACTGGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGata                                                                                                                                  >  1:1034589/1‑140 (MQ=255)      cacTTAAAAAATATCTACCCTGTTTTTCCTGAATCACTGGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGata                                                                                                                                  >  1:468370/1‑140 (MQ=255)           aaaaaaTATCTACCCTGTTTTTCCTGAATCACTGGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTATGTGACAGATATTATc                                                                                                                             >  1:1078281/1‑140 (MQ=255)           aaaaaaTATCTACCCTGTTTTTCCTGAATCACTGGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTATGTGACAGATATTATc                                                                                                                             >  1:1078239/1‑140 (MQ=255)           aaaaaaTATCTACCCTGTTTTTCCTGAATCACTGGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATc                                                                                                                             >  1:236441/1‑140 (MQ=255)            aaaaaTATCTACCCTGTTTTTCCTGAATCACTGGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATcc                                                                                                                            <  2:201631/140‑1 (MQ=255)              aaaTATCTACCCTGTTTTTCCTGAATCACTGGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTg                                                                                                                          <  2:896784/140‑1 (MQ=255)                   tCTACCCTGTTTTTCCTGAATCACTGGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAAt                                                                                                                     >  1:171014/1‑140 (MQ=255)                   tCTACCCTGTTTTTCCTGAATCACTGGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAAt                                                                                                                     >  2:426559/1‑140 (MQ=255)                                 ccTGAATCACTGGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAAt                                                                                                       >  1:927688/1‑140 (MQ=255)                                      aTCACTGGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCt                                                                                                  <  1:650474/140‑1 (MQ=255)                                      aTCACTGGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCt                                                                                                  <  1:650496/140‑1 (MQ=255)                                       tCACTGGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTa                                                                                                 <  2:738501/140‑1 (MQ=255)                                       tCACTGGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTa                                                                                                 <  2:740153/140‑1 (MQ=255)                                         aCTGGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTaca                                                                                               <  1:728804/140‑1 (MQ=255)                                                       tGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGa                                                                                 <  2:591605/140‑1 (MQ=255)                                                           ctctGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGAAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCATAAAgg                                                                             >  1:522919/1‑140 (MQ=255)                                                             ctGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAGTTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGt                                                                           >  1:5197/1‑140 (MQ=255)                                                              tGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGata                                                                                                                                  <  1:693379/84‑1 (MQ=255)                                                              tGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGata                                                                                                                                  >  2:693379/1‑84 (MQ=255)                                                                   tGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCAcc                                                                     >  1:184433/1‑140 (MQ=255)                                                                   tGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCAcc                                                                     >  1:819339/1‑140 (MQ=255)                                                                          aaaaaCATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAgt                                                              >  2:992433/1‑140 (MQ=255)                                                                               cATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTa                                                         >  2:985569/1‑140 (MQ=255)                                                                                 tACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTAtt                                                       <  1:947165/140‑1 (MQ=255)                                                                                 tACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTAtt                                                       <  1:348391/140‑1 (MQ=255)                                                                                      agaTGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCAcc                                                  >  2:1072245/1‑140 (MQ=255)                                                                                      agaTGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCAcc                                                  >  2:1072337/1‑140 (MQ=255)                                                                                         tgatgCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTga                                               <  2:687348/140‑1 (MQ=255)                                                                                                 gATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAAt                                       <  1:796789/140‑1 (MQ=255)                                                                                                  aTTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATc                                      >  1:54701/1‑140 (MQ=255)                                                                                                   ttAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCg                                     <  2:557196/140‑1 (MQ=255)                                                                                                      aatgaatgGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCaa                                  <  1:810868/140‑1 (MQ=255)                                                                                                              gATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTACAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTa                          <  2:470344/140‑1 (MQ=255)                                                                                                                tatCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAAt                        <  1:490753/140‑1 (MQ=255)                                                                                                                      aTAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTc                  <  2:1034589/140‑1 (MQ=255)                                                                                                                      aTAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTc                  >  1:796620/1‑140 (MQ=255)                                                                                                                      aTAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTc                  >  1:796615/1‑140 (MQ=255)                                                                                                                      aTAGAGTGCAGGAATTAGGTGACAGATATTAACCTGTAAAAAAGCACGTTAAAATCCGCTACACACTTTAGTTCAGACAGGGATCCACCTGCAAGTTGTTAGTCCACCTCAGACTGACTCGCAATGCAGGTAATCATTTc                  <  2:832501/140‑1 (MQ=255)                                                                                                                      aTAGAGTGCAGGAAGTAGGTGACAGTTATTATCCTGTAAATAAGCACGTTAAACTCCGCGACTCTCTTTAGTTCAGAATGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTc                  <  2:499066/140‑1 (MQ=255)                                                                                                                        agagTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCat                <  2:839259/140‑1 (MQ=255)                                                                                                                        agagTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTAGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCat                <  2:236407/140‑1 (MQ=255)                                                                                                                         gagTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCAACCTGAGAGTGAATAGCAATGCAGGTAATCGTTTCatc               >  1:143509/1‑140 (MQ=255)                                                                                                                          agTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCatct              >  2:853973/1‑140 (MQ=255)                                                                                                                          agTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCAGCTACACTCTTTAGTTCAGTAAGGGTTCCACCTGCAAGTTGTTATTCCAGCTGAGAGTGAATCGCAATGCAGGTAATCATTTCatct              >  2:141546/1‑140 (MQ=255)                                                                                                                            tGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCatctgt            >  2:273944/1‑140 (MQ=255)                                                                                                                            tGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCatctgt            >  1:926852/1‑140 (MQ=255)                                                                                                                            tGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCatctgt            >  2:1050295/1‑140 (MQ=255)                                                                                                                            tGCAGGAATTAGGTGACAGATATAATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAATGGGTTCCACCTGAAATTTTTTATTCCACCTGAGAGTGAATCTCTATACATGTAgtcatgtcataggt            >  2:1231025/1‑136 (MQ=255)                                                                                                                             gCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCatctgta           >  2:524189/1‑140 (MQ=255)                                                                                                                             gCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCatctgta           >  2:524218/1‑140 (MQ=255)                                                                                                                             gCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCatctgta           <  1:1175498/140‑1 (MQ=255)                                                                                                                              cAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCatctgtat          >  1:124147/1‑140 (MQ=255)                                                                                                                                ggAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCatctgtatat        >  1:428168/1‑140 (MQ=255)                                                                                                                                 gAATTATGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCatctgtatatc       <  2:1078281/140‑1 (MQ=255)                                                                                                                                 gAATTATGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCatctgtatatc       <  2:1078239/140‑1 (MQ=255)                                                                                                                                 gAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCatctgtatatc       <  2:1005748/140‑1 (MQ=255)                                                                                                                                 gAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCatctgtatatc       <  2:1005973/140‑1 (MQ=255)                                                                                                                                    ttAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCatctgtatatctgt    >  1:1036130/1‑140 (MQ=255)                                                                                                                                      aGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCatctgtatatctgtat  >  2:389147/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                                           TATTACACTTAAAAAATATCTACCCTGTTTTTCCTGAATCACTGGTGATGTTATTGCTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCATCTGTATATCTGTAT  >  NC_000913/344804‑345076

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