Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,586,671	C→T	D31N (GAT→AAT)	hsdR ←	type I restriction enzyme EcoKI endonuclease component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,586,671	0	C	T	100.0%	230.5 / NA	65	D31N (GAT→AAT)	hsdR	type I restriction enzyme EcoKI endonuclease component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/37);  total (28/37)

TCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCC  >  NC_000913/4586535‑4586807                                                                                                                                         |                                                                                                                                         tcaACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGattatc                                                                                                                                       >  1:1136735/1‑140 (MQ=255) tcaACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGattatc                                                                                                                                       <  1:1023316/140‑1 (MQ=255) tcaACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGattatc                                                                                                                                       <  1:57058/140‑1 (MQ=255)  caacaaAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATcc                                                                                                                                      <  1:134962/140‑1 (MQ=255)     caaAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCggg                                                                                                                                   >  1:197890/1‑140 (MQ=255)          gCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGtt                                                                                                                              >  1:561659/1‑140 (MQ=255)                 tGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTcc                                                                                                                       <  1:492476/140‑1 (MQ=255)                       gTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGca                                                                                                                 <  1:697942/140‑1 (MQ=255)                           aCGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTA‑TTTCCGCCGCACAGGc                                                                                                            >  2:966215/1‑140 (MQ=255)                            cGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGc                                                                                                            <  1:861290/140‑1 (MQ=255)                             gCAGGAGATCGTGTTTATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCg                                                                                                           <  2:1127112/140‑1 (MQ=255)                             gCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTA‑TTTCCGCCGCACAGGCGa                                                                                                          <  1:966215/140‑1 (MQ=255)                                  agaTCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGc                                                                                                      <  2:561659/140‑1 (MQ=255)                                    aTCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCAt                                                                                                    <  2:121998/140‑1 (MQ=255)                                    aTCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCAt                                                                                                    <  2:1136735/140‑1 (MQ=255)                                         gtTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGc                                                                                                            <  1:479947/127‑1 (MQ=255)                                         gtTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGc                                                                                                            >  2:479947/1‑127 (MQ=255)                                           tGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGaa                                                                                             >  2:535186/1‑140 (MQ=255)                                           tGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGaa                                                                                             >  1:625780/1‑140 (MQ=255)                                             aTTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGt                                                                                           >  2:38366/1‑140 (MQ=255)                                                   acaAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGa                                                                                     <  1:1102515/140‑1 (MQ=255)                                                   acaAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGa                                                                                     <  2:410478/140‑1 (MQ=255)                                                     aaGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACg                                                                                   <  2:1129166/140‑1 (MQ=255)                                                     aaGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACg                                                                                   <  1:90688/140‑1 (MQ=255)                                                     aaGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACg                                                                                   <  2:864/140‑1 (MQ=255)                                                          gggATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCtt                                                                              >  2:340588/1‑140 (MQ=255)                                                              tGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGa                                                                                     >  2:942720/1‑129 (MQ=255)                                                              tGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGa                                                                                     <  1:942720/129‑1 (MQ=255)                                                                  gAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAAtt                                                                      <  2:3798/140‑1 (MQ=255)                                                                   aGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTc                                                                     >  2:937072/1‑140 (MQ=255)                                                                   aGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTc                                                                     >  1:542251/1‑140 (MQ=255)                                                                   aGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTc                                                                     >  1:377584/1‑140 (MQ=255)                                                                     tAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCaa                                                                   >  2:777304/1‑140 (MQ=255)                                                                          gACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTg                                                              >  2:379893/1‑140 (MQ=255)                                                                           aCCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTgg                                                             >  1:710450/1‑140 (MQ=255)                                                                             cAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGAt                                                           <  2:753732/140‑1 (MQ=255)                                                                             cAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGAt                                                           <  1:473052/140‑1 (MQ=255)                                                                             cAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAACTGGAt                                                           <  1:937072/140‑1 (MQ=255)                                                                               aGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGAttt                                                         >  2:843597/1‑140 (MQ=255)                                                                               aGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGAttt                                                         >  1:418413/1‑140 (MQ=255)                                                                               aGATGTTTCGCTGTGGCTTCGCCAAACATACGCAGTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACACCCTTCAGGAATTCAAAATTGGAttt                                                         >  2:1065817/1‑140 (MQ=255)                                                                                     ttCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATtcat                                                   <  1:777304/140‑1 (MQ=255)                                                                                     ttCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATtcat                                                   >  2:33263/1‑140 (MQ=255)                                                                                     ttCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATtcat                                                   >  1:550388/1‑140 (MQ=255)                                                                                       cGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATtcatca                                                 <  2:625780/140‑1 (MQ=255)                                                                                       cGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATtcatca                                                 <  1:18165/140‑1 (MQ=255)                                                                                         cTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCAtt                                               <  1:740713/140‑1 (MQ=255)                                                                                              gCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGttat                                          <  2:197890/140‑1 (MQ=255)                                                                                              gCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGttat                                          <  1:379893/140‑1 (MQ=255)                                                                                              gCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGttat                                          <  2:1144518/140‑1 (MQ=255)                                                                                                  cGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAAt                                      <  1:535186/140‑1 (MQ=255)                                                                                                     cAAACATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCa                                   <  2:418413/140‑1 (MQ=255)                                                                                                       aaCATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCAtt                                 <  1:33263/140‑1 (MQ=255)                                                                                                        aCATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGTATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTg                                <  1:193861/140‑1 (MQ=255)                                                                                                        aCATACGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTg                                <  1:340588/140‑1 (MQ=255)                                                                                                             cGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCtgtg                           >  1:408815/1‑140 (MQ=255)                                                                                                             cGCATTTTAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCtgtg                           <  2:377584/140‑1 (MQ=255)                                                                                                                 ttttAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCa                                   >  1:46536/1‑128 (MQ=255)                                                                                                                 ttttAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCa                                   <  2:46536/128‑1 (MQ=255)                                                                                                                  tttAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGc                      >  2:817124/1‑140 (MQ=255)                                                                                                                    tAATCAGCGTCGTGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCt                    >  1:673122/1‑140 (MQ=255)                                                                                                                           cgtcgtGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCaaa             >  1:773841/1‑140 (MQ=255)                                                                                                                            gtcgtGTTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAat            >  2:850461/1‑140 (MQ=255)                                                                                                                                 gtTGGGATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTa       <  2:940967/140‑1 (MQ=255)                                                                                                                                    gggATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAgg    <  2:773841/140‑1 (MQ=255)                                                                                                                                     ggATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGc   >  2:479171/1‑140 (MQ=255)                                                                                                                                     ggATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGc   >  2:920085/1‑140 (MQ=255)                                                                                                                                      gATTATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGcc  <  1:669438/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                         TCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCC  >  NC_000913/4586535‑4586807

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