Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,447,891	G→A	*1260* (TAG→TAA)	tamB →	translocation and assembly module for autotransporter export, inner membrane subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,447,891	0	G	A	100.0%	256.8 / NA	80	*1260* (TAG→TAA)	tamB	translocation and assembly module for autotransporter export, inner membrane subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (43/37);  total (43/37)

TGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAG  >  NC_000913/4447770‑4448020                                                                                                                          |                                                                                                                                  tGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCt                                                                                                                 >  1:46047/1‑139 (MQ=255)             aCTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTAc                                                                                                     >  2:84050/1‑139 (MQ=255)               tctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCGCTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTc                                                                                                                  <  2:918221/124‑1 (MQ=255)               tctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCGCTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTc                                                                                                                  >  1:918221/1‑124 (MQ=255)               tctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGc                                                                                                   >  1:976393/1‑139 (MQ=255)               tctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGc                                                                                                   <  1:487319/139‑1 (MQ=255)                 tataGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCca                                                                                                 >  1:1026496/1‑139 (MQ=255)                         cacTCACGTTACGTTATCGCCTGATACCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCTCTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACCAACAAg                                                                                         >  2:69949/1‑139 (MQ=255)                             cACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                     <  2:976393/139‑1 (MQ=255)                             cACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                     <  1:84050/139‑1 (MQ=255)                                 ttacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGt                                                                                 <  2:77569/139‑1 (MQ=255)                                  tacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTc                                                                                <  1:745913/139‑1 (MQ=255)                                   acgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCa                                                                               >  1:817594/1‑139 (MQ=255)                                       taTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTgg                                                                           <  1:513297/139‑1 (MQ=255)                                              tGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGTCCa                                                                    >  2:142341/1‑139 (MQ=255)                                                aTGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                     >  1:130979/1‑120 (MQ=255)                                                aTGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                     <  2:130979/120‑1 (MQ=255)                                                   ccTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGcc                                                               <  2:51827/139‑1 (MQ=255)                                                       aGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCa                                                                               <  1:1010234/119‑1 (MQ=255)                                                       aGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCa                                                                               >  2:1010234/1‑119 (MQ=255)                                                            tatCTGGAAGCCGTGTCTGGTGTAGACCTGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTg                                                      >  2:444250/1‑139 (MQ=255)                                                            tatCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACa                                                                                                                                 >  2:114786/1‑64 (MQ=255)                                                            tatCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACa                                                                                                                                 <  1:114786/64‑1 (MQ=255)                                                            tatCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                     >  2:1095514/1‑108 (MQ=255)                                                            tatCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                     <  1:1095514/108‑1 (MQ=255)                                                            tatCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTg                                                      >  2:295307/1‑139 (MQ=255)                                                               cTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGc                                                   <  2:1026496/139‑1 (MQ=255)                                                                 ggAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGa                                                 >  1:449130/1‑139 (MQ=255)                                                                  gAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGAt                                                                         <  2:551960/114‑1 (MQ=255)                                                                  gAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGAt                                                                         >  1:551960/1‑114 (MQ=255)                                                                   aaGCCGTGTCTGGTGTAGACCAGGCACTGTA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCaa                                                                   <  2:269149/119‑1 (MQ=255)                                                                   aaGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCaa                                                                   >  1:269149/1‑119 (MQ=255)                                                                   aaGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGAtt                                               >  2:731034/1‑139 (MQ=255)                                                                    aGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                     <  2:1098846/100‑1 (MQ=255)                                                                    aGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                     >  1:1098846/1‑100 (MQ=255)                                                                      ccGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCa                                            >  1:142891/1‑139 (MQ=255)                                                                         tgtCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTa                                         >  1:122206/1‑139 (MQ=255)                                                                          gtCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATAttt                                                                                                                     >  2:773976/1‑62 (MQ=255)                                                                          gtCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATAttt                                                                                                                     <  1:773976/62‑1 (MQ=255)                                                                          gtCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGTCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTAt                                        <  1:391942/139‑1 (MQ=255)                                                                               gtgtAGACCAGGCTCTGCA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCCGTTACTGGGCGATTTCAGTATCGATa                                   <  2:46047/139‑1 (MQ=255)                                                                                tgtAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATaa                                  <  1:444250/139‑1 (MQ=255)                                                                                  tAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGTCCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACt                                <  1:142341/139‑1 (MQ=255)                                                                                    gACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTAc                              <  1:731034/139‑1 (MQ=255)                                                                                     aCCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                     <  1:1113615/83‑1 (MQ=255)                                                                                     aCCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                     >  2:1113615/1‑83 (MQ=255)                                                                                      ccAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCa                            <  2:817594/139‑1 (MQ=255)                                                                                         ggCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAgg                                                                                        >  1:957460/1‑76 (MQ=255)                                                                                         ggCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAgg                                                                                        <  2:957460/76‑1 (MQ=255)                                                                                         ggCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTgg                                                                           >  1:229097/1‑89 (MQ=255)                                                                                         ggCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTgg                                                                           <  2:229097/89‑1 (MQ=255)                                                                                         ggCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGATCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGtt                         >  1:263868/1‑139 (MQ=255)                                                                                         ggCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCaa                                                                   <  2:416439/97‑1 (MQ=255)                                                                                         ggCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCaa                                                                   >  1:416439/1‑97 (MQ=255)                                                                                         ggCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGtt                         >  2:454091/1‑139 (MQ=255)                                                                                          gCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTAcc                             >  1:479732/1‑134 (MQ=255)                                                                                          gCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTAcc                             <  2:479732/134‑1 (MQ=255)                                                                                          gCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTg                        >  1:546319/1‑139 (MQ=255)                                                                                                 gattGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGcc                                                               >  1:357630/3‑94 (MQ=255)                                                                                                 gattGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGcc                                                               <  2:357630/92‑1 (MQ=255)                                                                                                     gCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                     <  1:1177/68‑1 (MQ=255)                                                                                                     gCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                     >  2:1177/1‑68 (MQ=255)                                                                                                      ctctATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGAt                                                                         >  2:719552/1‑79 (MQ=255)                                                                                                      ctctATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGAt                                                                         <  1:719552/79‑1 (MQ=255)                                                                                                         tATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCa          >  2:1018615/1‑139 (MQ=255)                                                                                                           tCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACt        >  1:1114232/1‑139 (MQ=255)                                                                                                                tCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGt                                          <  1:82423/100‑1 (MQ=255)                                                                                                                tCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGt                                          >  2:82423/1‑100 (MQ=255)                                                                                                                tCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTAc                              >  2:494469/1‑112 (MQ=255)                                                                                                                tCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTAc                              <  1:494469/112‑1 (MQ=255)                                                                                                                tCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATcc    <  2:153154/138‑1 (MQ=255)                                                                                                                tCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATcc    <  2:1033212/138‑1 (MQ=255)                                                                                                                tCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATcc    >  1:1033212/1‑138 (MQ=255)                                                                                                                tCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATcc    >  1:153154/1‑138 (MQ=255)                                                                                                                 cGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCaa                                                                   <  2:506155/74‑1 (MQ=255)                                                                                                                 cGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCaa                                                                   >  1:506155/1‑74 (MQ=255)                                                                                                                 cGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAg                                           >  1:258619/1‑98 (MQ=255)                                                                                                                 cGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAg                                           <  2:258619/98‑1 (MQ=255)                                                                                                                 cGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAg  >  2:122472/1‑139 (MQ=255)                                                                                                                 cGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAg  >  1:106500/1‑139 (MQ=255)                                                                                                                          |                                                                                                                                  TGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAG  >  NC_000913/4447770‑4448020

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