Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,200,881	G→A	A51V (GCC→GTC)	xisE ←	e14 prophage, putative excisionase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,200,881	0	G	A	100.0%	229.6 / NA	72	A51V (GCC→GTC)	xisE	e14 prophage, putative excisionase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (38/34);  total (38/34)

CTACTATAAAGCGGATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTT  >  NC_000913/1200745‑1201010                                                                                                                                         |                                                                                                                                  ctactaTAAAGCGGATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACg                                                                                                                                 <  2:1126989/139‑1 (MQ=255) ctactaTAAAGCGGATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACg                                                                                                                                 <  1:198487/139‑1 (MQ=255)         aaGCGGATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACg                                                                                                                         <  2:32004/139‑1 (MQ=255)            cGGATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCAc                                                                                                                      >  2:96061/1‑139 (MQ=255)            cGGATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCAc                                                                                                                      >  2:563332/1‑139 (MQ=255)              gATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTc                                                                                                                    >  2:1048223/1‑139 (MQ=255)               atatAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCg                                                                                                                   >  1:10558/1‑139 (MQ=255)               atatAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCg                                                                                                                   >  1:627978/1‑139 (MQ=255)                 ataAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGcc                                                                                                                 >  1:1183014/1‑139 (MQ=255)                 ataAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGcc                                                                                                                 <  1:54001/139‑1 (MQ=255)                 ataAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGcc                                                                                                                 >  1:801527/1‑139 (MQ=255)                         gggACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCgg                                                                                                         >  1:428550/1‑139 (MQ=255)                         gggACTGAAACATTATTCTTAAGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGATGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCgg                                                                                                         >  2:53867/1‑139 (MQ=255)                          ggACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGc                                                                                                        >  1:1241017/1‑139 (MQ=255)                             cTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACg                                                                                                                         >  2:1071943/1‑119 (MQ=255)                             cTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACg                                                                                                                         <  1:1071943/119‑1 (MQ=255)                                          cTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCgg                                                                                        >  2:328462/1‑139 (MQ=255)                                             aCGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTG‑TTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGa                                                                                    >  2:786059/1‑139 (MQ=255)                                              cgcgGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGa                                                                                    >  1:585735/1‑139 (MQ=255)                                                cgGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCt                                                                                                      <  2:523010/119‑1 (MQ=255)                                                cgGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCt                                                                                                      >  1:523010/1‑119 (MQ=255)                                                cgGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACa                                                                                  >  1:144912/1‑139 (MQ=255)                                                  gTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACaaa                                                                                >  1:1216882/1‑139 (MQ=255)                                                   tCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACaaaa                                                                               >  1:1249339/1‑139 (MQ=255)                                                    cTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAaaat                                                                              <  1:96061/139‑1 (MQ=255)                                                     tGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGAc                                                                                   >  1:609341/1‑133 (MQ=255)                                                     tGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGAc                                                                                   <  2:609341/133‑1 (MQ=255)                                                        cTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTG‑TTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAAtt                                                                         <  1:786059/139‑1 (MQ=255)                                                         tGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAAtt                                                                         <  2:557605/139‑1 (MQ=255)                                                         tGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAAtt                                                                         <  1:328462/139‑1 (MQ=255)                                                           ccATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGAc                                                                                   >  2:202763/1‑127 (MQ=255)                                                           ccATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGAc                                                                                   <  1:202763/127‑1 (MQ=255)                                                           ccATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCt                                                                       >  2:933422/1‑139 (MQ=255)                                                             aTTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCAc                                                                                                                          >  2:1025790/1‑86 (MQ=255)                                                             aTTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCAc                                                                                                                          <  1:1025790/86‑1 (MQ=255)                                                             aTTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGc                                                                     >  2:1301416/1‑139 (MQ=255)                                                              ttACTCAAAATCCGTTGCAAGAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCTACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAGATTCTGcc                                                                    <  2:1216882/139‑2 (MQ=255)                                                              ttACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCt                                                                    <  2:1136815/139‑1 (MQ=255)                                                                 cTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAg                                                                 <  2:543232/139‑1 (MQ=255)                                                                   cAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCa                                                               <  2:428550/139‑1 (MQ=255)                                                                      aaTCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCgg                                                                                                         >  1:86023/1‑94 (MQ=255)                                                                      aaTCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCgg                                                                                                         <  2:86023/94‑1 (MQ=255)                                                                       aTCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAac                                                           >  2:887740/1‑139 (MQ=255)                                                                           gTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAacgacg                                                       >  1:99337/1‑139 (MQ=255)                                                                            ttGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGc                                                      >  1:903723/1‑139 (MQ=255)                                                                             tGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCa                                                     <  1:53867/139‑1 (MQ=255)                                                                              gCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCaa                                                    <  1:933422/139‑1 (MQ=255)                                                                                    gTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGa                                              <  2:627978/139‑1 (MQ=255)                                                                                    gTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGa                                              <  1:1081457/139‑1 (MQ=255)                                                                                    gTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGa                                              <  2:99337/139‑1 (MQ=255)                                                                                    gTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGa                                              <  1:1301416/139‑1 (MQ=255)                                                                                            tCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTg                                      >  1:63340/1‑139 (MQ=255)                                                                                             cATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTgg                                     <  2:1241017/139‑1 (MQ=255)                                                                                              aTTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTgga                                    >  1:1173509/1‑139 (MQ=255)                                                                                                        aTTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTc                          >  2:136949/1‑139 (MQ=255)                                                                                                                  ttACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTg                                                                                      <  1:383620/69‑1 (MQ=255)                                                                                                                  ttACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTg                                                                                      >  2:383620/1‑69 (MQ=255)                                                                                                                  ttACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTc                >  2:388925/1‑139 (MQ=255)                                                                                                                  ttACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAA‑CGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGAt                                >  2:846746/1‑122 (MQ=255)                                                                                                                  ttACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAA‑CGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGAt                                <  1:846746/122‑1 (MQ=255)                                                                                                                   tACCAACTCCCCAACTAACTCGACGTCCTCACGCa                                                                                                                       <  1:1091792/35‑1 (MQ=255)                                                                                                                   tACCAACTCCCCAACTAACTCGACGTCCTCACGCa                                                                                                                       >  2:1091792/1‑35 (MQ=255)                                                                                                                    aCCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTcag              >  2:1165984/1‑139 (MQ=255)                                                                                                                    aCCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTcag              <  2:801527/139‑1 (MQ=255)                                                                                                                    aCCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTcag              <  2:10558/139‑1 (MQ=255)                                                                                                                        aCTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCagcagc          <  2:144912/139‑1 (MQ=255)                                                                                                                          tCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCg                                                                                         >  2:745338/1‑58 (MQ=255)                                                                                                                          tCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCg                                                                                         <  1:745338/58‑1 (MQ=255)                                                                                                                          tCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAacgacg                                                       <  1:989885/92‑1 (MQ=255)                                                                                                                          tCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAacgacg                                                       >  2:989885/1‑92 (MQ=255)                                                                                                                           ccccAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCa       <  2:1173509/139‑1 (MQ=255)                                                                                                                                actaactCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCAttctt  >  1:43409/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                  CTACTATAAAGCGGATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTT  >  NC_000913/1200745‑1201010

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