Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,059,758	A→G	N168S (AAC→AGC)	torA →	trimethylamine N‑oxide reductase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,059,758	0	A	G	100.0%	222.5 / NA	63	N168S (AAC→AGC)	torA	trimethylamine N‑oxide reductase 1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (30/33);  total (30/33)

GATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGGAGATTACTCTACCGGTGCTGCGCAGGTGATCCTGCCGCGCGTAGTCGGTTCGATGGAAGTGTATGAACAGCA  >  NC_000913/1059622‑1059890                                                                                                                                         |                                                                                                                                     gATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAgcgc                                                                                                                                   <  1:441605/140‑1 (MQ=255) gATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAgcgc                                                                                                                                   <  1:326774/140‑1 (MQ=255)      ccGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCgg                                                                                                                              <  2:799908/140‑1 (MQ=255)        gTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGTAATCGACGGGGATGTTCCATAGCGCTTCgggg                                                                                                                            <  1:863806/140‑1 (MQ=255)          tttGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGAt                                                                                                                          >  2:77701/1‑140 (MQ=255)          tttGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGAt                                                                                                                          >  1:838275/1‑140 (MQ=255)          tttGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGAt                                                                                                                          >  2:429004/1‑140 (MQ=255)            tgtgCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGc                                                                                                                        >  1:471171/1‑140 (MQ=255)                        cTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTa                                                                                                            <  1:291179/140‑1 (MQ=255)                          gggATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTAtt                                                                                                          <  1:77701/140‑1 (MQ=255)                            gATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGc                                                                                                        >  2:84627/1‑140 (MQ=255)                             aTGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGcc                                                                                                       <  1:84627/140‑1 (MQ=255)                                aaGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTg                                                                                                    <  1:242529/140‑1 (MQ=255)                                 aGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGcc                                                                                                       >  2:584418/1‑136 (MQ=255)                                 aGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGcc                                                                                                       <  1:584418/136‑1 (MQ=255)                                 aGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGc                                                                                                   >  1:643804/1‑140 (MQ=255)                                  gCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCa                                                                                                  <  2:747939/140‑1 (MQ=255)                                  gCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCa                                                                                                  <  2:756952/140‑1 (MQ=255)                                  gCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCa                                                                                                  >  1:21461/1‑140 (MQ=255)                                     cTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTa                                                                                                            <  2:121974/127‑1 (MQ=255)                                     cTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTa                                                                                                            >  1:121974/1‑127 (MQ=255)                                         aCATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAAt                                                                                           <  1:597230/140‑1 (MQ=255)                                             gTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCg                                                                                       >  2:776423/1‑140 (MQ=255)                                              ttCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGACGAAAGCTATTG‑CTTGCATGGTaa                                                                                            <  1:3210/133‑1 (MQ=255)                                              ttCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGACGAAAGCTATTG‑CTTGCATGGTaa                                                                                            >  2:3210/1‑133 (MQ=255)                                              ttCTATGAAGAACTGGAACGCGGGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGt                                                                                      <  2:510435/140‑1 (MQ=255)                                                        aaCTGGAACGCGTGCAGAAAACTCATGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACgg                                                                              <  1:429041/138‑1 (MQ=255)                                                        aaCTGGAACGCGTGCAGAAAACTCATGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACgg                                                                              >  2:429041/1‑138 (MQ=255)                                                        aaCTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTa                                                                                 >  1:327529/1‑135 (MQ=255)                                                        aaCTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTa                                                                                 <  2:327529/135‑1 (MQ=255)                                                        aaCTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGc                                                                            >  1:52267/1‑140 (MQ=255)                                                        aaCTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGc                                                                            >  1:265938/1‑140 (MQ=255)                                                        aaCTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGc                                                                            <  2:52267/140‑1 (MQ=255)                                                        aaCTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGc                                                                            >  2:496700/1‑140 (MQ=255)                                                        aaCTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGc                                                                            >  2:263729/1‑140 (MQ=255)                                                        aaCTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGc                                                                            >  1:150531/1‑140 (MQ=255)                                                          cTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCgg                                                                          >  2:398266/1‑140 (MQ=255)                                                          cTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCgg                                                                          >  2:384872/1‑140 (MQ=255)                                                               aCGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGGAGAtt                                                                     >  2:651699/1‑140 (MQ=255)                                                                 gcgTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGGAGATTAc                                                                   <  1:429004/140‑1 (MQ=255)                                                                  cgTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGGAGATTAct                                                                  >  1:768114/1‑140 (MQ=255)                                                                    tGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGGAGATTActct                                                                >  2:356080/1‑140 (MQ=255)                                                                      cAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGGAGATTACTCTAc                                                              <  1:496700/140‑1 (MQ=255)                                                                           aaCTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGGAGATTACTCTACCGGtg                                                         <  1:269423/140‑1 (MQ=255)                                                                              tCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGTGGAGATTACTCTACCGGtgctg                                                      <  2:493249/140‑1 (MQ=255)                                                                               cACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGGAGATTACTCTACCGGtgctgc                                                     <  1:776423/140‑1 (MQ=255)                                                                                aCGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGGAGATTACTCTACCGGTGCTgcg                                                    >  1:395385/1‑140 (MQ=255)                                                                                 cGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGGAGATTACTCTACCGGtg                                                         >  2:769663/1‑134 (MQ=255)                                                                                 cGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGGAGATTACTCTACCGGtg                                                         <  1:769663/134‑1 (MQ=255)                                                                                    gCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGGAGATTACTCTACCGGTGCTGCGCAgg                                                <  2:643804/140‑1 (MQ=255)                                                                                     ccGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGGAGATTACTCTACCGGTGCTGCGCAGGt                                               >  2:102367/1‑140 (MQ=255)                                                                                               tGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGGAGATTACTCTACCGGTGCTGCGCAGGTGATCCTGCcg                                     >  2:2119/1‑140 (MQ=255)                                                                                                 cTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGGAGATTACTCTACCGGTGCTGCGCAGGTGATCCTGCcgcg                                   <  2:150531/140‑1 (MQ=255)                                                                                                  tGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGGAGATTACTCTACCGGTGCTGCGCAGGTGATCCTGCcgcgc                                  <  1:356080/140‑1 (MQ=255)                                                                                                       gCCAGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGGAGATTACTCTACCGGTGCTGCGCAGGTGATCCTGCCGCGCGTAGt                             <  1:651699/140‑1 (MQ=255)                                                                                                          aGTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGGAGATTACTCTACCGGTGCTGCGCAGGTGATCCTGCCGCGCGTAGTCgg                          <  1:2119/140‑1 (MQ=255)                                                                                                           gTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGGAGATTACTCTACCGGTGCTGCGCAGGTGATCCTGCCGCGCGTAGTCGGt                         >  2:260962/1‑140 (MQ=255)                                                                                                           gTGGTTGGCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTACCTTGCATGGTAATAGCGTTGGTACGGGCGGAGATTACTCTACCGGTGCTGCGCAGGTGATCCTGCCGCGCGTAGTCGGt                         >  2:851778/1‑140 (MQ=255)                                                                                                                  gCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGGAGATTACTCTACCGGTGCTGCGCAGGTGATCCTGCCGCGCGTAGTCGGTTCGATgg                  <  2:21461/140‑1 (MQ=255)                                                                                                                  gCAATCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGGAGATTACTCTACCGGTGCTGCGCAGGTGATCCTGCCGCGCGTAGTCGGTTCGATgg                  <  1:556957/140‑1 (MQ=255)                                                                                                                    aaTCGACGGGGATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGGAGATTACTCTACCGGTGCTGCGCAGGTGATCCTGCCGCGCGTAGTCGGTTCGATGGaa                <  2:395385/140‑1 (MQ=255)                                                                                                                           ggggATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGGAGATTACTCTACCGGTGCTGCGCAGGTGATCCTGCCGCGCGTAGTCGGTTCGATGGAAGTGTATg         <  2:838275/140‑1 (MQ=255)                                                                                                                             ggATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGGAGATTACTCTACCGGTGCTGCGCAGGTGATCCTGCCGCGCGTAGTCGGTTCGATGGAAGTGTATGaa       <  1:384872/140‑1 (MQ=255)                                                                                                                             ggATGTTCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGGAGATTACTCTACCGGTGCTGCGCAGGTGATCCTGCCGCGCGTAGTCGGTTCGATGGAAGTGTATGaa       <  1:398266/140‑1 (MQ=255)                                                                                                                                  ttCCATAGCGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGGAGATTACTCTACCGGTGCTGCGCAGGTGATCCTGCCGCGCGTAGTCGGTTCGATGGAAGTGTATGAACAGCa  >  1:119868/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                     GATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGGAGATTACTCTACCGGTGCTGCGCAGGTGATCCTGCCGCGCGTAGTCGGTTCGATGGAAGTGTATGAACAGCA  >  NC_000913/1059622‑1059890

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